PHOTOGRAPH  THIS  SHEET 

,-N.  p  n  ' 

f  <  i 

\j  V  A 

UUt 

LEVEL 

INVENTORY 


DOCUMENT  IDENTIFICATION 

ftP#  !9<fO 


DISTRIBUTION  STATEMENT 


NT1S  GRAJkl 

DT1C  TRAC  *L] 

UNANNOUNCED  D 

JUSTIFICATION 


H  V 

DISTRIBUTION’ _ 

AVAILABILITY  CODES 

d;s7r:3ut:o.n  availability  and/or  special 


DATE  ACCESSIONED 


DISTRIBLTION  STAMP 


DATE  RETURNED 


DATE  RECEIVED  LN  DTIC 


REGISTERED  OR  CERTIFIED  NUMBER 


PHOTOGRAPH  THIS  SHEET  AND  RETURN  TO  DTIC-FDAC 


DOCUMENT  processing  sheet 


STOCK  IS  EXHAUSTED 


WRDC-TR-90-4015 


AD-A219  828 


ELECTROMAGNETIC  MACHINES  WHICH  UTILIZE 
MICROGEOMETRY  FIELD  STRUCTURES 


Dr.  S.C.  Jacobsen 
Dr.  John  E.  Wood 

University  of  Utah 

Center  of  Engineering  Design 

Salt  Lake  City,  UT  84112 


April  1990 


Final  Report  for  Period  August  1984  -  June  1987 


Approved  for  Public  Release;  Distribution  is  Unlimited 


MATERIALS  LABORATORY 

WRIGHT  RESEARCH  AND  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT- PATTERSON  AIR  FORCE  BASE,  OHIO  45433-6533 


NOTICE 


When  Covemsenc  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  In  connection  with  a  definitely  Government-related 
procurement,  the  United  States  Government  Incurs  no  responsibility  or  any 
obligation  whatsoever.  The  fact  that  the  government  may  have  formulated  or 
in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not 
to  be  regarded  by  implication,  or  otherwise  in  any  manner  construed,  as 
licensing  the  holder,  or  any  other  person  or  corporation;  or  as  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  Invention 
that  may  In  any  way  be  related  thereto. 


This  report  Is  releasable  to  the  National  Technical  Information  Service 
(NTIS).  At  NTIS,  It  will  be  available  to  the  general  public,  including 
foreign  nations. 


This  technical  report  has  b.’en  reviewed  and  is  approved  for  publica¬ 
tion. 


i 


,  '  ,.A>:  .  '  ••''Var  °)D 

- - - c* — ‘ -  - 

CHARLES  R.  WRIGHT,  II  T.  USAF 
Program  Manager 


FOR  THE  COMMANDER 


STEVEN  R.  LECi.AIR,  MAO,  USAF 
Acting  3ranch  Chief 
Manufacturing  Research 
Air  Force  Materials  Laboratory 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization  please 
notify  krqc/mi.tm  »  WPAFB,  OH  45433-  ftsaa  to  help  ua  maintain  a  current 
mailing  list. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obligations,  or  notice  on  a  specific 
document . 


UNCLASSIFIED 


:&  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


RKINGS 


:d  OEClASSiFiCAHON/ DOWNGRADING  schedule 

N/A 


ERFORMiNG  ORGANIZATION  REPORT  NUMBER^) 


3  DISTRIBUTION/ AVAILABILITY  OF  REPORT 

Approved  for  Public  Release;  Distribution  is 
Unlimited . 


5  MONITORING  ORGANIZATION  REPORT  NuMBERiS, 
WRDC-TR-90-40 1 5 


6a  NAME  OF  PERFORMING  ORGANIZATION 
University  of  Utah 


6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGAN, ZATiON 

(if  applicable)  Materials  Laboratory  (WRDC/MLIM) 

Wright  Research  and  Development  Center 


7b  ADDRESS  (Cry.  Stitt.  and$ZlP  Code) 


Wright-Pa^terson  AFE,  Oh  >,5';22-  £522 


Sc.  ADDRESS  I  City.  Sta  re.  ana  ZIP  Code) 

Office  of  Research  Administration 

309  Park  Building 

Salt  Lake  City,  UT  84112 


8a  NAME  0=  "JND'NG  SPONSORING  8b  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NuVBt® 

ORGANIZATION  (II  applicable) 

DARPA  F33615-84-C-5092 


sc  ADDRESS  (Cry.  Stare,  ana  ZIP  Cooel 

Attn:  DARPA/ISTO  (Dr.  R.  Rosenfeld) 
Arlington,  VA  2209-2308 


10.  SOURCE  OF  FUNDING  NUVBE®S 


PROGRAM 
ELEMENT  NO 

62712E 


PROJECT 

"AS* 

NO 

NO 

4587 

_ _ 

00 

AO»*  uN'T 
ACCESS'ON  NO 
01 


ii  -iT.E  (mciuae  Security  Classification) 

Electromagnetic  Machines  Which  Utilize  Microgeometry  Field  Structures 


2  Dr,.S°sfc.A“JacobSsen,  Dr.  John  E.  Wood,  Center  for  Engineering  Design 


'3s 

from 


,  /-,n  /o-7  ia  ^TE  O*  3£pORT  (Year  Month.  Day)  '5  RAGE  COuN' 

36/30/87  1990  Aprii  137 


7  cQSA-i  CODES _  18  SU8jECT  TERMS  ( Continue  on  reverse  if  necessary  ana  identify  by  block  number) 

FA.D  i  group  I  Sv-B-GROup  Micro  Machining,  Micro  Electro  Mechanical  Systems,  Robotics, 

_ j _ j _  Sensors,  Actuators  .  .ectric  Field  Devices,  FET  Sensors, 

1  i  Position  Sensing  ______ _ _ 

9  ABSTRACT  [Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  Microfield  Project  at  the  Center  for  Engineering  Design  was  originally  motivated  by 
a  quest  to  understand  how  to  improve  the  performance  of  actuators  in  demanding  applications 
such  as  robots  and  artificial  limbs.  It  was  determined  that  by  making  the  internal 
interactive  elements  of  the  system  (actuator)  numerous,  closely  packed,  and  very  small  (on 
the  order  of  one  to  ten  microns),  sizeable  force  outputs  and  power  densities  could  be 
produced,  with  improvements  in  impedance  characteristics  over  conventional  actuators.  A 
by-product  of  these  studies  was  the  generation  of  a  number  of  very  interesting  spin-off 
applications  for  microelectromechanical  systems  in  general,  including  mechanical  sensors, 
optical  devices,  actuators,  and  chemical  sensors.  Thus,  for  the  Microfield  Project, 
multiple  exploratory  investigations,  both  experimental  and  analytical,  were  undertaken 
aimed  at:  (1)  understanding  important  issues  necessary  for  the  systematic  design,  analysis, 
fabrication,  and  testing  of  micro  electro-mechanical  sensors  and  actuators,  and  (2)  estab¬ 
lishing  the  feasibility  of  various  concepts.  (Continued) 


20  DISTRIBUTION/  AVAILABILITY  OF  ABSTRACT 

□  UNCuASSiF'EDMJNUMrED  □  SAME  AS  RPT  □  DTiC  USERS 

21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 

22a  NAVE  O'  RESpONS'BlE  INDIVIDUAL 

Cnarles  R.  Wright,  1LT 

22b  TELEPHONE  (Include  Area  Cude)  1 

(313)  255-8786  j 

22c  OFFICE  SYMBOL 
WRDC/MLIM 

DO  FORM  1473,  8«  mar 


83  APR  edition  may  oe  used  until  e»n»usteo 
All  otner  editions  are  obsolete 


SECURITY  CLASSIFICATION  Qf  thiS  PAGE 

fU.l  G«*vrftn«Mt  PnntiM  ttl#— 

UNCLASSIFIED 


UNCLASSIFIED 


Block  J 9  (Cone  inued) 


Activities  reported  herein.  include:  (!)  the  design  and  fabrication  of  micron-scale' 
field  generating  structures,  including  conductors  and  st at lor.ary-charged  electrets, 
(2)  the  analysis  c f  ferces  produced  by  the  resulting  complex  electromagnetic  fields 
produced  between  structures,  (3'  the  development  of  micro  position  FEI-based  sensors, 
and  { **  the  implementation  of  ccmrreher.sive  systems  which  utilize  actuating  elements, 
sensors,  and  feedback  centre!  principles  ir.  order  to  manipulate  small  objects. 


As  a  prototype  system  employing  small  actuac 
constructed  and  evaluated  a  Micro: leld-based 
Cantilevered  Fiber  Optic  Serve  System  (SCOFS 


r,  optical  and  sensor  elements, 
system,  referred  to  as  a  Small 


we 


MicroFields  Final  Report 


v 


Table  of  contents 

1 .  Introduction  and  Summary  1 

1 . 1  Origin  of  the  Project  2 

1 .2  Discussion  of  Important  Issues  and  Initial  Activities  4 

1.2.1  Electromagnetic  Fields  4 

1.2.2  Scaling  5 

1.2.3  Materials  which  Generate  and/or  Modify  Fields  5 

1.2.4  Controlling  the  Movement  of  Objects  -  Stabilization  7 

1 .2.5  Sensing  the  Position  of  Small  Elements  8 

1.2.6  The  First  Comprehensive  Experiments  -  a  Small  Cantilevered 

Fiber  Optic  Servo  System  (SCOFSS)  8 

1 . 3  Summary  of  Microfield  Experimental  Results  9 

1.3.1  Comments  9 

1.3.2  Micro  Servo  Systems  -  9 

1.3. 2.1  Objectives  9 

1.3. 2. 2  System  Configuration  10 

1.3. 2. 3  Analysis  of  Force  on  Armature  10 

1 . 3 . 2 . 4  Field  V alidation  Experiments  1 0 

1.3. 2. 5  Feedback  Control  ll 

1 . 3 . 2 . 5 . 1  Approach :  Controller  Formulation  1 1 

1.3. 2. 5. 2  Simulations  11 

1.3. 2. 5. 3  Experimental  Results  ll 

1.3. 2. 6  AC  Stabilization  12 

1.3. 2. 6.1  Approach  12 

1.3. 2. 6.2  Simulations  12 

1.3. 2. 6. 3  Experimental  Results  10 

1.3.3  Micro  Position  Sensing  Systems  13 


MicroFields  Final  Report 


vl 


1.33.1  Objective  of  Micro  Sensors 

13 

1.3  3.2  Open -Gated  FET  Positions  Sensors 

13 

1 .3.3.3  Optical  Position  Sensor 

14 

1.3. 3. 4  Capacitive  Position  Sensor 

14 

1.3.3  5  Hall  Effect  Position 

15 

1.3.4 

Electrct  Experiments 

15 

1 .3.4. 1  Geometry,  Materials,  Fabrication 

15 

1.3. 4. 2  Charging  Methods 

15 

1.3. 4. 3  Charge  Densities  and  Stability 

16 

Resources 

17 

1.4.1 

Important  Areas  of  Development 

17 

1.4.2 

Project  Personnel 

17 

1.4.3 

Equipment  and  Facilities 

20 

2.  Basic  Research  Studies 

41 

2 . 1  Electrets  as  Field  Emitters 

41 

• 

2.1.1  Fabrication 

41 

2. 1.1.1  Sheets 

41 

2. 1.1.2  Filaments 

42 

• 

2. 1 .2  Compensated  vs.  Uncompensated 

43 

2.1.3  Patterns 

45 

2. 1 .4  Non-uniformities 

46 

2.2  FETs  as  Sensor  Elements 

2.2.1  Chem-FETs 

2.2.2  SG-Chips 

2.2.3  K-Chips 

2.2.4  Octa-FETs 


47 

47 

48 

49 


SO 


MicroFields  Final  Report 


vil 


Mathematical  Studies  and  Simulations 

51 

2.2.5. 1 

Theory  of  Operation 

51 

2.2. 5. 2 

Sensing  of  Long  Conductive  Emitters 

51 

2.2.5. 3 

Dielectric  Effects  on  Elec  tret  Lines 

52 

2. 2. 5. 4 

Characteristics  of  Annular  Sensing  Gates 

54 

2.3  Simulation  and  Design 

2.3.1  Conceptual  Framework 

2.3.2  The  Organization  of  the  Design  System 

2.3.3  Field  Emitter  Solutions 


2.3.3. 1  Fiber-Fiber  Interactions 

2. 3. 3. 2  Fiber-Plate  Interactions 

2. 3. 3. 3  Sensor  Interactions 

2. 3. 3. 4  Numerical  Integration  Issues 

2.3.4  The  Model  Definition  Programs 

2.3.5  The  Mapping  Display  Manager 


2 . 3 . 5 . 1  Two  Dimensional  Field  Maps 

2. 3. 5. 2  Feature  Extraction 

2. 3.5. 3  Sensor  Sensitivity  Displays 

2. 3. 5. 4  Three  Dimensional  Force  Field  Displays 

2.3.6  The  Dynamic  Display  Manager 

2.3.6. 1  What  are  Animated  Displays 

2. 3. 6. 2  Solid  Modelling 

2. 3. 6. 3  'Real-Time' Dynamics 

2. 3. 6. 4  Computational  Restrictions 

2. 3. 6. 5  Two-Degree-of-Freedom  Triad-Fiber  Systems 

2. 3.6. 6  Six-Degree-of-Freedom  Triad-Fiber  Systems 


55 

55 

55 

56 

56 

57 
57 

57 

58 
58 

58 

59 

60 
60 
61 
61 
62 
62 
63 
63 


65 


MicroFields  Final  Report 


viii 


2.3.7  The  Static-Display  Manager 

3.  A  Small  Cantilevered  Optical  Fiber  Servo  System:  SCOFSS  in  Detail 

3. 1  Introduction 

3.2  Experimental  Apparatus 

3.3  System  Forces 

3.3. 1  Passive  System  Forces 

3.3.2  Active  System  Forces 

3.3.3  Stability  Considerations 

3.3.4  Beam  Alignment 

3.4  Closed  Loop  Servo  Control 

3.4.1  State  Estimator 

3.4.2  Control  Law 

3.4.3  Driver 

3.4.4  Mechanical  System 

3.5  Servo  Operation 

3.5.1  Trajtc  lor,  Opera  Li 

3.5.2  Regulator  Operation 

3.6  Experiments 

3.6. 1  Linear  Charge  Density 

3.6.2  Control  Space 

3.6.3  Step  Input 

3.6.4  Sinusoidal  Input 

3.7  Higher  Modes 

3.7.1  Mode  Amplitudes 

3.7.2  Stability  Considerations 

3.7.3  Controller  Delays 


3.8  Ehscussion 


MicroFields  Final  Report 


ix 


• 

4. 

Oscillatory  Stabilization  of  a  Small  Cantilevered  Optical  Fiber 

95 

4. 1  Introduction 

95 

4.2  Oscillatory  Stabilization 

95 

• 

4.3  SCOFSS:  In  Oscillatory  Mode 

96 

4.4  Results 

98 

4.4.1  Computer  Simulations 

98 

• 

4.4.2  Experimental  Results 

99 

4.5  Multimode  Systems 

100 

4.6  Conclusions 

103 

• 

5. 

Actuators 

105 

5.1  Background 

105 

5.2  Forces  on  Conductors,  Dielectrics  and  Ferroelectric  Materials 

105 

• 

5.2.1  Force  Between  Parallel  Conductive  Plates 

105 

5.2.2  Forces  on  Dielectric  Materials 

107 

5.2.3  Forces  on  Ferroelectric  Materials 

109 

5 . 3  Forces  Generated  by  Multi -Element  Planar  Actuators 

110 

5.3.1  Introduction 

110 

5.3.2  Force  Analysis 

113 

• 

5.4  Fiber-Based  Actuators 

117 

5.4.1  Introduction 

117 

• 

5.4.2  Fiber-to-Fiber  Systems  (FTF) 

117 

BIBLIOGRAPHY 

127 

ATTACHMENTS 


148 


MicroFields  Final  Report 


1 


1  .  IN  r  .ODUCTION  AND  SUMMARY 

This  Final  Report  documents  work  done  by  the  Center  for  Engineering  Design  on  the  Contract 
"Electromagnetic  Machines  which  Utilize  Microgeometry  Field  Structures",  for  the  period  2d 
August  1984  through  30  June  1987.  This  report  continues  from  the  Final  Report  for  Phase  1. 
dated  15  December  1984.  That  Phase  (from  10  January  1983  through  8  June  1984)  covered 
preliminary  investigations  on  scaling,  actuator  designs,  filament-based  sensors,  disk -based  sensoi 
and  spatial  light  modulators,  proportional  and  multi-state  systems,  and  the  use  of  electrets  and 
FETs  as  interacting  sensor  emitter  and  detector  elements.  A  technical  summary  of  the  research  and. 
the  basic  results  of  the  Microfields  project,  covering  all  Phases  through  30  June  1987,  is  presented 
below.  Various  sections  of  this  report  have  been  extracted  from  papers  published  or  to  be 
published,  and  from  theses,  by  various  personnel  within  the  CED  (see  section  entitled' 
"Attachments",  at  the  end  of  this  report). 

Technical  Summary: 

The  Microfield  Project  (MF)  was  originally  motivated  by  a  quest  for  better  servo  actuators  lot  io. 
in  demanding  applications  such  as  robots  and  prosthetic  limbs.  As  a  result  of  initial  wore 
actuators,  a  number  of  additional,  very  interesting  applications  for  micro  elevtro-mc.  h.mi. 
systems  became  apparent.  In  particular,  it  was  determined  that  by  making  the  internal  mtera.  to.  c 
elements  of  the  system  (actuator)  numerous,  closely  packed,  and  very  small  ton  the  order  o!  o:  e  to 
ten  microns),  sizeable  force  outputs  and  power  densities  could  be  produced,  vv  ith  improv  cmcn>  c 
impedance  characteristics  over  conventional  actuators.  The  internal  charged  structure'  vvould 
consist  of  conductive  elements  and  fixed-charge  dielectrics  (electrets).  An  SEM  can  appropriated, 
pattern  the  fixed-charge  elements,  while  present  1C  technology  is  capable  of  making  the  cotid.x  to  . 
elements  on  the  appropriate  scale  and  geometry.  Moreover,  1C  technology  could  also  be  used  to 
fabricate  any  necessary  internal  sensing  elements  (to  determine  actuator  length)  and  also  to  generate 
electronic  logic  elements  (such  as  for  performing  commutation  switching),  thus  putting  all  clement' 
on  the  same  scale.  Thus,  for  the  Microfield  Project,  multiple  exploratory  investigations,  bom 
experimental  and  analytical,  were  undertaken  aimed  at:  (1)  understanding  important  issue-- 
necessary  for  the  systematic  design,  analysis,  fabrication,  and  testing  of  micro  electro-mechaniv.d 
sensors  and  actuators,  and  (2)  establishing  the  feasibility  of  various  concepts. 

Our  most  important  activities  have  included:  (1)  the  design  and  fabrication  of  micron-scale  field 
generating  structures,  including  conductors  and  stationary-charged  electrets,  (2)  the  analysis  <>: 
forces  produced  by  the  resulting  complex  electromagnetic  fields  produced  between  structures. » ;  i 
the  development  of  micro  position  FET-based  sensors,  and  (4)  the  implementation  oi 
comprehensive  systems  which  utilize  actuating  elements,  sensors,  and  feedback  control  principle* 
in  order  to  manipulate  small  objects. 

Technical  Results: 

Results  of  the  initial  investigations  provided  by  the  Microfield  Project  have  been  promising  It  o 
clear  that  small  scale  can  provide  improvements  in  many  systems  including  mechanical  sensors, 
optical  devices,  actuators,  and  chemical  sensers.  We  have  gained  a  better  understanding  of  settling 
issues  and  have  developed  a  library  of  analytical  procedures  capable  of  producing  approximate 
predictions  of  system  behavior.  We  have  produced  electrets  of  significant  intensity  in  Teflon  b> 
electron  beam  implantation.  We  have  measured  the  forces  produced  on  electrets  by  surrounding 
fields,  and  have  used  open-gated  Field  Effect  Transistors  to  determine  the  position  of  an  elec  tret 
based  on  the  local  field  strength.  These  have  been  difficult  undertakings  since  the  characteristic- 
size  and  detail  of  the  electrets  and  FETs  have  been  on  the  order  of  tens  of  microns. 

Finally,  in  an  effort  to  focus  the  activities  of  the  group  on  a  prototypical  system,  we  constructed 
and  evaluated  a  Microfield-based  system,  referred  to  as  a  Small  Cantilevered  Fiber  Optic  Ser\o 
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System  (SCOFSS).  The  system  includes  a  very'  small,  cantilevered  optical  fiber  covered  with  an 
electret  materiaJ  (polycarbonate).  The  optical  fiber,  which  is  situated  about  500  microns  above  ten 
200-micron  wide  conductor  driver  strips,  has  been  stabilized  and  manipulated  by  clean,  fields 
controlled  by  classical  feedback  techniques  and  by  the  use  of  appropriately  configured  ahernatinc 
fields. 


1 . 1  Origin  of  the  Project 

The  Microfield  Project  in  the  Center  for  Engineering  Design  began  with  attempts  to  understand 
how  to  improve  the  performance  of  actuators  in  demanding  applications  such  as  robots  and 
artificial  limbs.  After  a  decade  of  work,  studying  and  applying  high  performance  actuato:  v  it 
was  concluded  that,  short  of  fundamental  alterations  in  the  wav  actuators  operate,  all  that 
will  be  e.sentially  exists  now.  That  is,  without  a  different  approach  to  the  design  a 
actuators  (which  are  really  energy  transformation  systems)  maximum  performance  linin'-  ii. im¬ 
probably  been  reached.  Even  if  new  materials  are  utilized  in  existing  configurations  (hydra;:!:. . 
electric  and  pneumatic  actuators)  only  slight  (factors  of  two)  variations  in  performance  may  be 
expected. 

As  a  result  of  these  disappointing  discoveries,  radically  different  alternatives  were  explored.  In 
personnel  in  the  CED.  A  review'  was  made  of  the  many  seemingly  innovative  approach.es  m 
actuation  that  appeared  in  the  literature.  In  the  end/most  turned  out  to  be  only  peculiar 
concepts  which  were  unable  to  deal  successfully  w  ith  the  real  function  of  actuators  -  that  is, 
the  efficient  transformation  of  energy  important  variables  such  as  torque, 
speed,  and  impedance  under  dynamic  control. 

As  a  result  of  various  studies  on  actuation  systems,  it  became  apparent  that  radical 
improvements  in  the  performance  of  actuators  could  be  achieved  by  a  change  in  basic 
approach.  For  example,  in  the  case  of  electric  motors,  higher  interactive  field  strength  is  really 
the  only  way  by  which  higher  torque  can  be  generated  at  lower  speed  without  the  need  for 
increased  gear  reduction  (while  maintaining  size,  weight,  impedance,  and  efficiency 
constraints).  An  important  way  in  which  field  strengths  can  be  intensified  is  by  placing 
smaller,  electrically  active  elements  in  closer  proximity. 

In  the  case  of  hydraulic  and  pneumatic  actuation  systems,  more  attention  must  be  paid  to  fluid 
valving  methods,  the  characteristics  of  the  variable  volume  elements  (i.e.,  cylinder  volumes) 
and  how  proportional  operation  is  achieved  (flow  dividing,  pulse  width  modulation,  etc  ). 
These  factors  dominate  system  efficiency  and  dynamic  performance.  Again,  many  of  these 
issues  can  be  more  appropriately  addressed  by  the  use  of  multiple,  small  elements. 

Since  biological  actuators  (muscle)  display  such  impressive  capabilities  they  were  also  studied 
early  in  the  project.  It  was  discovered  that  muscle  uses  a  number  of  interesting  features  and 
elements,  exquisitely  arranged,  to  achieve  its  rather  impressive  performance. 

These  include: 

1 .  Small  elements. 

2 .  Large  numbers  of  comp  onents. 

3 .  Complex  geometry  in  structure  and  in  generated  fields. 

4 .  Distributed  sensors  and/or  local  control  of  processes. 

5 .  Direct  energy  conversion  -  chemical  to  mechanical  work. 

In  general,  muscle  exhibits  striking  performance,  being  able  to  produce:  (1)  contraction 
intensities  of  over  100  pounds  per  square  inch  of  muscle  cross  sectional  area,  and  (2) 
contraction  speeds  of  over  two  lengths  per  second.  Muscle  is  also  quite  responsive,  due  to  the 
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distributed  nature  of  its  electrochemical  control  system,  which  can  generate  maximum  tension 
within  less  than  20  milliseconds  of  an  activation  command.  In  contrast  to  man-madt  machines 
such  as  DC  motors,  commonly  used  in  robots,  muscle  is  least  power  consumptive  near  stall 
(DC  motors  are  most  consumptive  at  slow  speeds).  Moreover,  w  rut  the  need  lor  gearing, 
the  mechanical  output  impedance  of  muscle  is  suostantially  lower  than  mat  of  man -made 
actuators  if  other  characteristics  are  equivalent.  1 

As  a  result  of  our  muscle  studies  it  was  decided  to  pursue  the  design  of  actuation  systems 
which  utilize  large  numbers  of  small  elements  appropriately  combined  to  produce  force  Past 
efforts  on  the  original  Microfteld  Project  are  reviewed  in  this  document  in  terms  ol  tin- 
following  three  sections.  Section  1 .2  reviews  a  number  of  issues  important  to  the  development 
of  micro  electro-mechanical  systems.  These  include: 

1 .  Electromagnetic  fields. 

2.  Scaling  issues. 

3 .  Materials  which  generate  and/or  modify  fields. 

4.  Strategies  for  controlling  the  movement  of  small  objects. 

5.  Methods  to  sense  the  position  of  small  elements. 

6.  Comprehensive  experiments  aimed  at  illuminating  important  issues 

Section  1.3  introduces  selected  experimental  results  including  work  with: 

1 .  A  small  cantilevered  optical  fiber  servo  system. 

2.  Work  with  micro  position  sensors  -  field  effect  transistors.  Hall  effect 
devices,  and  lateral  effect  photodiodes. 

3.  Issues  important  to  the  fabrication  and  handling  of  electrets 


1  The  rather  impr.  ssive  performance  of  muscle  is  made  possible  via  its  use  of  large 
numbers  of  small  functional  elements  for  the  production  of  force.  The  use  of  over  10I!' 
cross-bridges  per  muscle  bundle  is  made  possible  by  the  fact  that  muscle  is  really  an  area 
machine.  The  hexagonal  arrangement  of  sliding  filaments  generate  large  areas  across 
which  many  small,  weak,  cross-bridges  combine  to  produce  significant  output  forces. 
(There  is  on  the  order  of  106  cm2  of  area  per  cm3  of  volume  generated  within  muscle 
tissue). 

We  very  briefly  sketch  here  the  nature  of  muscle  physiology:  Muscle  produces 
actuation  forces  via  cross-bridge  linkages  between  actin  and  myosin  fibers.  The  myosin 
cross-bridges  (CB)  exist  on  a  100  angstrom  size  scale.  One  presently  accepted  theory  is 
that  the  cross-bridges  electrostatically  seek  helically  arrayed  sites  on  actin  fibers  to  which 
they  covalently  bind.  Due  to  the  close  proximity  between  bridges  and  sites,  actuation 
involves  strong  field  forces. 

As  cross-bridges  sequence  through  their  biochemical  cycle,  utilizing  adenosine 
triphosphate  (ATP),  they  trade  chemical  free  energy  for  mechanical  work  at  upwards  of 
50%  efficiency.  These  cross-bridge  units,  whose  operation  is  controlled  by  chemical 
transmitters  are  locally  self  regulated  (commutated)  to  sequentially  repeat  their  operation  in 
concert  with  adjacent  elements.  Approximately  1014  cross-bridges  are  present  pei  cm2  of 
sarcomere  cross-section. 

The  actin  and  myosin  filaments  are  structured  such  that  individual  forces  add  in 
parallel  to  produce  muscle  contraction  intensities  of  approximately  100  pounds  per  in2. 
Sarcomeres  are  then  added  in  series  to:  (1)  determine  maximum  possible  contractile 
displacement,  (2)  increase  velocity,  and  (3)  decrease  muscle  stiffness. 
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Section  1.4  briefly  reviews  those  capabilities  necessary  to  work  in  the  Microfieids  area  These 
include: 

1 .  Technical  areas  where  important  understanding  must  be  generated 

2 .  The  acquisition  and  training  of  personnel. 

3.  Equipment  and  facilities. 

1.2  Discussion  of  Important  Issues  and  Initial  Activities  •  Scientific  and 
Engineering  Basis  for  the  Project 

1.2.1  Electromagnetic  Fields 

Electric  fields  with  millimeter  or  smaller  spatial  length  scales  can  be  considered  accurate!) 
as  electrostatic  fields  so  long  as  they  do  not  change  on  a  time  scale  shorter  titan  a 
nanosecond.  A  crucial  step  towards  using  these  fields  for  microsystems  is  to  understand 
the  forces  they  apply  on  objects.  This  might  appear  to  be  quite  straightforward  since 
electrostatics  is  presented  in  textbooks  as  a  long-ago  completed  exercise  in  applied 
mathematics.  But  the  real  world  is  not  so  straightforward. 

Real  devices  do  not  have  the  high  degree  of  symmetry  that  characterizes  idealized  tevtb  >ok 
problems.  Even  in  problems  limited  to  conducting  surfaces  and  known  static  charge 
distributions,  numerical  methods  must  be  used.  Thus,  both  analytical  and  numerK.il 
methods  for  studying  complex  field  structures  are  being  used  in  the  Microfield  Project, 
along  with  sophisticated  computer  graphics  techniques  to  help  permit  the  visualization  ol 
field  structures. 

The  real  problem,  however,  is  not  in  finding  a  numerical  solution  to  a  problem,  but  m 
arranging  a  model  to  be  a  sufficiently  accura.  •  representation  of  the  physical  system  being 
designed.  Useful  systems  typically  contain  a  much  richer  set  of  ingredients  than  simple 
charges  and  conductors.  Dielectrics,  semiconductor  materials,  and  active  elements  (signal 
flow  elements  such  as  transistors,  etc.)  are  all  important  determinants  of  electrostatic  fields 
in  situations  of  interest.  In  addition,  microfield  electromagnetic  systems  are  particularly 
difficult  to  model,  since  fields  are  influenced  substantially  by  nearby  environmental  objects 

The  complications  do  not  end  with  the  determination  of  field  configuration  Forces  on 
electrical  elements  are  more  complex  than  simple  Coulombic  coupling.  In  particular, 
uncharged  dielectric  materials  experience  other  forces,  such  as  the  dielectrophoretic  or  DEP 
force  when  the  material  is  immersed  in  a  nonuniform  electric  field.  Furthermore,  the 
distribution  of  charge  on  a  conductor's  surface  varies  under  the  influence  of  nearby  charges 
thus  leading  to  interactions  between  conducting  and  charged  elements  which  are  often, 
somewhat  simplistically,  called  image  forces.  Charges  near  dielectrics  similarly  induce 
polarization  and  lead  to  image-like  interactions.  A  closed-form,  analytical  solution  for 
forces  in  any  interesting  real-world  device  is  exceedingly  implausible.  Thus,  other 
approaches  are  being  utilized  to  facilitate  the  design  of  Micro  Electro  Mechanical  Systems 
(MEMS).  Approximate  solutions  to  approximate  models  are  being  used  to  understand 
behaviors  along  with  carefully  formulated  experimental  procedures.  Efforts  to  aid  the 
visualization  of  forces  produced  by  fields  are  continually  pursued.  The  use  of  computer 
graphics  for  visualizing  the  nature  of  dynamical  systems  is  a  significant  part  of  the 
microsystems  effort. 

An  example  of  methods  developed  for  the  examination  of  field  characteristics  is  the 
apparent  potential  surface  illustrated  in  Fig.  1.1.  The  slope  of  the  surface  indicates  the 
direction  of  the  force  on  a  cantilevered  fiber  system  described  in  Figs.  1.8  and  1.9. 
Apparent  potential  surfaces  can  be  constructed  for  any  two-dimensional  force  field. 
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whether  or  not  it  is  conservative.  Unfortunately  this  method  cannot,  in  general,  be  applied 
to  systems  with  higher  dimensionality. 

1.2.2.  Scaling 

Understanding  how  size  scale  affects  the  behavior  of  a  system  can  greatly  enhance  intuition 
and  thus  judgment  in  the  design  of  physical  systems.  Scaling  permits  the  knowledge 
gained  from  the  study  of  one  system  to  contribute  to  understanding  another.  A  number  of 
simple  scaling  arguments  can  be  used  to  indicate  tha.  the  performance  of  an  electro¬ 
mechanical  system  can  be  improved  by  reducing  its  size.  More  specifically,  if  conducting 
elements  are  maintained  at  a  fixed  potential  difference  V,  and  if  L  is  the  length  scale  for  the 
fields  (the  characteristic  distance  over  which  the  fields  change  significantly),  then  as  L 
decreases,  electric  field  strength  E  should  be  typically  V/L  and  hence  should  depend 
inversely  on  the  size  of  the  device,  i.e.,  as  (size)-l.  Scaling  also  affects  electrically 
chargeable  polymer  materials  which  can  retain  surface  charge  densities  typically  of  the 
order  of  30  nC/cm2  (nanocoulombs  per  square  centimeter).  The  charge  that  can  be 
deposited  on  an  object  is  then  proportional  to  its  surface  area  and  will  scale  as  (size)-. 

Since  an  object's  mass  scales  as  (size)3,  acceleration  in  a  fixed  E  field  scales  as  (size 1  and 
by  the  previous  argument  for  scaling  of  E,  in  a  field  produced  by  conductors  at  fixed 
potential,  achievable  acceleration  should  scale  as  (size)-2.  Each  type  of  interaction 
(Coulombic,  image,  and  DEP)  scale  somewhat  differently,  but  for  each  a  similar  analysis 
can  be  done. 


The  tendency  for  system  behavior  to  improve  with  decreasing  size  is  quite  general  and  is  an 
important  pan  of  the  motivation  for  the  Microfields  Project,  but  the  glib  argument  just 
presented  for  specific  scaling  laws  are  naive.  In  actual  situations,  voltage  differences 
cannot  be  held  constant  for  arbitrarily  small  size  scales.  Once  E  fields  become  too  large, 
dielectric  strengtn  will  be  exceeded,  polarization  of  field-generating  materials  can  degrade, 
etc. 

The  actual  scaling  of  performance  with  size  depends  strongly  on  a  number  of  issues. 
Thus,  we  ore  continuously  pursuing  rigorous  analyses  of  micro  systems  and  determining 
answers  to  some  difficult  questions  regarding  materials.  These  studies  should  produce  a 
more  accurate  capability  to  predict  the  performance  of  systems  as  their  size  is  reduced. 


1.2.3.  Materials  Which  Generate  and/or  Modify  Fields 

Materials  of  interest  in  connection  with  MEMS  include  conductors,  chargeable  polymers, 
dielectrics,  and  ferroelectrics.  Both  the  mechanical  and  electrical  properties  of  these 
materials  must  be  understood,  as  well  as  appropriate  fabrication  methodologies  if 
successful  systems  are  to  be  developed. 

The  central  component  of  many  proposed  micro  systems  is  the  electret,  which  consists  cf  a 
material  capable  of  maintaining  a  fixed  and  permanent  charge  distribution.  These  dielectric 
materials  can  emit  external  fields  via  two  mechanisms.  First,  they  can  be  internally 
polarized  (zero-net-charge),  and  second,  they  can  contain  implanted  charges  (non-zero- net- 
charge).  Of  course,  combinations  of  the  two  are  possible  as  shown  in  Fig.  1 ,2a. 

Electrets  are  of  interest  in  MEMS  primarily  because,  unlike  conductive  elements,  they  do 
not  need  to  be  connected  to  voltage  or  current  supplies  in  order  to  respond  to  local  fields 
Thus  they  may  be  used  as  free  floating  micro  armatures.  They  can  also  serve  as  field 
generators  capable  of  developing  forces  on  elements  equivalent  to  those  produced  by 
conductors  operating  at  high  voltages. 
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Using  appropriate  methods,  electrets  can  include  spatially  varying  charge  distributions  so 
that  nonuniform  (textured)  external  fields  can  be  produced.  The  class  of  materials  which 
seems  best  suited  to  this  application  are  polymers,  in  particular  Teflon®  FEP.  Studies 
report  that  negative  charge  densities  of  several  hundred  nanocoulombs  per  square 
centimeter  can  be  reliably  maintained  in  Teflon®  for  long  periods  of  time.  Electrets  also 
respond  differently  in  fields  than  conductors  since  charges  are  relatively  immobile  and  the 
electrets’  surface  is  not  an  equipotential. 

Studies  in  the  literature  almost  all  deal  with  compensated  charged  polymers.  These  are 
sheets  of  polymer  with  attached  conducting  surfaces  which  can  acquire  a  positive  charge 
which  neutralizes  the  polymer's  negative  charge  resulting  in  a  zero  net  charge  system. 
Compensated  electrets  emit  only  weak  external  fields  and  are  therefore  of  limited  interest 
for  micro  systems.  Unfortunately,  virtually  all  strongly  charged  electrets  described  in  the 
literature  are  compensated.  To  our  knowledge,  we  are  the  only  group  actively  pursuing  the 
fabrication  of  strong  uncompensated  electrets  with  spatially  varying  charge  densities.  As 
will  be  shown,  such  systems  will  have  a  wide  variety  of  application. 

An  example  of  an  actual  electret  is  shown  in  Fig,  1.2b  which  is  further  described  in  Section 
3.4.  In  this  case,  a  scanning  electron  microscope  beam  has  implanted  electrons  into  a 
Teflon®  FEP  sheet  of  25  microns  in  thickness.  The  field  is  visualized  via  sprinkled 
Xerox®  toner  particles  which  are  aggressively  attached  to  the  surface  due  to  the  field 
emitted  by  the  Teflon®.  Line  widths  on  Fig.  1.2b  are  approximately  20  microns  and  the 
height  of  the  smallest  letters  is  250  microns. 

A  number  of  questions  remain  regarding  materials  issues  which  will  be  pursued  in 
upcoming  work.  These  are: 

1.  What  materials  are  best  as  uncompensated  electrets9  Reported 
results  indicate  that  several  polymers.  Teflon®  TFE,  Teflon® 

FEP,  polycarbonate,  polystyrene,  and  others  are  effective  in 
trapping  electrons.  Our  experiments,  however,  have  been  spotty 
and  information  is  needed  on  how  to  process  these  materials  more 
efficiently. 

2.  What  are  the  best  charging  methods?  A  number  of  methods  have 
been  reported  in  the  literature.  Of  particular  interest  is  charging  via 
a  scanning  electron  microscope  beam  which  allows  considerable 
control  and  the  ability  to  "paint"  fine  textured  charge  patterns. 

3.  What  charge  densities  can  be  achieved?  Answers  are  needed  both 
for  surface  charge  density  on  electret  sheets  and  for  linear  charge 
density  on  electret  fibers.  These  answers  require  the  development 
of  instrumentation  for  measuring  charge  densities  and  other 
parameters. 

4.  How  stable  is  an  uncompensated  electret's  charge?  That  is,  do  the 
electrets  retain  their  charge  for  usefully  long  periods?  A  number 
of  papers  review  related  issues,  but  factors  which  affect  charge 
longevity  are  not  really  understood. 

5.  What  is  the  behavior  of  implanted  charge  as  length  scales  shrink? 

Studies  of  compensated  polymers  show  that  electrons  tend  to  be 
implanted  at  a  level  5-10  microns  deep.  What  happens  when 


MicroFields  Final  Report 


7 


polymer  thickness  is  less  than  5  microns?  The  answer  to  this 
question  will  govern  the  scaling  behavior  of  MEMS  using  charged 
polymers. 

6.  What  are  the  behaviors  of  other  interesting  materials?  Ferroelectric 
materials  placed  in  strong  electric  fields  can  become  strongly  and 
permanently  polarized,  i.e.,  they  can  acquire  a  very  high  density 
of  electric  dipole  moments.  In  principle  internal  polarization  can 
produce  enormous  surface  charge  densities  and  electric  fields  with 
very  high  strength.  For  example,  more  than  105  nC/cm?  and  10* 
volt/cm  can  be  produced  in  principle  with  lithium  niobate.  In 
practice,  effective  surface  charge  densities  and  external  fields  are 
generally  not  encountered  in  studies  of  ferroelectrics;  surface 
charges  are  masked  by  compensating  charges  in  attached 
conductive  layers.  Even  if  there  were  no  such  conductive 
field  strengths  of  10*  volt/cm  are  implausible.  Such  fields  would 
surely  cause  breakdown  and  surface  emission,  surface 
conduction,  etc.,  and  would  allow  for  compensation  via  surface 
charge  collection.  It  is  not  known  to  what  extent  this 
compensation  is  complete  and  whether  successfully  strong  and 
long  lasting  fields  may  remain. 

1.2.4,  Controlling  the  Movement  of  Objects  •  Stabilization 

One  of  the  most  important  tasks  of  the  Microfield  Project  (and  a  task  that  demands  some 
accuracy  in  analysis)  is  to  stabilize  elements  in  the  presence  of  fundamentally  unstable 
electrostatic  forces.  According  to  Eamshaw’s  theorem  (Maxwell,  1873),  a  charged  object 
in  an  electrostatic  field  produced  by  a  fixed  charge  distribution  cannot  be  in  stable 
equilibrium.  A  related  theorem  by  Jones  and  Bliss  (1977)  states  that  a  dielectric  object 
cannot  be  in  stable  equilibrium  (unless  the  object  is  suspended  in  a  fluid  of  higher  dielectric- 
constant)  under  the  action  of  static  DEP  forces.  These  two  fundamental  truths  do  not 
encompass  all  possible  electrostatic  elements,  but  as  a  practical  matter  electrostatic 
equilibrium  is  unstable. 

Our  efforts  now  include  exploration  of  three  methods  by  which  fields  can  be  augmented  to 
produce  stable  systems.  The  most  generally  applicable  of  these  is  the  control  of  fields 
through  sensing,  control  logic,  and  active  feedback;  a  classical  feedback  approach.  The 
general  scheme  is  illustrated  in  Fig.  1.3.  A  movable  object  with  electric  features  (an 
armature  which  may  be  charged  or  polarized)  moves  in  an  electric  field  created  by  field 
generators  in  its  base  (sets  of  charges,  dielectrics,  conducting  elements,  etc.).  Sensors  read 
a  signal  via  optical  or  electrical  means  from  which  a  state  estimator  infers  the  actual  state  of 
the  armature  (position,  velocity,  etc.).  The  controller  logic  then  determines  the  resultant 
force  necessary  to  move  the  armature  along  some  trajectory  toward  a  desired  state.  The 
controller  logic  includes  a  movement  strategy  as  well  as  a  model  of  the  dynamics  of  the 
system.  The  field  controller  computes  and  produces  the  voltages,  which  when  applied  to 
field  generating  elements  in  the  base,  generate  the  necessary  force.  Field  controller 
behavior  is  also  based  on  a  model  of  the  fields  and  the  interaction  of  fields  with  the  object. 

A  second  method  of  stabilizing  electrostatic  systems,  AC  stabilization,  applies  an 
oscillatory  electrical  field  to  an  object  which  is  otherwise  in  unstable  equilibrium.  For  a 
certain  range  of  field  strengths  and  frequencies,  the  oscillation  stabilizes  the  motion  of  the 
object.  Farther  discussion  of  AC  stabilization  is  contained  in  Sections  1 .3.2.6.  and  4.0. 


MicroFields  Final  Report 


8 


A  third  method  of  stabilizing  system  elements  is  to  add  non  electrical  forces  via  mechanical 
connection,  i.e.,  elastic  forces. 

1.2.5.  Sensing  the  Position  of  Small  Elements 

A  number  of  MEMS  will  require  the  use  of  micro  sensors  for  the  determination  of  system 
state.  Four  approaches  are  currently  being  explored  include:  (1)  the  use  of  field  effect 
transistors  together  with  electric  field  generators  to  determine  position,  (2)  optical  position 
sensors  including  lateral  effect  photodiodes  and  optical  occlusion  systems,  (3)  capacitance 
based  position  sensing,  and  (4)  Hall  effect  devices  together  with  magnetic  field  emitters  to 
sense  position. 

The  FET-based  system,  in  its  simplest  form,  consists  of  a  charged  object  which  moves  in 
the  proximity  of  a  field  effect  transistor  gate  (see  Fig.  1.4).  Our  work  initially  utilized 
chemical  sensitive  field  effect  transistors  for  the  experiments  shown  in  Fig.  1.5.  In  this 
case,  the  gate  is  moved  relative  to  a  charged  Teflon®  fiber  to  produce  quite  reliable  position 
measurements  shown  in  Fig.  1.6.  We  recently  produced,  in  the  HEDCO  laboratory  at  the 
University  of  Utah,  field  effect  transistors  more  appropriate  for  this  application,  as  shown 
in  Fig.  1.7. 

1.2.6.  The  First  Comprehensive  Experiment  -  The  Very  Small  Servo  System 

Scaling  arguments  suggest  that  the  realm  in  which  electrostatic  forces  are  most  interesting 
(i.e.,  when  electric  forces  dominate  over  inertial  forces)  is  that  of  micro  chip  size  scales 
(i.e.,  below  50  microns).  However,  procedures  with  such  small  systems  at  this  early  stage 
in  our  work  had  obvious  disadvantages.  Instead,  it  was  decided  to  perform  preliminary 
experimental  studies  on  small  but  "laboratop'  size"  apparatus  suitable  to  study  small  scale 
fields,  and  which  had  the  possibility  of  rapid  fabrication  and  manual  modification  and/or 
adjustment.  The  first  system  is  shown  schematically  in  Fig.  1.8  and  more  realistically  in 
Fig.  1.9.  The  actual  experimental  apparatus  is  shown  in  Fig.  1.10.  An  additional  system 
with  only  one  degree-of-freedom,  used  for  characterization  of  fibers,  is  shown  in  Fig. 
1.11. 

The  system  studied  consists  of  a  number  of  components  including  a  stator  or  base  which  is 
a  system  of  strip  drivers  fabricated  in  the  HEDCO  laboratory  (see  Fig.  1.12).  The  ten 
adjacent  planar  insulated  conducting  strips,  each  200  microns  by  2  cm.,  are  attached  to 
amplifiers  which  set  or  sense  strip  voltages  individually.  A  long,  thin,  movable  fiber  is 
introduced  into  the  region  of  the  strips  as  shown  in  Fig.  1.13.  Due  to  the  100  to  1  aspect 
ratio  of  the  strips  and  the  armature,  the  system  may  be  approximated  by  a  two  dimensional 
system  with  the  degrees-of-fireedom  x  and  y. 

The  basic  strip  driver  and  fiber  system  can  be  arranged  to  include  almost  all  of  the  elements 
of  a  total  microsystem  including  an  armature,  a  stator,  electrets,  dynamic  field  generators, 
field  sensors,  and  position  sensors  of  several  types.  The  system  was  first  used  for  field 
validation  as  shown  in  Fig.  1.14.  A  thin  wire  was  used  as  the  probe  (armature)  and  was 
connected  to  a  very  special,  high  input  impedance,  intermittent  sensing  voltmeter.  In  this 
way  the  field,  generated  by  strip  drivers,  in  a  number  of  voltage  configurations,  could  be 
measured  and  compared  to  mathematical  models  as  shown  in  Fig.  1.15. 

A  number  of  experiments  were  performed  with  a  cantilevered  armature  consisting  of  a  140 
micron  diameter  optical  fiber  as  shown  in  Fig.  1.16a-c.  These  optical  fibers  were  coated  or 
wrapped  with  electret  materials  as  shown  in  Figure  1.1 6d  (20  micron  Teflon®  wrapped  on 
the  fiber).  The  fibers  were  then  electrically  charged  to  a  linear  charge  density  on  the  order 
of  0.1  nCVcm  (nanocoulombs  per  centimeter). 
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Under  the  action  of  fields  produced  by  the  driver  strips,  the  fibers  deflect  and  the  deflection 
is  detected  with  a  lateral  effect  photo  diode  (LEPD)  monitoring  a  laser  beam  passing 
through  the  fiber.  This  basic  configuration  has  been  used  for  several  related  studies 
including: 

1.  Static  deflection  of  a  fiber  was  measured  for  various  strip  driver 
voltages,  with  the  results  used  to  confirm  measurements  of  linear 
charge  density  and  field  strength. 

2.  The  apparatus  was  used  to  test  an  active  feedback  system  for 
controlling  fiber  motion  and  stabilizing  it  in  otherwise  unstable 
positions  (positions  close  to  the  drivers  where  image  forces 
dominate).  Fig.  1.17  shows  motion  of  the  fiber  under  servo  control 
being  commanded  to  generate  a  five-pointed  star  and  the  letters 
DARPA. 

3.  The  apparatus  was  used  to  test  various  concepts  in  AC  stabilization. 

Fig.  1.19  illustrates  the  trajectory  of  a  cantilevered  fiber  in  an  AC 
field  for  stable  and  unstable  cases. 


1.3  Summary  of  Experimental  Results 
1.3.1  Comments 

The  earliest  work  on  the  Microfield  project  involved  literature  searching  and  analytic 
modeling  of  idealized  configurations.  Real  progress  was  limited  by  the  fact  that  many 
important  experimental  questions  were  not  addressed  in  the  literature,  and  that  analysis  was 
futile  without  a  clearer  picture  of  what  constraints  would  be  imposed  by  practical 
configurations.  Initial  experimental  efforts  were  therefore  chosen  to  begin  answering  the 
many  outstanding  questions,  and  to  foster  progress  on  a  broad  front.  Three  areas  of 
experimentation  were  developed:  (1)  a  microservo  system  as  a  testbed  for  very  many  of  the 
microsystem  ideas,  (2)  microposition  sensing,  and  (3)  experiments  with  electret  materials. 
These  three  areas  were  not  totally  disjoint,  e.g.,  the  microservo  system  used  microposition 
sensing  and  electrets,  the  electret  studies  used  some  measuring  techniques  developed  with 
the  micro  servo  system,  the  microposition  studies  used  electrets,  etc.  Nevertheless,  it  is 
convenient  to  separate  the  three  areas  in  the  descriptions  that  follow. 


1.3.2  Micro  Servo  System  -  Two  Degree-of-Freedom  Cantilevered  Fiber 

The  CED  has  pursued  the  development  of  a  two  degree-of-freedom  cantilevered  fiber 
system  which  allows  integrated  examination  and  evaluation  of  features  central  to  our 
proposed  microdevices:  model-based  servo  control,  position  sensing,  and  field  generation 
(actuation).  This  system  serves  as  an  intermediate  testbed  prior  to  making  systems  which 
are  smaller,  and  which  contain  more  degrees  of  freedom. 

1.3. 2.1  Objectives 

This  project  allows  for  the  comparison  of  analytical  and  experimental  results  of  field 
effects,  dynamic  effects,  control  system  performance  and  sensing  system  sensitivity 
and  accuracy,  and  provides  experience  with  microsystem  materials  (especially  with 
electrets). 
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1.3. 2. 2  System  Configuration 

At  present,  the  experimental  system  consists  of  a  cantilevered  quartz  fiber-  optic 
armature,  which  is  free  to  move  in  two  dimensions  (see  description  in  Section  1 .2.6). 
This  fiber  is  coated  with  a  charged  electret  material.  The  charged  region  is  acted  upon 
by  the  fields  emitted  from  the  base  structure,  which  consists  of  parallel  silicon 
(conductive)  strips  (of  about  200  microns  in  width)  which  may  be  independently 
connected  to  voltage  sources  or  sensors.  Voltages  applied  to  the  strips  through  driver 
amplifiers  are  used  to  create  desired  control  fields  and  control  forces. 

The  position  of  the  fiber  is  determined  by  projecting  a  laser  beam  through  the  fiber  optic 
onto  a  lateral  effect  photodiode  (LEPD),  which  then  outputs  a  reading  of  the  (y,z ) 
position  of  the  end  of  the  fiber  optic.  Alternate  studies  of  position  sensing  have  been 
made  using  the  driving  strips  as  sensing  strips.  This  is  done  by  monitoring  the  amour, . 
of  charge  which  is  induced  in  the  strips  by  the  presence  of  the  electreted  fiber. 
However,  these  are  intermediate  methods  of  position  sensing.  The  method  which  is 
most  promising,  and  for  which  we  have  an  ongoing  development  effort,  is  the  use  of 
open-gated  FETs.  FETs  designed  for  this  purpose  will  be  incorporated  into  our  next 
microchip  test  base. 

1.3.2.3  Analysis  of  Force  on  Armature 

The  forces  acting  on  the  armature  are  numerous  and  complex  (see  Fig.  1.18).  These 
include  the  gravitational  force,  the  image  force  between  the  charged  fiber  and  the 
conductive  base  strips,  the  elastic  restoring  force  (if  the  filament  is  cantilevered),  and 
the  control  forces  applied  by  the  nonzero  strip  potentials.  The  control  forces  are  further 
divided  into  monopolar  forces  (which  are  proportional  to  the  embedded  charge)  and  the 
dielectrophoretic  (DEP)  forces  (which  are  dependent  on  the  dielectric  nature  of  the  fiber 
and  the  gradient  of  the  surrounding  electric  fields).  For  the  present  system,  which  uses 
a  140  micron  diameter  fiber  optic,  the  DEP  forces  are  small  enough  to  be  ignored,  as  a 
reasonable  approximation.  (However,  as  devices  get  smaller,  the  DEP  forces  will 
become  more  dominant.)  All  of  the  other  forces  listed  above  are  included  in  the  model- 
based  controller  formulation  used  in  feedback  control. 

1.3.2.4  Field  Validation  Experiments 

In  order  to  evaluate  the  validity  of  the  mathematical  model  of  the  fields  generated  by  the 
control  strips,  experiments  using  a  potential  probe  were  performed.  The  probe 
consisted  of  a  fine  wire  (25  micron  diameter)  set  parallel  to  the  control  strips  (see  Fig. 
1.14).  The  potential  of  the  probe  was  then  monitored  at  a  variety  of  (y,z)  points;  the 
very  low  (picofarad  range)  probe  capacitance  required  special  measurement  techniques 
(computer-triggered  intermittent  readings,  driven  shields,  ultra  high  input  impedance 
electronics).  Measurements  were  compared  with  the  simple  mathematical  model  of  the 
field  which  assumed  strips  were  infinite  in  the  x-direction  and  that  the  z=0  plane  was 
grounded  everywhere  except  for  the  driven  strips. 

Measured  results  were  generally  in  good  agreement  with  the  theory  (Fig.  1.15)  and 
they  helped  to  evaluate  the  extent  to  which  experimental  techniques  and  analytical 
approximations  introduced  errors.  In  particular,  because  its  diameter  was  relatively 
large  (to  allow  mechanical  stability),  the  probe  could  not  resolve  fine  scale  features  of 
fields.  More  importantly,  except  fairly  close  to  the  strip  drivers,  measurements 
revealed  that  nearby  conducting  surfaces  significantly  modified  the  field,  and  that  the 
simple  mathematical  model  was  not  fully  adequate. 
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1.3. 2. 5  Feedback  Control 

1 .3.25.1  Approach:  Controller  Formulation 

The  feedback  approach  to  control  of  the  fiber  is  an  example  of  the  general  feedback 
scheme  outlined  in  Section  1.2.4  and  in  Fig.  1.3.  The  three  major  components  of 
the  control  system  are:  (1)  the  slate  estimator  which  utilizes  information  from 
sensors  to  determine  the  position  and  velocity  of  the  armature;  (2)  the  control  block 
which  defines  the  force  to  be  applied  to  the  armature,  depending  on  particular 
strategies;  and  (3)  the  stator  strip  driver  which  determines  base  strip  potentials 
required  to  impress  desired  forces  on  the  armature.  The  actual  system  also  includes 
the  electronics  of  driving  amplifiers,  sensing  amplifiers,  etc. 

The  total  force  on  the  armature,  due  to  passive  (gravity,  image  and  elastic)  and 
active  (control-related)  sources,  is  chosen  by  the  control  block  to  match  the  desired 
force  of  the  controller  strategy  and  is  a  function  of  positional  error  with  respect  to 
some  set  point,  and  of  the  velocity  of  the  fiber.  The  actual  acceleration  of  the 
armature  (fiber)  follows  from  the  combined  effect  of  all  of  the  above  forces  and 
should  follow  the  desired  acceleration  since  the  controller  chooses  base  strip 
potentials  which  compensate  for  the  passive  forces. 

At  present,  to  achieve  good  control  it  is  necessary  to  know  the  linear  charge  density 
of  the  fiber  reasonably  accurately.  In  future  systems,  adaptive  control  algorithms 
might  be  derived  which  compensate  for  such  parameter  uncertainties.  It  is  expected 
that  the  above  control  principles  will  be  applicable  to  systems  of  smaller  scale  and 
more  degrees  of  freedom. 

13.25.2  Simulations 

Numerous  simulations  of  the  cantilevered  fiber  system  have  been  conducted  in  an 
attempt  to  test  a  priori  the  controller  formulations  and  expected  system  behavior. 
The  first  simulations  used  only  two  base  driver  strips,  with  the  armature  assumed  to 
be  an  infinitely  thin  electret  filament,  and  the  base  constructed  of  a  dielectric  material 
which  influenced  the  armature  through  image  forces.  An  error-based  controller, 
with  velocity  damping,  was  used.  The  simulation  also  assumed  that  the  controller 
would  feed  driver  amplifiers  subject  to  saturation.  As  seen  in  the  controller 
apparent  potential  surface  of  Fig.  1.1,  an  equilibrium  point  exists  at  the  bottom  of 
the  surface.  When  the  fiber  deviates  from  that  point,  it  experiences  restoring  forces 
in  the  direction  of  the  gradient  of  the  surface. 

More  recent  simulations  have  employed  multiple  (more  than  two)  driver  strips.  In 
anticipation  of  saturation  limitations,  and  to  resolve  redundancies  (such  as  three 
applied  potentials  and  two  controlled  degrees-of-freedom),  the  sum  of  the  square  of 
the  plate  potentials  was  minimized.  It  is  this  procedure  which  is  presently 
implemented  in  the  experimental  setup.  For  sample  results  of  simulated  trajectories, 
see  Fig.  1.19. 

1.3.253  Experimental  Results 

A  variety  of  experimental  tests  have  been  conducted  in  order  to  validate  the  behavior 
of  the  model-hased  controller.  The  first  set  of  tests  were  conducted  without 
feedback,  in  fiber  positions  which  were  inherently  stable  under  combined  electric, 
gravitational,  and  elastic  forces.  The  driver  potentials  were  set,  and  the  static 
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displacement  of  the  fiber  was  then  monitored.  For  improved  signal-to-noise  ratio 
and  to  avoid  drift  in  the  electronics,  this  "static"  procedure  was  usually  earned  out 
at  very  low  frequency  (1  Hz).  The  displacement  results  agreed  very  well  with 
theory. 

A  second  set  of  tests  used  feedback,  with  a  commanded  set-point.  The  controller 
selected  the  strip  voltages  necessary  to  hold  the  fiber  statically  at  that  set-point.  The 
observed  positional  error  between  desired  and  actual  positions  was  on  the  order  of 
2%. 

A  number  of  dynamic  tests  have  been  conducted.  Although  the  controller  seems  to 
be  performing  as  in  its  computer  simulations,  a  number  of  technical  details 
introduced  errors.  For  instance,  the  system  has  a  low  intrinsic  damping.  Adding 
damping  by  increasing  the  velocity  gain  can  reduce  the  oscillations  around  a  set 
point,  however,  the  noise  of  the  velocity  signal  adds  noise  to  the  force  signal, 
which  leads  to  instability  at  the  higher  velocity  gains  required  for  fast,  crisp 
response.  Also,  higher  mode  oscillations  introduce  errors.  Despite  the  above 
problems,  a  number  of  tests  have  proved  promising.  These  include  the  generation 
of  a  five-pointed  star  (Fig.  1.17a),  the  writing  of  the  letters  "DARPA"  (Fig.  1  17b  ), 
the  tracking  of  a  sine  wave,  the  response  to  a  square  wave,  and  the  drawing  of  a 
box  (in  which  the  horizontal  motion  is  correct,  with  the  vertical  motion  showing 
some  error,  possibly  due  to  driver  saturation). 

1.3. 2.6  AC  Stabilization 


1.32.6.1  Approach 

An  alternate  method  for  controlling  the  stable  location  of  a  microdevice  might  be 
with  AC  stabilization,  which  does  not  require  a  model-based  controller  or  sensors 
to  feed  back  state  information.  Instead,  AC  stabilization  produces  a  local  stability 
by  adding  a  sinusoidal  force  to  the  fiber.  The  sinusoidal  force  field  must  be  formed 
so  that  it  does  not  change  the  position  of  the  fibe:  equilibrium  point,  but  simply 
serves  to  "stiffen"  (i.e.,  stabilize)  the  equilibrium.  This  has  required  a  voltage 
configuration  in  which  one  strip  is  used  to  create  a  steady  state  equilibrium  (gravity, 
elastic  and  electrostatic  forces)  and  two  symmetric  pairs  of  strips  are  driven  with 
AC  signals,  with  voltage  ratios  chosen  so  that  the  AC  electric  field  vanishes  at  the 
equilibrium  position  (see  Fig.  1.20). 

1. 3.2.62  Simulations 


For  small  oscillations  of  the  fiber,  the  equations  for  fiber  motion  can  be  linearized  in 
the  displacement  from  equilibrium.  In  the  linearized  approximation,  the  vertical  and 
horizontal  motions  act  independently  and  for  either  degree  of  freedom  the  stability 
can  be  analyzed  from  the  theory  of  Mathieu  equations.  It  was  not  clear  at  the  outset 
that  the  results  of  such  an  analysis  would  be  applicable  to  realistic  experimental 
situations  in  which  deviations  from  equilibrium  may  be  large.  For  this  reason,  a 
simulation  of  the  action  of  the  fiber  was  done  according  to  the  full  (nonlinear)  set  of 
dynamical  equations  (see  Fig.  1.21).  Generally  good  agreement  was  found  with 
Mathieu  equation  predictions  of  stability. 


1.3.2.63  Experimental  Results 

Demonstrations  of  AC-stabilizing  effects  have  progressed  despite  difficulties  with 
poor  signal-to-noise  ratios  caused  by  acoustical  oscillations  in  the  experimental 
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setup.  The  following  AC-stabilization  phenomena  have  been  observed:  (1)  a 
"stiffening"  (i.e.,  increase  in  frequency)  of  the  system's  fundamental  mode  when 
the  AC-stabilizing  signal  was  added,  as  shown  in  Fig.  1.22,  and  (2)  a  greatly 
decreased  sensitivity  of  a  marginally  stable  fiber  to  perturbations.  The  sudden 
transition,  as  fiber  height  is  decreased,  from  an  inherently  stable  fiber  configuration 
to  an  inherently  very  unstable  one,  has  so  far  prevented  the  demonstration  of  AC- 
stabilization  in  an  inherendy  unstable  position.  However,  experimental  technique, 
especially  with  improving  the  signal-to-noise  ratio,  is  evolving. 

1.3.3  Micro  Position  Sensing  Systems 

1.3.3. 1  Objective  of  Micro  Sensors 

A  means  to  determine  the  position  of  a  dynamic  element  bearing  a  pattern  of  charges  is 
a  central  requirement  for  most  of  the  applications  for  micro  systems.  In  the 
experiments  with  the  two  degree-of-freedom  cantilevered  fiber,  an  external  lateral  effect 
photodiode  has  been  used,  for  convenience,  to  monitor  the  fiber  position.  Future  micro 
systems  will  require  position  sensors  which  can  be  fabricated  in  the  base  structure  and 
which  provide  precise  position  information  for  feedback  control  systems  or  for 
generation  of  the  useful  output  signal  from  the  device.  Several  methods  of  position 
sensing  which  are  compatible  with  semiconductor  processing  have  been  considered: 

1 .  Open-gated  FETs  in  arrays  appropriate  for  one  to  six  degrees  of 
freedom. 

2.  Optically  active  arrays  in  the  base  responding  to  light  from  internal 
or  external  sources. 

3 .  Capacitive  sensors  measuring  potential  or  induced  charge. 

4.  Hall  effect  devices  in  a  system  with  magnetic  thin  films. 

The  area  of  application  of  the  different  methods  will  probably  depend  on  the  particular 
system.  For  example,  a  system  incorporating  a  source  of  light  emission  from  a  light 
emitting  diode  or  a  solid  state  laser  in  the  base  would  be  able  to  utilize  light  interactions 
more  conveniently  than  a  system  not  already  involving  light.  The  experiments 
conducted  with  open-gated  FETs  have  been  very  encouraging  in  that  the  response  has 
been  strong  enough  to  be  easily  measured  and  a  relatively  simple  model  has  agreed  well 
with  the  data  (Fig.  1.6). 

1.3.3.2  Open-Gated  Field  EfTect  Transistor  Position  Sensors 

The  sensitivity  of  the  current  to  electric  fields  at  the  gate  of  an  open-gated  FET 
suggested  that  a  FET  would  make  a  good  microfield  meter  (see  Figs.  1.4  -  1.7).  In 
early  experiments,  FETs  manufactured  in  the  HEDCO  laboratory  were  used  to  scan 
lines  of  electrons  implanted  by  a  scanning  electron  microscope  in  a  Teflon®  sheet.  It 
was  found  that  the  FET  current  was  stable,  relatively  free  from  noise,  and  could  be  read 
with  0.1  microamp  resolution.  The  current,  furthermore,  was  in  very  good  agreement 
with  a  simple  model  of  FET  behavior.  In  further  experiments,  small  SEM-charged 
spots  on  Teflon®  were  located  with  FETs,  although  the  FET  response  was  weak. 

These  experimental  results  were  the  basis  for  designing,  prior  to  fabrication  at  the 
HEDCO  laboratory,  a  FET  optimized  to  do  microposition  sensing  in  conjunction  with  a 
fine  textured  field  generating  electret.  From  the  experimental  results  it  has  been 
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calculated  that  this  device  will  perform  microposition  sensing  with  resolutions  better 
than  0. 1  micron.  Even  higher  precision  may  be  possible  through  the  use  of  fabrication 
facilities  offering  higher  resolution  photolithography. 

1.3. 3. 3  Optical  Position  Sensor 

A  lateral  effect  photodiode  (see  Fig.  1.23  and  Figs.  1.9-1.11)  was  used  in  the  two 
degree-of-freedom  experiments  described  earlier  to  monitor  the  position  of  a  quartz 
optical  Fiber.  In  this  use,  the  device  showed  excellent  stability,  reliability  and 
resolution,  as  configured,  in  the  micron  range.  The  micro  systems  being  considered 
now  require  much  smaller  detectors,  probably  fabricated  in  the  base  structure. 
Optically  sensitive  areas  such  as  photoconductive  areas  or  photodiodes  can  be 
fabricated  in  a  silicon  base  by  conventional  procedures. 

1.3. 3. 4  Capacitive  Position  Sensor 

The  first  studies  with  the  drive  strips  used  in  the  two  degree-of-freedom  experiments 
investigated  the  potentials  above  the  drive  strips  using  a  Fine  metal  wire  probe  (see  Fig. 
1.14).  The  probe  was  connected  to  the  input  of  a  very  high  input  impedance  amplifier 
with  which  the  voltage  on  the  probe  was  measured.  This  arrangement  is  an  example  of 
how  capacitive  interactions  have  been  used  to  characterize  a  position-dependent 
parameter.  The  system  acted  as  a  capacitive  voltage  divider  formed  by  the  probe 
capacitance  and  input  capacitance  of  the  amplifier  which  reduced  the  potential  measured 
by  the  amplifier. 

An  experimental  study  was  made  of  the  use  of  drive  strip  capacitance  for  position 
sensing.  A  charged  fiber  was  positioned  over  the  drive  strips  using  a  precision 
mechanical  stage.  A  low  leakage  op  amp  was  connected  with  a  capacitor  to  form  a 
charge  integrator.  The  current  flowing  to  and  from  the  drive  strip  was  integrated  so  that 
the  charge  on  the  strip  could  be  known  at  all  times. 

Two  measurement  configurations  were  used.  In  one,  the  drive  strip  was  held  at  ground 
potential  at  all  times.  In  the  second,  the  strip  could  be  raised  to  a  set  potential  and  used 
to  produce  forces  while  the  charge  was  monitored.  The  potential  of  the  other 
conductors  on  the  base  were  known,  so  the  charge  induced  by  them  in  the  strip  being 
monitored  could  be  added  or  subtracted  from  the  integrator  signal.  It  was  found  that 
the  charged  fiber  produced  induced  charge  levels  which  could  be  measured  with  about 
1%  accuracy  (see  Fig.  1.24).  The  position  of  the  fiber  could  then  be  calculated.  The 
speed  of  response  was  not  adequate  for  control  of  the  fiber,  however. 

The  results  of  this  experiment  provide  an  interesting  comparison  with  the  open-gated 
FET  in  measuring  induced  charge.  Both  methods  need  FETs,  but  one  locates  the  FET 
at  the  site  of  the  measurement,  while  the  other  locates  it  at  a  distance  (within  a  high 
impedance  amplifier)  and  connects  it  to  a  wire  probe.  The  integrator  suffered  from 
current  leakage  which  was  much  greater  than  that  of  the  FET  gate,  and  after  a  few 
seconds  was  die  major  source  of  error  in  the  measurement.  The  integrator  used  the 
whole  length  of  the  charged  area  to  make  its  measurement  while  the  FET,  when  close  to 
the  fiber,  monitored  only  about  \%  of  the  charged  length.  The  results  showed  the 
advantages  of  the  FET  measurement  and  supported  the  idea  that  local  measuring 
devices  are  superior. 
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1.3.3.S  Hall  Effect  Position  Sensor 

Hall  effect  devices  have  been  used  at  the  Center  for  Engineering  Design  to  measure 
joint  angles,  linear  position,  and  force.  In  each  application,  the  Hall  effect  device 
detected  variation  in  magnetic  field  intensity.  The  magnetic  field  was  modulated  by  a 
moving  element  connected  to  the  rotating  or  translating  body.  These  systems  have  been 
successful  in  the  Utah/MIT  Dextrous  Hand  (Jacobsen,  et  al.,  1984)  and  will  continue  to 
be  evaluated  for  use  in  micro  systems. 

1.3.4  Electret  Experiments 

1.3. 4.1  Geometry,  Materials,  Fabrication 

Sheets  of  Teflon®  (primarily  TFE  rather  than  FEP)  were  found  to  hold  fairly  high 
charge  densities  over  long  periods  of  time,  even  for  Teflon®  foils  as  thin  as  12.5 
microns.  The  Teflon®,  furthermore,  was  easy  to  work  with;  for  planar  geometries  no 
other  material  was  studied. 

For  the  two  degree-of-freedom  micro  servo  system  (described  above  in  Section  1 .2.6 1 
it  was  necessary  to  apply  thin  electret  coatings  to  140  micron  diameter  quartz  optical 
fibers.  Materials  studied  for  this  purpose  included  Teflon®  TFE,  polycarbonate, 
polyimide,  polystyrene,  and  acetophenone-doped  polystyrene.  None  of  these  materials 
had  clearly  and  repeatably  superior  electrical  properties  (maximum  charge  density, 
charge  stability,  etc.),  but  the  materials  differed  significantly  in  the  methods  that  were 
needed  to  apply  them  to  the  fibers.  Two  methods  for  producing  electret-coated  quartz 
fibers  proved  successful:  (1)  winding  a  20  micron  diameter  filament  of  Teflon®  TFE 
around  the  quartz  fiber;  and  (2)  dip  coating  the  fiber  with  a  3  to  10  micron  thick  coating 
of  polycarbonate  or  polystyrene  (see  Fig.  1.16).  Fibers  were  also  prepared  with  50()A- 
thick  vacuum  deposited  paladium-gold  alloy  surfaces  (Figure  1.16).  Their  behavior  in 
the  micro  servo  system  was  compared  to  the  fixed-charge  (electret)  fibers. 

1.3. 4. 2  Charging  Methods 

A  number  of  methods  have  been  reported  for  charging  electrets  (see  Bibliography, 
Section  I.B).  We  have  investigated  coronal-discharge  and  liquid-contact  methods,  but 
have  concentrated  on  electron  bombardment  with  a  scanning  electron  microscope 
(SEM)  which  allows  the  greatest  control  over  electron  kinetic  energy  (via  the  SEM 
accelerating  voltage)  and  the  rate  of  electron  implantation  (via  the  SEM  beam  current). 
In  the  SEM,  Teflon®r/i^erj  were  put  in  contact  with  a  grounded  metal  surface,  w  hich 
provided  charge  compensation  during  charging  (see  Fig.  1.25a).  Electron  beams  w  ith 
lOkV  to  30kV  provided  good  results  for  charge  density  and  stability  on  sheets.  The 
case  of  electret  fibers  was  much  more  delicate.  Successful  charging  required  the 
empirical  discovery  of  a  number  of  handling  techniques,  among  them:  (1)  while 
charging  in  the  SEM,  the  fiber  was  kept  in  a  V-groove  in  a  grounded  metal  plate  (which 
apparently  provided  temporary  charge  compensation,  see  Fig.  1.25b);  (2)  the  fiber 
while  charging  was  rotated  by  a  small  "rotisserie"  built  from  a  mechanical  watch  (see 
Fig.  1.26);  (3)  the  density  and  stability  of  charge  turned  out  to  be  enhanced  if  the  SEM 
charging  process  was  cycled,  with  the  SEM  vacuum  broken  between  charging  cycles; 
and  (4)  charge  stability  was  enhanced  with  appropriate  thermal  treatment  prior  to 
charging. 
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1.3. 4. 3  Charge  Densities  and  Stability 

•  Charge  densities  in  Teflon®  sheets,  averaged  over  lmm?  or  so,  were  measured  with  a 

Monroe  electrostatic  voltmeter.  Finer  spatial  resolution  of  charge  density  was  studied 
with  an  open-gated  FET  sensor  developed  at  CED  and  described  in  Section  1.3. 3. 2 
above.  Charge  densities  of  around  20  nC/cm2  were  found  (see  Figs.  1.27-1.29), 
equivalent  to  an  electric  field  of  order  105  volts/cm.  These  densities,  furthermore,  were 
stable  over  periods  of  at  least  months  (see  Fig.  1.30).  The  measured  charge  densities 

•  were  in  reasonable  agreement  with  published  charge  densities  (of  order  100  nC/crrP). 

The  FET  was  used  for  measuring  fine  scale  linear  charge  density  on  fibers.  Averaged 
linear  charge  densities  were  determined  by  measuring  the  elastic  deflection  of  the 
cantilevered  fiber  in  a  known  electric  field  produced  by  parallel  plate  electrodes.  Fiber 
linear  charge  density  as  high  as  0.4  nC/cm  was  measured,  equivalent  to  a  fiber  voltage 

9  of  a  few  hundred  volts.  Studies  of  the  stability  of  fiber  charge  are  not  yet  complete. 

There  is  evidence  that  the  charge  longevity  on  fibers  may  be  very  sensitive  to  ambient 
humidity,  and/or  to  mechanical  flexing.  Initial  results  suggest  that  significantly  greater 
charge  density  and  stability  can  be  achieved  with  optimized  thermal  treatment  recipes. 
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1.4  Resources 

1.4.1  Important  Areas  of  Development 


Since  the  inception  of  this  project,  it  was  understood  that  progress  in  the  development  of 
micro  systems  would  involve  the  capability  to  design,  fabricate,  and  evaluate  micro 
systems.  Specifically,  the  following  capabilities  have  been  identified  as  important: 


1 .  Developing  the  analytical  and  intuitive  capabilities  required  to  permit  the 
design  of  systems  which  have  at  least  some  chance  of  working. 
Activities  include  classical  analysis,  simulations,  and  the  graphic 
representation  of  fields  and  forces. 


2.  Developing  the  capability  to  manufacture,  handle,  and  evaluate  electrets, 
and  other  electrically  active  materials  (such  as  ferroelectrics).  These 
include  SEMs  for  electron  implantations,  and  deposition  of  electron- 
retentive  materials  such  as  Teflon®  or  poly  carbon  ate. 


3.  Developing  and/or  accessing  fabrication  systems  necessary  for  the 
production  of  micro  systems.  These  include  IC  techniques,  such  as 
mask  making,  etching,  deposition,  ion  implantation,  etc. 


4.  Developing  techniques  for  the  experimental  evaluation  of  micro 
systems.  This  includes  high  performance  electronics,  micro  probes, 
controlled  environments,  etc. 


1.4.2  Project  Personnel  •  Areas  of  Expertise 


The  following  personnel,  grouped  as  faculty,  staff,  students  and  consultants,  where 
involved  in  the  research  reported  herein. 


Faculty 


NAME 

RANK 

DEPT. 

AFFILIATION 

AREA  OF 

EXPERTISE 

S.C.  Jacobsen* 

Professor 

ME,  BE,  CS 

Electro-mechanical 
systems  design 

J.E.  Wood** 

Associate 

Professor 

BE,  ME 

Biomechanics,  physio¬ 
logy,  muscle  mechanics 

R.E.  Price 

Professor 

Physics 

Physics,  electromagnetics, 
black  holes 

S.G.  Meek 

Assistant 

Professor 

ME 

Electro-mechanical 
systems  design 

*  PRINCIPAL  INVESTIGATOR 
**  CO-PRINCIPAL  INVESTIGATOR 


I 
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Staff 


NAME 

DEGREE 

DEPT. 

AFFILIATION 

AREA  OF 
EXPERTISE 

C.R.  Bowier 

CED 

Electronics 

Technician 

N.W.  Clayton 

M.S.  Degree 
Physics 

CED 

Physics, 

electronic  materials 

Jack  Corey 

Micro  circuits 
fabrication 

B.K.  Hanover 

M.E.  Degree 

Mech.  Engr. 

CED 

Electro- mechanical 
design 

T.E.  Hansen 

M.S.  Degree 

Elect.  Engr. 

CED 

Micro  circuits 
fabrication 

R.T.  Johnson 

M.S.  Degree 

Elect.  Engr. 

CED 

Electronics  systems 
&  analog  design 

D.F.  Knutti 

M.S.  Degree 

Mech.  Engr. 

CED 

Electro-mechanial 
systems  design 

Max  Levy 

B.S.  Degree 
Biology 

HEDCO 

Laboratory 

Micro  circuits 
fabrication 

E.M.  Simon 

M.S.  Degree 
Bioengineering 

CED 

Polymer 

fabrication 

K.  Wilhelmsen 

CED 

Electronics 

technician 
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Students 


NAME 

DEGREE 

(DATE) 

DEPARTMENT 

AREAS  of 

INTEREST 

E.J.  Cetron 

Ph.D. 

(1988) 

Bioengineering 

Electro- magnetic 
field  displays 

B.  Eastham 

Computer 

Science 

CAD 

K.W.  Grace 

M.S. 

(1987) 

Mechanical 

Engineering 

Electro- mechanical 
systems  control 

P.S.  Khanwilkar 

M.S. 

(1987) 

Bioengineering 

AC  stabilization 
modelling 

R.P.  Phillips 

Ph.D. 

(1988) 

Bioengineering 

Solid-state  sensors 

T.S.  Shelby 

0988) 

Materials 

Science  Engineering 

FET  drift 

D.  Trembath 

M.S. 

(1986) 

Materials 

Science  Engineering 

Electret  fabrication 

Consultants 


NAME 

title; 

AFFILIATION 

AREA  OF 

EXPERTISE 

J.D.  Andrade 

Professor 

BE,  MSE 

Materials,  interfaces 

J.  Hansen 

Scientist 

AD 

Microfabrication 

F.  Harris 

Professor 

Physics 

E/M  field  computations 

R.J.  Huber 

Professor 

EE 

microdevice  physics 
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1.4.3  Equipment  and  Facilities 

Having  identified  the  capabilities  required  to  design  and  produce  micro  system  devices, 
used  equipment  and  personnel  to  meet  the  need. 

By  way  of  equipment,  we  acquired: 

Computation: 

The  Microfield  Project  utilized  the  following  digital  computation 
resources  for:  (1)  Data  acquisition,  (2)  Management  of  experiments,  (3) 
Implementation  of  control  functions,  (4)  Simulation,  and  (5)  Display  of 
information. 

1 .  Three  MicroVAX  Os  running  Micro  VMS  4.4 

2.  One  VAX  11/750  running  4.3  BSD  Unix 

3.  Two  MINC  1 1  computers  with  analog  interfaces;  using  1 1/73  and 
1 1/23  CPU's  running  RSX-1 1M+ and  RT-1 1 

4.  One  PDP  1 1/44  with  analog  interfaces,  running  RSX-1 1M+ 

5.  ComputerVision  CAD  system,  CIT  467  color  display  terminal,  HP 
7220C  8  pen  plotter 

6.  HP  9000/300  "Bob-Cat"  processor/display  running  HP-UX 

7 .  Wasatch  high  resolution  (1024  x  780)  color  display 

8 .  VS1 1  high  resolution  color  display 

Mask  Design: 

a.  CADDS-2  software  for  Computer  Vision  design  station. 

b.  Our  own  software  to  edit/modify  pattern  generator  files. 

1C  Testing: 

a.  Access  to  the  Semiconductor  Device  Analyzer  in  the  Integrated 
Circuits  Instrumentation  Laboratory  in  the  Computer  Science 
Department . 

b .  Access  to  JEOL  SEM  facilities  at  Hercules  for  inspection  of  wafers. 

1C  Fabrication: 

a.  Diffusion  furnaces  (pad  oxidation) 

i)  Nitride  stress  reduction  for  local  oxidation  process 

b.  Diffusion  furnaces  (nitride  deposition) 

i)  Nitride  deposition  to  prevent  oxidation  of  active  areas  during 
field  oxidation. 

c.  Diffusion  furnaces  (field  oxidation) 

i)  Used  for  isolation  of  active  circuits. 

d .  Diffusion  furnaces  (gate  oxidation) 

i)  To  form  the  gate  of  the  MOS  transistors. 
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e.  Diffusion  furnaces  (phosphorus  doping) 

i)  To  form  the  N+  source/drain  the  MOS  transistor  to  dope  the 
polysilicon  layer  N+,  and  to  dope  the  contacts  N+  so  that  the 
contacts  are  omhic. 

f .  Diffusion  furnaces  (Anneal  cycle) 

i)  To  Anneal  the  Aluminum  to  make  the  circuit  active. 

g.  Diffusion  furnaces  (polysilicon) 

i)  Used  to  deposit  polysilicon  for  the  formation  of  the  gates  of 
the  transistors. 

h.  Ion  implant 

i)  Boron:  to  drive  the  field  areas  more  p-type  to  prevent 
formation  of  field  transistors. 

ii)  Phosphorus:  Implanted  in  gate  to  form  depletion  mode 
transistors. 

i.  CVD  deposition 

i)  Used  to  isolate  layer  of  the  integrated  circuit. 

j.  Evaporation 

i)  Use  to  deposit  aluminum  onto  the  wafers. 

k.  Aligners 

i)  Use  to  mask  the  layer  of  the  design  layout  which  will 
eventually  for  the  integrated  circuits. 

l.  Pattern  generator 

i)  Used  to  form  the  lOx  reticle  of  the  circuit  design. 

m.  Step  and  repeat  camera 

i)  Used  to  reduce  the  lOx  reticle  and  to  form  the  mask  set  that 
will  be  used  to  fabricate  the  integrated  circuit. 


Electret  Fabrication: 

a.  Access  to  the  JEOL  SEM  facilities  at  Hercules  for  implantation  and 
inspection  of  electrets. 

b.  Access  to  the  Hitachi  SEM  facilities  at  the  Department  of  Geology. 

c.  Access  to  wet-lab  facilities  for  the  chemical  fabrication  of  polymers  to 
be  made  into  charged  electrets. 

d.  Availability  of  plasma  deposition  (Teflon)  facilities  at  the  University 
of  Washington,  and  Crest  Coatings,  Inc.  in  Los  Angles. 
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Experimentation: 

a.  Optical  stages  and  LEPD 

b.  Faraday  cage 

c.  High-voltage  amplifiers 

d.  Storage  oscilloscope  (Tektronix) 

e.  3  degree-of-freedom  stage  for  FET  tests 

f .  Monroe  electrostatic  voltmeter 

g.  Standard  electronic  apparatus  (power  supplies,  DVMs,  function 
generators,  etc.). 

h.  Field  validation  apparatus 
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Fig.  1.1.  An  apparent  potential  surface.  For  a  two  degree-of-freedom  electrostatic 
system  with  velocity-independent  active  feedback  the  force  is  nonconservative  and  cannot 
be  represented  by  a  potential  surface.  It  can,  however,  be  represented  by  an  apparent 
potential  surface  with  a  slope  that  gives  the  correct  direction,  but  not  magnitude,  of  the 
force  field.  [  1986] 


SURFACE  CHARGES 


Fig.  1.2a.  A  cross-section  of  an  electret  showing  different  types  of  charge.  1 19X5 


i 
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Fig.  1.2b.  An  electret  generated  by  implantation  of  electrons  in  a  25-micron  thick 
Teflon®  sheet  by  a  computer-controlled  SEM  beam  tracing  the  pattern  of  the  letters.  Xerox 
toner  (particle  size  5-20  microns),  which  is  attracted  to  charges,  was  used  to  indicate  the 
location  of  the  charged  regions.  Smallest  letters  are  approximately  250  microns  in  height. 
Smallest  line  width  shown  is  less  than  35  microns.  Actual  width  of  line  charge  is  probably 
somewhat  less  than  this.  [1985] 


Driver  Voltages 


Fig.  1.3.  An  illustration  of  a  classical  feedback  approach  to  control  fields  through 
sensing,  control  logic,  and  active  feedback.  [1985] 
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Fig.  1.4.  Schematic  of  electret  element  above  gate  of  FET.  Induced  Field  on  gate 
modulates  current  output  of  FET.  [1985] 
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Fig.  1.6.  A  plot  of  drain  current  I<j  of  a  field  effect  transistor  as  a  function  of  the  height 
of  a  charged  electret  fiber  above  the  sensitive  gate  region.  [1984] 


Fig.  1.7.  A  specialized  field  effect  transistor  developed  by  the  CED  group  and 
manufactured  in  the  HEDCO  Laboratory  for  measuring  electric  fields.  [1985] 
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Fig.  1.8.  A  schematic  representation  of  an  experiment  set-up  to  validate  fields  (potentials 
and  forces)  and  to  control  a  straight  filament-triad  system.  Potentials  Ej  can  be  applied  to 
base  driver  strips,  to  move  filament  or  to  vary  fields.  The  system  will  also  monitor  induced 
potentials  on  sensor  strips  Sj.  The  vector  quantities  are  also  listed.  [1985] 
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Fig.  1.9.  A  realistic  sketch  of  a  "laboratory  size"  apparatus  that  is  suitable  to  study  small 
scale  fields  that  allow  rapid  fabrication  and  manual  adjustment.  [1985] 
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Fig.  1.10.  A  photograph  of  the  actual  apparatus  used  to  study  small  scale  fields  on  an 
electret  fiber  with  two  degrees-of-freedom.  Note  the  conductive  strips  below  the  fiber  and 
the  LEPD  light  sensor  for  measuring  the  position  of  the  end  of  the  fiber.  [1986] 


c'< 


Fig.  1.11.  A  photograph  of  an  apparatus  with  one  degree-of-freedom  actuated  electret 
fiber  used  to  determine  the  charge  density  of  the  electret  fiber  between  the  two  parallel 
plates.  The  LEPD  senses  the  end  point  position  of  the  electret  fiber.  [1986] 


Fig.  1.13.  An  illustration  of  the  base  strips  with  volt  ge  drivers  (D)/senso  amplifiers 
(S),  with  an  electret  filament  above.  [  1985] 
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Fig.  1.15.  A  mathematical  model  of  field  intensities  generated  by  the  base  strips.  [1985] 
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Fig.  1.16a.  An  illustration  of  the  three  different  fibers  used  in  the  various  experiments. 
Figure  (a)  shows  a  140pm  quartz  fiber  with  a  vapor  deposited  aluminum  coating  to  provide 
a  conductive  sheath;  (b)  shows  a  quartz  fiber  coated  with  a  polycarbonate  sheath  that  when 
bombarded  with  electrons  from  an  SEM  forms  the  electret;  and  (c)  shows  a  Teflon® 
wrapped  quartz  fiber.  [1986] 


Fig.  1.16b.  SEM  photograph  of  an  actual  Teflon®  wrapped  quartz  fiber.  [1985] 
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Fig.  I.j7a.  Photograph  of  a  five-pointed  "star"  taken  from  an  oscilloscope  generated  by 
a  two  degree-of -freedom  electret  fiber  under  feedback  control.  [1987] 


Fig.  1.17b.  Photograph  of  the  letters  "DARPA"  off  an  oscilloscope  screen  generated  by 
a  two  degree-of-freedom  electret  fiber  under  feedback  control.  [1987] 
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Fig.  1.18.  Schematic  of  forces  acting  on  cantilevered  filament  set  above  conductive 
strips.  [1986] 


Fig.  1.19.  Simulated  fiber  trajectories  with  active  control  feedback.  Part  (a)  shows  the 
result  of  underdamped  control,  and  Part  (b)  shows  overdamped  control.  [1986] 


VFRT  DIST 


Fig.  1.20.  Voltage  configuration  for  AC  stabilization.  [1985] 
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Fig.  1.21.  Simulation  of  AC  stability:  Case  (a)  shows  trajectory  for  inherently  unstable 
initial  fiber  position,  with  inadequate  AC  stabilization:  case  (b)  shows  AC-stabilized 
trajectory.  [1986] 
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Fig.  1.22.  The  stiffening  of  a  fiber  by  an  AC  driving  force.  Curve  (a)  shows  the  natural 
oscillation  of  the  fiber  with  no  AC  force  applied.  Curve  (b)  represents  the  same 
experimental  conditions  except  that  an  AC  force  has  been  added;  the  AC  force  slightly 
increases  the  natural  oscillation  frequency.  [1987] 


(a) 


(b) 


Fig.  1.23.  Schematic  of  LEPD.  Part  (a)  shows  relationships  of  voltages  and  positions; 
Part  (b)  shows  LEPD  cross-section  with  p-i-n  layers.  [1985] 


MicroFields  Final  Report 


37 


Induced  Charge  (pico  Ceil) 


A'ctan  (W/2#h)  (degrees) 


Fig.  1.24.  Induced  charge  in  a  conducting  sensing  strip,  measured  by  an  integrator  as  a 
function  of  h,  the  height  of  a  charge  filament,  and  W,  the  strip  width.  [1985] 


(a)  0») 


Fig.  1.25.  Electron  beam  charging  configuration  for  (a)  electret  sheet,  and  (b)  electret 
coated  fiber.  [1984] 


Fig.  1.26.  Mechanical  rotisserie  for  rotating  fiber  during  charging  inside  the  SEM. 


[1985] 


Fig.  1.27.  Relative  linear  charge  density  for  lines  of  charge  implantation  in  50  micron 
nonmetallized  Teflon®  TFE  foil,  shown  as  a  function  of  beam  sweep  rate,  with  a  30kV 
electron  beam  current  of  0.2nA.  [1986] 
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Fig.  1.28.  Linear  charge  density,  X,  vs.displacement,  x,  across  Teflon®  sheet  having 
lines  of  charge,  each  line  implanted  at  various  beam-sweep  rates.  Data  taken  with  FET. 
[1986] 


Fig.  1.29.  Linear  charge  density,  X,  vs.  displacement,  x,  along  length  of  charged 
Teflon®- wrapped  fiber.  Data  taken  with  FET.  [1986] 


Siftn*  (nC/cm) 


MicroFields  Final  Report 


40 


Fig.  1.30.  One  dimensional  scan  with  an  electrostatic  voltmeter  across  a  Teflon®  foil 
charged  in  two  regions:  (a)  1  day  after  charging;  (b)  45  days  after  charging.  [1986] 
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2 .  Basic  research  Studies 

The  following  Sections  provide  expanded  detail  on  some  of  the  topics  presented  in  Section 
1.  These  basic  research  studies  where  meant  to  study  and  resolve  some  of  the  technical 
issues  necessary  to  be  overcome  in  order  that  Micro-Field  systems  could  be  fabricated. 
Some  of  these  "studies"  where  also  involved  in  the  design  and  fabrication  of  system 
subelements,  such  as  the  generation  of  a  number  of  different  FET  devices  for  different 
sensor  studies  and/or  system  requirements.  The  "Designer  System",  discussed  in  Section 
2.3,  is  a  CAD-based  system  which  allows  a  "Micro-Field"  system  designer  to  analyze  and 
visualize  the  complex  3D  vector  and  potential  fields  generated  by  any  proposed  micro¬ 
geometry  system.  It  is  important  to  have  faithful  analytical  representations  of  small-scale 
systems  since  they  are  typically  too  small  to  "tinker"  with,  and  too  costly  to  make 
numerous  trial-and-error  prototypes. 

2.1  Electrets  as  Field  Emitters 

Electrets,  which  are  intended  to  be  a  primary  source  of  non-uniform  fields  to  be 

detected  by  FETs  or  to  be  "force-source"  elements,  where  investigated  extensively. 

Some  of  our  results  are  discussed  below. 

2.1.1  Fabrication 
2 . 1 . 1 . 1  Sbeets 

Sheet  electrets  were  formed  by  implanting  thin  films  of  Teflon  with  30-kY 
electrons  from  the  beam  of  a  scanning  electron  microscope  (SEM).  Teflon  was 
selected  because  of  its  high  bulk  resistivity,  as  well  as  the  fact  that  its  structure 
provides  a  high  density  of  stable  trapping  sites,  where  implanted  electrons  can 
remain  with  a  half-life  of  decades.  Sheet  electrets  fabricated  du^ng  the 
Microfield  project  fall  mosdy  into  one  of  the  following  categories: 

1)  Uncompensated  unpatterned  sheets.  These  were  usually  formed  by 
mounting  a  non-metallized  film  of  Teflon,  usually  on  the  order  of  one 
mil  (25  microns)  in  thickness,  on  a  polished  metal  surface  which  was 
then  grounded  in  the  SEM  chamber.  The  electron  beam  was  then 
rastered  over  a  rectangular  patch  of  the  surface,  covering  anywhere  from 
10  to  100  mm2  (see  Fig.  1.25a).  The  rastering  was  continued  until 
saturation  was  achieved;  it  was  found  that,  regardless  of  how  high  the 
dose  was,  surface  charge  densities  tended  to  have  an  upper  limit, 
estimated  at  100  nC/cm2.  After  implantation,  the  Teflon  could  be 
removed  from  the  metal. 

2)  Compensated  unpatterned  sheets  (see  Fig.  2 .la).  Almost  all  of  the 
electrets  describe  in  the  literature  are  of  this  type.  For  the  Microfield 
project,  we  formed  such  electrets  by  rastering  a  30-kV  electron  beam,  as 
described  below,  except  that  the  rear  surface  of  the  Teflon  was 
aluminized  and  grounded,  thus  forming  an  integral  and  permanent 
source  of  compensation  charge. 

3)  Patterned  sheets  (see  Fig.  2.1b).  Patterned  electrets  were  formed  by 
stretching  a  sheet  of  Teflon  (usually  aluminized  on  the  back)  over  a 
special  holder  which  maintains  surface  tension  on  the  material  (similar  to 
a  drumhead).  Control  of  the  SEM  beam  scanning  is  then  taken  over  by 


an  external  drive  circuit.  The  beam  is  swept  over  the  surface  of  the 
Teflon  in  spot  mode,  as  opposed  to  the  usual  raster  mode.  The  electron 
beam  can  be  made  to  "write"  a  specified  pattern  of  charge  on  the  electret, 
in  a  manner  analogous  to  E-beam  lithography.  The  linear  charge  density 
implanted  in  the  Teflon  can  be  controlled  by  varying  the  beam  current, 
the  scan  speed  of  the  beam,  or  both.  Most  of  our  fabrication  effort 
consisted  of  making  parallel  lines  of  charge  using  a  simple  analog  beam 
driver.  This  was  sufficient  for  studying  the  effects  of  various 
fabrication  techniques  on  such  electret  characteristics  as  charge 
distribution,  charge  stability,  and  the  effect  of  non-uniform  fields  on 
FET  field  sensors.  A  more  sophisticated  computer-driven  beam-driving 
package  was  developed  and  used  on  occasion  to  write  more  complicated 
patterns,  such  as  alphanumerics. 

4)  Hybrid  Materials.  Some  of  the  design  ideas  conceived  during  the 
Microfield  project  are  best  implemented  by  directly  integrating  a  field 
emitter  (such  as  an  electret)  with  the  field-sensing  device.  For  example, 
a  strain  sensor  might  be  fabricated  by  attaching  an  electret  directly  to  an 
FET  chip  by  "gluing"  the  two  together  with  a  thin  layer  of  elastomer 
(see  Fig.  2.1c).  We  successfully  fabricated  such  an  electret  by  thus 

attaching  a  sheet  of  uncharged,  unaluminumized  Teflon®  to  a  silicon 
chip,  then  implanting  patterns  of  charge.  The  silicon  acts  essentially  as 
a  compensating  conductive  layer,  hence  these  are  compensated  electrets, 
with  strong  fields  occurring  in  the  elastomer  between  the  sheets  where 
they  can  be  readily  detected  by  the  FET. 

It  should  be  noted  that  the  fabrication  of  patterned  electrets  requires  that  careful 
attention  be  paid  to  certain  issues  which  are  not  important  in  the  manufacture  of 
unpattemed  electrets,  such  as  focusing  of  the  beam  on  the  surface  of  the  Teflon 
(which  must  be  done  without  the  benefit  of  SEM  imaging,  as  this  would 
implant  unwanted  electrons  over  a  large  area)  and  alignment  of  the  pattern 
relative  to  specific  reference  points. 

2. 1.1. 2  Filaments 

For  certain  Microfields  applications,  such  as  SCOFSS,  we  required  charged 
fibers  or  filaments.  In  the  case  of  SCOFSS,  we  needed  an  optical  fiber  which 
had  net  electrical  charge  along  a  certain  portion  of  its  length  (see  Fig.  2. Id). 
Considerable  work  was  done  to  find  a  method  of  coating  part  of  a  quartz  fiber 
with  an  electretable  material.  The  major  types  are  listed  below: 

1)  Teflon-wound  fiber.  A  quartz  optical  fiber  with  a  diameter  of  140 
microns  was  used  (see  Fig.  1.16c).  A  thin  (12-25  micron)  filament  of 
pure  TFE  Teflon  was  wound  tightly  around  the  optical  fiber  over  a 
length  of  about  1  cm.  The  filament  was  then  heat-treated  and  exposed  to 
a  rastered  30-kV  electron  beam  in  an  SEM.  Charge  densities  on  the 
order  of  0.4  nC/cm  were  achieved.  However,  this  method  was  very 
labor-intensive  and  time-consuming. 

2)  Teflon-coated  fiber.  A  quartz  fiber  was  dipped  in  a  Teflon  emulsion 
consisting  of  sub-micron  spheres  of  pure  Teflon  suspended  in  liquid 
(see  Fig.  1.1b).  The  liquid  was  allowed  to  evaporate,  then  the  Teflon 
coating  was  heat-treated  to  produce  sintering  (in  the  case  of  TFE  Teflon) 
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or  fusion  (in  the  case  of  FEP  Teflon).  Such  coatings  are  tricky  to  work 
with,  however,  and  we  were  unable  to  achieve  the  same  degree  of 
charge  retention  as  with  the  Teflon  fiber. 

3)  Polycarbonate  and  polyimide-coated  fibers.  Methods  were 
developed  to  obtain  quickly  and  with  a  minimum  of  labor  a  thin,  high- 
purity  coating  of  polycarbonate  or  polyimide  on  a  quartz  fiber.  Both  of 
these  materials  are  good  for  making  electrets,  although  the  charge 
trapping  mechanism  becomes  rather  unstable  in  the  presence  of 
moisture;  hence  great  care  must  be  taken  to  fabricate  and  store  such 
electrets  in  a  dry  environment.  Charge  densities  as  as  high  as  0.4 
nC/cm  were  obtained. 

In  the  process  of  developing  fabrication  technology,  we  discovered  several 
interesting  phenomena,  two  of  which  which  turned  out  to  be  of  great  practical  use 
in  the  manufacturing  of  uncompensated  fiber  electrets.  The  first  of  these  is  the  fact 
that  a  much  higher  charge  density  can  be  obtained  by  providing  a  means  of 
temporary  compensation  charge  during  the  charging  process.  Early  attempts  at 
charging  fibers  produced  charge  densities  on  the  order  of  0.01  nC/cm,  which  is  an 
order  of  magnitude  less  than  what  was  later  achieved.  These  early  electrets  were 
charged  simply  by  suspending  the  polymer-coated  fiber  in  space  and  hitting  it  with 
the  rastered  electron  beam.  Theoretical  considerations  led  us  to  suspect,  however, 
that  the  electrostatic  field  resulting  from  the  buildup  of  a  net  negative  charge  on  the 
fiber  might  be  capable  of  deflecting  subsequently  incident  electrons.  Hence,  efforts 
were  made  to  design  a  fiber  holder  which  would  greatly  reduce  the  strength  of  this 
field.  The  most  successful  design  consists  of  a  grounded  metal  plate  with  a  90- 
degree  V-shaped  groove  milled  in  its  surface  to  a  depth  approximately  equal  to  the 
diameter  of  the  fiber.  The  fiber  is  then  laid  in  this  groove  and  exposed  to  the 
electron  beam  (see  Fig.  1.25b).  The  resulting  charge  density  was  an  order  of 
magnitude  greater  for  fibers  of  all  materials. 

Another  important  discovery  was  that  the  charge  retention  of  polycarbonate-coated 
fibers  could  be  greatly  enhanced  by  charging  the  fiber  once  in  the  vacuum 
environment  of  the  SEM,  then  exposing  it  to  air  at  atmospheric  pressure,  followed 
by  another  charging  in  vacuum.  This  can  boost  the  charge  density  by  a  factor  of  3 
or  more,  and  also  results  in  much  better  longevity  of  the  charge.  Polycarbonate 
fibers  thus  treated  and  subsequently  stored  and  used  in  a  normal  environment  (i.e. 
normal  room  air)  retained  their  net  negative  charge  with  a  half-life  of  several  weeks. 
Such  fibers  were  adequate  for  use  in  SCOFSS  and  other  experimental  apparatus. 
Further  testing  will  be  required  to  determine  the  longevity  of  such  a  fiber  when 
exposed  to  an  oscillating  electric  field  in  an  ideal  environment,  such  as  vacuum. 

2.1.2  Compensated  vs.  Uncompensated 

It  is  noted  that  virtually  all  of  the  thin-film  electrets  described  in  the  literature  were 
coated  with  metal  on  one  or  both  faces.  In  the  most  common  configuration,  the 
backside  of  the  polymer  film,  i.e.  the  surface  opposite  that  on  which  electrons  were 
to  be  deposited,  is  coated  with  aluminum,  which  is  grounded  during  the  charging  of 
the  electret.  As  charge  is  implanted  in  the  front  surface,  charges  of  the  opposite 
polarity  flow  into  the  metal  from  ground  and  form  a  compensating  layer  (see  Fig. 
2.2a).  The  result  is  that  the  electrostatic  field  outside  the  electret  is  weakened, 
analogous  to  a  dipole.  At  the  same  time,  the  internal  field  between  the  two  layers  of 
charge  is  strengthened.  Such  an  electret  is  referred  to  as  "compensated"  because  of 
the  compensating  layer  of  charge  which  forms  in  the  metal.  A  fully  compensated 
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electret  has  zero  net  charge,  because  the  compensating  charge  is  of  the  same 
magnitude  as  the  charge  deposited  on  the  polymer. 

"Uncompensated"  electrets,  i.e.  ones  with  a  net  charge  and  hence  a  strong 
monopolar  external  field,  are  given  very  little  treatment  in  the  literature  (recently,  the 
German  audio  company  Senhizer  has  introduced  a  line  of  electret  headphones  that 
use  uncompensated  electret  sheets  in  a  hermetically  sealed  cell).  Unlike  their 
compensated  cousins,  these  electrets  have  no  grounded  metal  backing  to  provide 
compensation  charge  (see  Fig.  2.2b).  We  have  had  experience  with  the  fabrication 
of  both  types  as  part  of  the  Microfield  project. 

Compensation  affords  certain  advantages  during  the  fabrication  and  storage  of  an 
electret.  Our  electrets  were  fabricated  by  implanting  charge  in  the  polymer  surface 
with  an  electron  beam  in  a  scanning  electron  microscope  (SEM).  When  the  back 
side  of  the  polymer  is  metallized  and  grounded,  the  field  above  the  front  surface  of 
the  film  is  kept  relatively  weak  by  the  canceling  effect  of  the  compensation  charge. 
This  allows  the  electrons  in  the  SEM  beam  to  traverse  this  space  without  being 
deflected  or  backscattered  by  strong  E-fields.  Hence,  more  charge  can  be  implanted 
in  the  polymer.  If  no  source  of  compensation  charge  is  provided,  the  net  charge 
buildup  occurring  in  the  film  produces  a  strong  external  field  which  can  deflect 
incoming  electrons,  thus  reducing  the  maximum  attainable  charge  density  in  the 
electret. 

During  storage  and  actual  use,  an  uncompensated  electret  attracts  charged  particles 
from  the  environment,  which  can  attach  themselves  to  the  surface  of  the  polymer 
and  thus  cancel  the  implanted  charge.  A  compensated  electret,  on  the  other  hand,  is 
much  less  prone  to  such  cancellation,  due  to  the  fact  that  the  "ground-based” 
cancellation  substantially  weakens  any  particle-attracting  fields. 

A  compensated  electret  standing  alone  in  free  space  will  produce  a  weak  external 
field,  which  may  cause  such  an  electret  to  appear  poorly  suited  to  Microfield 
applications.  However,  the  situation  changes  when  a  conducting  surface  is  brought 
close  to  the  charged  surface  of  the  electret.  If  the  conductor  is  electrically  connected 
to  the  metallized  rear  surface  of  the  electret,  a  strong  field  is  generated  between  the 
surface  of  the  electret  and  the  surface  of  the  external  conducting  plane,  as  simple 
analysis  using  image  charges  readily  shows.  Thus,  for  example,  we  can  get  a 
strong  reading  from  one  of  our  FET  field  detecting  chips  simply  by  grounding  the 
substrate  of  the  chip  to  the  back  of  a  compensated  sheet  electret,  and  then  bringing 
the  chip  close  to  the  charged  polymer  surface.  The  actual  sensing  element  occupies 
only  a  tiny  fraction  of  the  surface  area  of  the  chip;  thus,  it  views  the  surrounding 
surface  as  an  infinite  grounded  conducting  plane  and  can  sense  the  strong  field 
occurring  between  the  chip  and  the  electret  However,  the  chip  essentially  shields 
the  outside  world  from  this  field.  Hence,  the  external  E-field  is  still  comparatively 
weak,  while  the  e-field  between  the  interacting  Microfield  elements  -  in  this  case  the 
FET  field  sensor  chip  and  the  electret  -  is  strong.  This  is  a  lever-law  relationship, 
as  discussed  in  Jackson  (1975). 

Uncompensated  electrets  can  best  be  fabricated  by  holding  the  polymer  temporarily 
against  a  conductive  backing  during  the  charging  process,  thus  imitating  a 
compensated  electret.  However,  this  conductor  may  be  removed  after  charging  to 
leave  an  uncompensated  electret  with  a  net  charge.  In  Microfieids  applications, 
such  an  electret  is  likely  to  be  used  in  close  proximity  to  conductors,  which  provide 
a  source  of  compensation  charge  and  thus  shield  the  outside  world  from  the  strong 
E-field  which  is  present  close  to  the  electret.  Such  an  electret  is  thus  actually  a 
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compensated  electret,  after  a  fashion.  As  long  as  the  space  in  which  the  E-field  is 
strong  is  kept  clean  and  dry  (preferably  vacuum),  the  possibility  exists  that  such  an 
electret  will  be  as  stable  as  a  true  compensated  electret.  However,  in  our  research 
we  found  that  uncompensated  electrets  were  more  difficult  to  fabricate  in  a 
reproducible  fashion  than  were  truly  compensated  ones.  Nonetheless,  we  did  in 
fact  fabricate  a  large  number  of  fiber  electrets  (see  Section  3  and  Fig.  1.16b)  which 
were  of  the  true  uncompensated  variety  and  which  were  not  used  in  a  well- 
compensated  environment.  These  were  employed  with  great  success  in  the  Small 
Cantilevered  Optical  Fiber  Servo  System  (SCOFSS)  described  in  Section  3  below. 
Such  electrets  exhibited  charge  half-lives  of  several  weeks,  when  kept  in  room  air. 

2.1.3  Patterns 

Virtually  all  of  the  electrets  describe  in  the  literature  are  of  the  "unpattemed"  variety, 
meaning  that  the  surface  of  interest  is  charged  in  a  gross  (uniform)  fashion.  For 
example,  sheet  electrets  charged  with  an  electron  beam  in  an  SEM  are  typically 
charged  by  rastering  the  beam  over  a  large  area  of  polymer  film.  Such  electrets  are 
useful  for  applications  in  w-hich  a  uniform  electrostatic  field  is  required,  for  example 
in  electret  microphones  or  headphones.  However,  the  Microfield  concept 
encompasses  a  variety  of  designs  which  require  the  generation  of  non-uniform 
electric  fields  on  a  small  spatial  scale.  One  means  of  doing  so  is  to  make  an  electret 
with  a  non-uniform  charge  distribution,  such  as  a  patterns  of  lines  or  grids  (see  Fig 
1.2b).  Consequently,  a  considerable  effort  was  directed  toward  learning  how  to 
fabricate  and  handle  patterned  sheet  electrets;  at  the  same  time,  a  better 
understanding  of  the  characteristics  of  such  electrets  was  achieved. 

Most  of  the  patterned  electret  generation  was  done  by  "drawing",  with  the  SEM 
beam,  series  of  parallel  lines  of  charge  several  mm  in  length  and  separated  by 
distances  ranging  from  100  microns  to  5  mm  (see  Fig.  2.3).  The  material  most 
commonly  used  was  TFE  Teflon  film  aluminized  on  one  side,  with  a  thickness  of 
25  microns;  hence,  the  electrets  were  mostly  of  the  compensated  variety.  Linear 
charge  dosages  ranged  from  0. 1  nC/cm  to  c  nC/cm. 

The  thus  generated  patterns  were  studied  by  placing  the  electret  in  a  specially 
designed  stage  which  could  be  used  to  bring  the  charged  surface  within  10  micron’s 
of  a  floating-gate  FET  field  sensor  developed  by  CED  for  this  verv  purpose. 
Rectangular  field  sensors  as  small  as  2  microns  by  40  microns  are  avaiLK1e  on  the 
sensor  chip.  The  sensor  was  then  passed  over  the  parallel-line  charge  pattern,  and 
its  response  was  recorded  as  a  function  of  position.  The  result  was  a  series  of 
sharp  spikes  corresponding  to  the  location  of  the  line  charges.  When  tne  sensor 
was  located  over  one  of  the  line  charges,  displacements  as  small  as  0.1  micron 
resulted  in  an  observable  change  in  FET  source-dram  current  over  a  bandwidth  of 
20  MHz,  due  to  the  change  in  electric  field.  Reduction  of  bandwidth  by  means  of 
low-pass  filtering  could  allow  for  improved  resolution. 

One  of  the  interesting  advantages  of  electrets  as  field  generators  was  emphasized 
during  our  studies  of  patterned  compensated  electrets.  We  constructed  a  device 
which  was  in  every  way  identical  to  a  patterned  electret,  except  that  instead  of  lines 
of  implanted  charge,  we  used  metal  wires  25  microns  in  diameter  stretched  in 
parallel  over  the  surface  of  the  Tefion  film  and  driven  by  a  high-voltage  DC  power 
supply.  We  found  that  field  strengths  similar  to  those  achieved  with  an  electret 
could  be  produced,  but  required  the  use  of  at  least  600  volts  DC.  Suci.  high 
voltages  would  be  impractical  in  most  of  the  micro  devices  which  we  envision. 
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2.1.4  Non-Uniformities 

While  we  were  studying  compensated  Teflon  electrets  implanted  with  parallel  lines 
of  charge,  we  found  a  that  the  E-field  sensed  by  the  FET  field  sensor  was  10%- 
20%  stronger  at  the  lines  on  either  end  of  the  series  than  at  the  lines  in  between. 
The  resulting  plot  of  field  strength  as  a  function  of  position,  produced  as  the  FET 
was  scanned  across  the  lines,  has  the  appearance  of  an  inverted  set  of  fangs  (see 
Fig.  2.1  below).  Consequendy,  we  dubbed  this  phenomenon  the  "fang  effect”. 

The  fang  effect  turned  out  to  be  highly  reproducible,  given  that  the  various 
fabrication  parameters  (beam  current,  beam  voltage,  beam  scan  rate,  etc.)  were 
reproduced,  and  given  that  the  FET  was  scanned  over  the  same  cross-section  of  the 
electret  pattem.each  time.  A  considerable  experimental  effort  was  aimed  at 
understanding  and  controlling  the  effect.  A  brief  summary  of  results  is  presented 
here: 


1 )  The  existence  of  the  effect  does  not,  to  any  detectable  degree,  depend  on 
the  order  in  which  the  lines  of  charge  are  implanted. 

2)  The  effect  is  not  due  to  the  geometry  of  the  experimental  arrangement, 
i.e.  a  pattern  in  which  all  lines  of  charge  were  of  identical  linear  charge 
density  would  not  produce  a  fang  effect  of  the  magnitude  observed.  Hence, 
the  effect  reflects  the  charge  distribution  in  the  electret  as  well  as  the  E-field 
close  to  the  surface. 

3)  The  effect  begins  to  die  out  significantly  only  when  the  lines  are 
implanted  several  millimeters  apart.  Conversely,  the  closer  the  line  spacing, 
the  more  pronounced  the  effect. 

4)  The  effect  appears  to  be  pioduced  during  or  very  shortly  after  the  actual 
implantation  process.  Thereafter,  it  is  very  stable  with  time. 

5)  The  effect  can  be  reduced  or  even  made  to  disappear  by  changing  the 
electron  dose  from  the  SEM  beam. 

The  last  observation  is  of  primary  importance,  as  we  shall  likely  find  desirable  the 
manufacture  of  patterns  of  line  charge  having  a  common  uniform  linear  charge 
density.  We  observed  that  the  effect  appeared  most  prominently  in  electrets  which 
were  at  the  saturation  limit,  i.e.  ones  in  which  the  maximum  attainable  linear  charge 
density  had  been  achieved.  Increasing  the  dose  beyond  this  point  did  not  result  in 
increased  linear  charge  density,  but  invariably  produced  a  strong  fang  effect.  As 
the  dose  was  lowered  below  the  saturation  value,  the  effect  lessened,  and 
completely  disappeared  below  0.2  nC/cm  to  0.4  nC/cm.  Interestingly,  the  dose  at 
which  the  effect  disappeared  was  found  to  depend  somewhat  on  whether  that  dose 
was  achieved  via  adjustment  of  beam  current  or  via  adjustment  of  beam  scan  rate. 
The  linear  charge  density,  it  should  be  noted,  is  not  a  linear  function  of  dose,  as 
many  of  the  incident  electrons  are  not  stably  trapped  in  the  Teflon  and  therefore 
escape  quickly.  The  linear  charge  density  resulting  from  the  dose  at  which  fangs  no 
longer  appeared  was  about  half  to  two-thirds  the  maximum  attainable  linear  charge 
density  (see  Fig.  2.4). 

Ultimately,  more  work  must  be  done  in  order  to  fully  understand  the  fang  effect. 
We  have  encountered  no  reference  to  such  an  effect  in  the  literature,  but  this  is  not 
surprising;  few  people  have  studied  patterned  electrets,  and  to  our  knowledge 
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nobody  has  characterized  a  sheet  eleccret  in  the  way  we  have,  i.e.  by  measuring  the 
E-field  very  close  to  the  surface  with  a  high-resolution  FET  probe  whose  critical 
dimension  is  on  the  order  of  a  few  microns.  In  the  meantime,  the  fang  effect  can  be 
effectively  controlled  by  controlling  charge  dose. 

2.2  FETs  as  Sensor  Elements 

The  majority  of  devices  envisioned  in  the  Microfield  project  require  that  we  either 
procure  or  develop  very  small  devices  which  can  sense  the  strength  of  electric  or 
magnetic  fields.  We  arrived  early  on  at  the  conclusion  that  such  devices  ought  to  be 
implementable  via  integrated  circuit  technology.  One  reason  for  this  is  the  requirement 
that  the  devices  be  cheaply  fabricated  in  large  numbers;  in  many  device  designs,  we 
envisioned  having  arrays  of  many  field  sensors  spaced  microns  apart.  A  further 
motivation  was  our  desire  to  eventually  incorporate  signal  processing  circuitry  on  the 
same  chip  as  the  sensor  elements.  Such  designs  would  be  impossible  to  implement 
with  discreet  components. 

While  Hall-effect  devices  may  be  suitable  for  magnetic  measurements,  no  such  off-the- 
shelf  micro-device  was  found  to  fill  the  need  for  electrostatic  sensing.  We  therefore 
developed  a  design  for  a  floating-gate  FET;  i.e.  one  in  which  the  gate  is  electrically 
isolated.  Such  a  configuration  will  cause  the  source-drain  current  in  the  FET  to 
respond  almost  linearly  to  the  strength  of  the  E-field  at  the  gate  surface,  assuming  that 
the  device  is  fabricated  with  the  correct  threshold  voltage  (see  Fig.  1.4).  Temperature 
compensation  can  then  be  achieved  by  incorporating  a  second  FET  on  the  chip  (referred 
to  as  a  "reference"  device),  with  its  gate  connected  to  a  constant  voltage.  The  difference 
between  the  source-drain  currents  of  the  two  FETs  is,  to  first  order,  insensitive  to 
temperature  changes  in  the  chip,  while  remaining  sensitive  to  external  E-fields.  Finally, 
we  concluded  that  FETs  with  a  high  aspect  ratio  were  desirable,  by  virtue  of  their  being 
more  sensitive  to  E-field  strength. 

Unfortunately,  off-the-shelf  FET  designs  do  not  offer  the  above  features. 
Furthermore,  the  well-calibrated  processes  most  commonly  used  for  custom  integrated- 
circuit  fabrication  do  not  provide  the  threshold  voltage  (approximately  -1.5  volts)  which 
we  required.  We  therefore  relied  initially  on  an  experimental  FET  device  called  a 
"ChemFET",  which  happened  to  meet  most  of  our  requirements  (see  section  2.2.1 
below).  Later,  we  designed  a  series  of  field-sensing  FET  chips  which  were  custom 
fabricated  by  Analog  Devices,  Inc.,  using  a  trial-and-error  method  to  obtain  the  unusual 
threshold  voltage  required.  Descriptions  of  these  devices  are  given  below. 

2.2.1  ChemFETs 

Early  experiments  in  the  use  of  floating-gate  FETs  as  a  means  of  sensing  the 
strength  of  an  electrostatic  field  were  performed  with  a  device  known  as  a 
ChemFET.  The  ChemFET  was  designed  at  the  University  of  Utah  Department  of 
Bioengineering  for  use  as  a  chemical  sensor,  and  was  fabricated  in-house  at  the 
HEDCO  microelectronics  laboratory  at  the  University.  It  is  an  N-channel  depletion¬ 
mode  device.  A  few  of  the  lot  were  found  to  possess  the  correct  threshold  voltage 
for  our  application  (-1.5  volts).  These  were  given  to  us  without  the  final  layer  of 
aluminum  applied;  we  made  them  work  as  field  sensors  by  depositing  the  top  layer 
of  aluminum  over  the  FET  channel  regions  in  such  a  way  that  the  metal  was  not 
electrically  connected  to  any  voltage  sources.  In  this  way,  a  floating-gate  device 
was  produced;  slight  changes  in  die  gate  potential  about  zero  (with  respect  to 
ground)  could  be  induced  by  an  external  electric  field,  and  these  changes  would 
produce  a  reasonably  linear  change  in  the  source-drain  current. 
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A  single  ChemFET  chip  has  two  parallel  floating-gate  FETs  at  one  end  (see  Fig. 
1.7).  The  active  sensing  area  of  each  is  approximately  20  microns  by  400  microns. 
Two  similar  FETs  lie  at  the  other  end  of  the  die,  but  these  have  gates  which  are 
conductively  connected  to  bonding  pads,  and  are  therefore  useful  as  reference 
devices  for  temperature  compensation.  Temperature  compensation  was  found  to  be 
so  good  that  the  compensated  signal  did  not  change  when  the  entire  device  was 
heated  to  approximately  120  degrees  Fahrenheit. 

Most  of  the  early  studies  of  electrets  were  done  using  the  ChemFET.  We  were  able 
to  obtain  a  response  to  a  line  charge  suspended  in  space  over  the  floating  gate  (see 
Fig.  1.5)  which  corresponded  very  well  with  that  predicted  by  a  simple  model  of 
such  a  system  (see  Fig.  1.6).  Noise  levels  in  the  device  were  on  the  order  of  0.1 
microamps  over  a  bandwidth  of  20  MHz.  By  passing  the  device  over  the  surface  of 
an  electret  with  lines  of  implanted  charge,  we  were  able  to  sense  the  change  in  the 
E-field  with  displacement  with  such  sensitivity  that  movements  as  small  as  0.1 
microns  were  detectable  over  the  same  high  bandwidth.  Short-term  repeatability 
was  also  found  to  be  excellent.  However,  we  noticed  that  the  source-drain  current 
in  the  device  had  a  tendency  toward  long-term  drift  of  the  total  current  of  as  much 
as  30%,  characterized  by  an  average  period  of  several  days.  Furthermore,  any 
object  coming  in  contact  with  the  gate  metal  (which  was  not  covered  with  an 
insulator)  immediately  produced  a  sudden  change  in  FET  current,  due  to  a  change 
in  the  gate  potential.  Also,  we  were  able  to  find  only  a  small  number  of  devices 
suited  to  our  peculiar  needs  in  the  lot  given  us  by  the  ChemFET  group.  In  addition, 
the  floating  gates  of  the  field-sensing  FETs  were  too  large  to  allow  us  to  see  details 
of  E-field  structure  with  as  much  spatial  resolution  as  we  would  have  liked. 
Finally,  only  the  "off-the-shelf'  device  geometry  was  available,  whereas  we  had 
designed  a  number  of  Microfield  devices  which  incorporated  geometrically  novel 
arrays  of  FETs.  Thus,  while  the  ChemFET  was  of  immense  utility  in  confirming 
theoretical  predictions  and  testing  certain  design  ideas  for  Microfield  devices,  it  was 
of  limited  use  in  constructing  actual  prototypes  of  microsensors;  nor  could  it  be 
used  to  conduct  the  entire  array  of  experiments  in  which  we  were  interested.  With 
these  issues  in  mind,  we  set  out  to  design  our  own  family  of  floating-gate  FET 
chips  to  be  custom-fabricated  by  an  outside  facility.  These  are  described  in  the 
following  sections. 

2.2.2  SG-Chips 

One  of  the  goals  of  the  Microfield  project  was  to  lay  the  groundwork  for  the  design 
and  construction  of  a  prototype  strain  sensor  which  uses  the  relative  displacement 
of  a  small  field-emitting  structure  and  a  field-sensing  structure  in  order  to  infer 
strain.  One  method  of  achieving  this  objective  is  to  use  a  field  emitter  (either  an 
electret  or  a  driven  conductor)  consisting  of  a  single  narrow  strip  located  a  few 
microns  away  from  one  or  more  rectangular  floating-gate  FET  E-field  sensors. 

To  this  end,  we  designed  an  FET  chip  which  we  call  the  "SG"  chip  (for  "strain 
gauge").  The  design  consists  of  three  pairs  of  floating-gate  FETs  (see  Fig.  2.8). 
The  gates  within  each  pair  are  identical.  All  FETs  on  the  chip  have  the  same 
channel  length  (the  short  dimension  of  the  FET),  but  different  pairs  have  different 
widths  (the  long  dimension)  and  different  spacings.  The  FETs  in  a  given  pair  are 
parallel  and  are  separated  by  no  more  than  a  few  tens  of  microns.  All  gates  have 
high  aspect  ratios.  This  feature  is  of  importance  in  that  the  sensitivity  of  a  floating- 
gate  FET  to  E-field  strength  is  proportional  to  the  aspect  ratio  of  the  device. 
Ideally,  one  wants  a  gate  which  is  small  in  the  direction  of  current  flow  (i.e.  the 
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"length"  of  the  channel)  and  large  in  the  orthogonal  dimension  (i.e.  the  "width"  of 
the  channel). 

The  SG  chips  were  designed  to  be  used  in  a  "triad"  configuration  in  which  a  single 
field-emitting  strip  is  held  above  the  gates  and  parallel  to  them.  If  one  then 
subtracts  the  currents  of  the  two  FETs,  one  obtains  a  single  current  which  is  related 
uniquely  to  the  horizontal  position  of  the  field  emitter  between  the  two  gates.  The 
greatest  sensitivity  to  lateral  displacement  of  the  field  emitter  occurs  when  the 
emitter  lies  roughly  halfway  between  the  two  gates  at  a  height  equal  to  half  the 
separation  of  the  FETs.  Maximum  resolution  of  displacement  is  of  great 
importance  with  the  SG  chip,  because  a  strain  sensor,  in  order  to  compete  with 
current  commercially  available  strain  gauges,  must  be  capable  of  detecting  a 
microstrain;  this  is  equivalent  to  a  displacement  of  5  nm  for  a  gauge  length  of  5 
mm. 

The  SG  chip  contains  three  reference  FETs,  one  corresponding  to  each  size  of 
floating-gate  field-sensing  FET. 

2.2.3  K-Chips 

A  major  thrust  of  the  Microfield  project  is  the  study  of  electrets  as  field  emitters. 
Early  studies  were  done  with  a  ChemFET  of  the  type  described  in  Sec.  2.2.1 
above.  Unfortunately,  the  ChemFET  had  the  major  drawback  that  its  gates  were 
rather  large  (20  microns  by  400  microns)  and  were  thus  not  capable  of  resolving 
smaller  details  of  field  structure.  Furthermore,  only  one  ChemFET  chip  with  the 
correct  threshold  voltage  was  available  to  us.  Therefore,  we  designed  a  chip  to  use 
specifically  for  the  purpose  of  studying  electrets  and  their  interactions  with  FETs. 
The  result  was  the  so-called  K-chip  (arbitrarily  named),  which  has  four  different 
sizes  of  floating-gate  FETs  on  board,  along  with  reference  FETs  of  the  same  sizes 
(see  Fig.  2.7). 

The  smallest  FET  on  the  K  chip  is  approximately  2  microns  by  40  microns.  The  2- 
micron  length  is  the  lower  limn  of  the  fabrication  technology  used  by  the 
manufacturer.  The  40-micron  wi,  ih  is  intended  to  provide  a  high  aspect  ratio,  and 
thus  a  high  sensitivity  to  external  E-fields.  The  three  other  sizes  are  approximately 
5  by  100  microns,  10  by  200  microns,  and  20  by  400  microns  respectively.  A 
second  2-by-40-micron  device  is  located  at  the  other  end  of  the  array.  By  using  the 
two  small  FETs  simultaneously,  one  can  achieve  optimum  parallel  alignment  of  a 
linear  field  emitter  -  such  as  a  line  or  fiber  implanted  with  charge  -  with  the  FET 
array  on  the  K-chip. 

The  K-chip  has  its  floating-gate  FETs  located  close  to  the  end  of  the  die,  with  the 
bonding  pads  on  the  opposite  end.  This  is  to  allow  the  FETs  to  be  brought 
arbitrarily  close  to  the  electret  under  study.  The  wires  bonded  to  the  pads  protrude 
at  least  25  microns  above  the  surface,  even  when  held  flat  against  the  chip.  This 
would  have  prevented  the  FETs  from  being  brought  any  closer  to  an  electret  than  25 
microns,  had  they  been  located  close  to  the  pads.  By  tilting  the  entire  chip  at  a 
small  angle,  however,  the  FETs  at  the  far  end  can  be  brought  into  direct  contact 
with  the  surface  of  an  electret,  while  the  wires  bonded  at  the  other  end  are  given 
clearance. 

Shown  in  Fig.  2.6  is  a  plot  of  FET  differential  drain  current  as  a  function  of 
position.  This  data  was  obtained  from  one  of  the  2-by-40-micron  FETs  of  a  K-chip 
by  placing  it  about  20  microns  away  from  the  surface  of  a  sheet  of  Teflon  implanted 
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with  parallel  lines  of  charge.  The  FET  was  then  moved  across  the  pattern  in  a 
direction  perpendicular  to  the  lines  of  charge,  with  the  long  dimension  of  the  gate 
parallel  to  them.  The  distance  between  lines  is  1000  microns.  The  plot  is  actually 
composed  of  multiple  scans,  indicating  excellent  repeatability. 

2.2.4  Octa-FETs 

Another  device  envisioned  in  the  Microfield  project  is  a  very  small  multi-axis  sensor 
which  can  be  used  to  detect  small  displacements  in  multiple  degrees  of  freedom. 
These  displacements  may  then  be  used  to  infer  such  quantities  as  strain  or  force, 
depending  on  how  the  device  is  coupled  to  the  outside  world.  In  order  to  provide  a 
testbed  for  such  a  device,  we  designed  a  chip  which  has  eight  annular-gate  FETs 
arranged  in  a  square  pattern  as  shown  in  Fig.  2.9  below.  This  design  was 
arbitrarily  named  the  "Q-FET". 


An  E-field  emitter,  such  as  an  electret  or  a  conductive  pattern,  may  be  fabricated  in 
the  form  of  a  cross,  such  that  one  branch  of  the  cross  lies  between  the  gates  of  each 
pair  of  FETs.  The  center  of  the  cross  is  aligned  with  the  center  of  the  FET  pattern 
(see  Fig.  2.10).  The  field-emitting  device  may  be  rigidly  attached  to  or  integrated 
with  a  larger  armature  which  is  coupled  to  the  medium  of  measurement.  A 
displacement  of  the  field  emitter  in  any  of  three  translational  degrees  of  freedom  or 
three  rotational  degrees  of  freedom  will  result  in  a  change  in  the  output  of  at  least 
some  of  the  FETs.  For  a  given  pair  of  FETs  with  a  linear  field  emitter  located 
between  them  and  a  few  microns  above  the  surface  of  the  chip,  the  difference  of  the 
two  FET  currents  is  a  function  of  the  horizontal  displacement  of  the  field-emitting 
line.  At  the  same  time,  the  sum  of  the  two  FET  outputs  is  a  function  of  the  height 
of  the  emitter  above  the  surface  of  the  chip.  Thus,  each  of  the  four  pairs  provides 
two  pieces  of  information,  for  a  total  of  eight  numbers.  From  these  eight  numbers, 
the  position  of  the  armature  in  all  six  degrees  of  freedom  can  be  determined  by 
algebraic  means  for  small  displacements. 

A  unique  feature  of  the  Q-chip  design  is  the  incorporation  of  annular  gates  instead 
of  the  conventional  rectangular  geometry.  Such  a  design  provides  an  output  which 
is  uniformly  sensitive  to  the  horizontal  displacement  of  a  nearby  linear  field  emitter 
for  all  angles.  In  other  words,  a  linear  field  emitter  located  a  fixed  distance  above 
the  chip  may  be  moved  toward  or  even  over  a  single  FET,  and  the  response  of  the 
device  will  depend  only  on  the  horizontal  distance  between  the  center  of  the  FET 
gate  and  the  line;  it  will  not  depend  on  the  horizontal  angle  of  the  emitter.  This 
considerably  simplifies  the  inversion  mathematics  required  to  extract  small- 
displacement  6-degree-of-freedom  information  from  the  outputs  of  the  eight  FETs, 
as  the  algorithm  for  doing  so  relies  on  knowing  the  displacements  of  the  field 
emitting  fines  relative  to  the  various  FETs  regardless  of  the  horizontal  angle. 

Two  Types  of  Q-chips  were  produced.  The  two  differ  only  in  the  size  scale  and 
spacing  of  the  annular  gates.  The  Q1  chip  has  gates  with  a  diameter  of  50  microns, 
lying  along  the  sides  of  a  square  of  length  800  microns.  For  the  Q2  chips,  these 
dimensions  are  80  microns  and  800  microns  respectively.  Center-to-center  spacing 
of  adjacent  gates  is  75  microns  for  the  Q1  chip  and  125  microns  for  the  Q2  chip. 
Each  type  of  Q  chip  has  a  single  reference  FET  located  near  the  bonding  pads. 
Only  one  reference  device  is  needed,  as  all  FETs  on  the  chip  possess  identical 
geometry. 
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2.2.5.  Mathematical  Studies  and  Simulations 
2. 2. 5.1.  Theory  of  Operation 

The  FETs  used  in  our  initial  work  had  thermal  oxide  layers  on  the  order  of  500 
A  between  the  gate  and  the  substrate.  The  electric  fields  (typically  104  volts/cm) 
in  which  the  FETs  are  to  be  used  can  produce  gate-substrate  voltage  differences 
of  order  10*2  volts.  In  studying  the  electric  field  configuration  produced  by 
field  generators  (electrets,  driven  conductors,  etc.)  and  measured  by  FETs  this 
10" 2  volt  difference  is  negligible.  In  its  effect  on  the  field  the  FET  therefore  can 
be  idealized  as  being  located  on  a  conducting  plane  (the  FET  chip).  The  FET 
design  and  external  electronics  serve  to  monitor  with  great  precision  the 
integrated  normal  E  field  (or  equivalently  the  induced  charge)  on  one  portion  of 
the  FET,  the  gate. 

Mathematical  investigations  of  FET  behavior  in  field  sensing  involve  the 
solution  of  the  field  problem  with  the  field  emitter  in  the  neighborhood  of  a 
conducting  plane,  and  the  integration  of  the  normal  E  field  on  a  region  of  that 
plane.  The  approach  used  has  been  to  find  an  analytic  solution  (a  closed  form, 
a  series,  an  integral,  etc.)  of  the  field  problem,  and  to  use  a  numerical  evaluation 
of  that  solution  where  appropriate.  A  fully  numerical  approach  to  solving  field 
problems  has  not  been  necessary  due  to  the  relatively  simple  field  geometries 
dealt  with.  Furthermore  a  fully  numerical  approach  would  require  highly 
specialized  software,  which  has  not  yet  been  acquired,  or  the  very  time- 
consuming  development  of  our  own  finite  element  or  variational  code  for 
solving  a  field  problem. 

As  examples  of  the  application  of  analysis  to  the  use  of  FETs  as  field  sensors 
we  report  below  three  areas  where  mathematical  modelling  was  done  to  support 
experimentation  and  design:  (i)  The  prediction  of  the  sensitivity  of  gates  in  a 
strain-gauge  configuration  with  conducting  emitters;  (ii)  An  analysis  of 
dielectric  effects  on  a  charged  line  embedded  in  an  electret,  to  determine  whether 
such  effects  could  explain  "fangs"  (see  §  2.1.4);  (iii)  A  study  of  the 
performance  of  annular  gates  as  sensors  of  electreted  lines. 


2. 2. 5. 2.  Sensing  of  Long  Conductive  Emitters 

The  most  basic  configuration  for  FETs  as  one-dimensional  position  sensors  is 
shown  in  Fig.  2.1 1.  The  grounded  conducting  plane  at  y  =  0  represents  the 
FET  chip.  Parallel  to  it  at  a  height  y  =  h  (typically  15  microns)  is  a  conducting 
plane  which  is  grounded  except  for  two  conducting  strips  which  are  maintained 
at  potentials  Vi,  V2,  V3...  These  conducting  strips,  and  the  FET  gates  are 
considered  to  be  very  long  (typically  larger  than  200-300  microns)  in  the  z- 
direction  so  that  the  field  problem  can  be  approximated  as  two  dimensional.  A 
relative  displacement  of  the  two  planes  results  in  a  change  of  the  signal  (i.e.,  the 
integrated  normal  E  field)  on  the  gate,  and  provides  the  basis  for  sensing  that 
displacement. 

This  two-dimensional  field  problem  is  straightforward  to  solve  with  the 
Dirichlet  Green  function  for  the  geometry,  i.e.,  the  solution  to  the  problem 
shown  in  Fig.  2.12:  The  field  of  a  line  charge  of  charge-per-unit  length  1  at 

position  x’,  y’.  The  solution  for  the  potential  $  is  found  e.g.,  with  a  Schwarz- 
Christoffel  transformation,  to  be 
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X  ,  . -i 

sin  (*£)  sin  (k  J) 

tann 

2n£o 

.cosh  (*  ip)  ~  cos  (*  h)  cos  (n  h). 

(In  addition  to  its  importance  in  connection  with  FET  sensing,  this  solution  has 
direct  experimental  applicability:  Electreted  fibers  were  placed  in  the  field 
between  parallel  conducting  plates,  with  a  large  potential  difference.  The 

deflection  of  the  fibers  gave  a  measure  of  their  linear  charge  density  X). 

The  above  solution  (with  X  =  1)  can  now  be  used  in  a  standard  Green  function 
method  to  find  the  solution  to  the  problem  pictured  in  Fig.  2.13: 


<D(x,y)  =  Jo(x',y’)  ^7  3(x,x';  y,y')  dA\ 


4tc  s' 


9n’ 


Here  the  integration  is  ever  the  bounding  surfaces  at  y’  =  0,  h,  and  9G/3n’ 
represents  the  outward  normal  derivative.  The  number,  location,  and  potentials 

of  the  conducting  strips  is  contained  in  the  specification  of  <J>(x’,y’)  in  the 

integral.  In  the  case  of  a  single  driven  conducting  strip,  at  O  =  Vi,  extending 
from  x  =  xi,  to  x  =  X2  the  potential  is  found  to  be: 


<t>(x,y)  = 


Xi 

jt 


|  tan’1  [tanh  n  (^“)tan  71  ^J-tan*1  [tanh 


tan  rc 


From  this  result  it  is  straightforward  to  find  the  charge  density 

•  oOO  --TT  [tanh  (TE1)-  ,anh  *  (TIT1)] 

induced  in  the  y  =  0  FET  plane.  The  FET  signal  is  simply  the  integral  of  o(x) 
over  the  extent  in  x  of  the  gate,  multiplied  by  the  width  of  the  gate  (i.e.,  the  z- 

•  extent  of  the  gate).  That  is,  the  FET  signal  is  proportional  to  the  total  induced 

charge. 

2. 2.5.3.  Dielectric  Effects  on  Electreted  Lines 

Unexpected  effects  were  found  when  FETs  were  used  to  measure  the  linear 

•  charge  density  X  of  each  of  a  set  of  parallel  lines  of  charge  embedded  with  an 

SEM  in  Teflon®  with  metallized  backing  (see  §  2.1.4.).  The  outermost  lines 

were  measured  to  have  higher  X  than  inner  lines.  This  motivated  a  closer  look 
at  the  field  of  electreted  lines  and,  especially,  whether  dielectric  effects  could 
play  a  significant  role. 


The  typical  problem  to  be  studied  is  pictured  in  Fig.  2.13,  idealized  as  a  two- 
dimensional  configuration  in  the  x,y  plane.  The  grounded  plane  at  y  =  c 
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represents  the  FET  chip,  and  that  at  y  =  -b  represents  the  metallized  backing  of 
the  Teflon®.  The  Teflon®  slab,  with  relative  dielectric  constant  e  *  2,  extends 

from  y  =  a  to  y  =  -b  with  the  electreted  line,  of  linear  charge  X,  located  at  the 
coordinate  origin. 

To  solve  this  problem  a  technique  was  developed  that  applies  to  any  problem 
containing  a  line  charge  and  parallel  dielectric  or  conducting  boundaries.  The 
problem  represents  each  horizontal  slab  (i.e.,  each  region  bounded  by  y  = 
constant  surfaces),  in  terms  of  a  cosine  transform,  as 


<b(x,z)  =  Je'tycosCkxWk^  +  Je+kycos(kx)\jr(k)dk 
o  o 

In  the  slab  containing  the  line  charge  (at  the  origin)  the  potential, 


^  X  fe'k'ytoskx 

<D  =  -  - r - dk 

2rce  J  k 

o 

of  a  line  charge  in  an  unbounded  medium  of  dielectric  constant  E  is  added.  For 

each  slab  region  there  are  two  unknown  functions  <J>(h)  and  \|/(k).  These  are 
determined  by  the  boundary  conditions  at  each  interface  (e.g.,  V  =  0  at 

grounded  surface,  V  continuous  at  dielectric  interface,  edV/dy  continues  at 
dielectric  interface).  The  problem  is  thereby  reduced  to  an  algebraic  problem  for 

the  <j>  and  y  functions.  As  the  number  of  interfaces  grow,  the  algebra  becomes 
more  intricate,  but  the  method  remains  applicable.  As  an  illustration  we  can 
consider  the  problem  shown  in  Fig.  2.13  with  the  conductive  layers  removed, 
representing,  e.g.  (a)  line  charge  embedded  in  an  isolated  Teflon®  sheet 
without  metal  backing.  The  solution  above  the  sheet  (that  is,  for  y  >  a)  is  found 
to  be 


<D  = 


2K£n 


(1-P)  fc°sfae.ty(1  +  Mtb) 
H  J  k  i_R2e-2k(a+b) 


dk 


where 


The  solution  can  in  principle  be  used  directly  in  this  cosine  transform  but  for 
most  purposes  it  is  more  useful  to  interpret  the  solution  more  physically.  The 

solution  is  expanded  in  the  parameter  (5  (which  is  guaranteed  to  be  less  than 
unity).  Each  term  in  this  expansion  is  interpreted  as  an  image  line  charge.  For 
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example,  there  is  an  image  line  charge  of  charge  per  unit  length  X  (l-(3)  at  y  = 

0,  one  of  charge  per  unit  length  X  (l-(3)p  at  y  =  -2b,  etc.  All  effects  can  be 
viewed  as  due  to  die  superposition  of  the  actual  line  charge  and  the  infinite  set 
of  image  lines.  The  cosine  transform  becomes,  in  effect,  a  method  of 
organizing  an  image  line  charge  solution,  and  the  image  solution  is  very  useful; 
for  most  results  only  a  few  of  the  images  need  to  be  included.  When  this 
method  is  applied  to  the  problem  pictured  in  Fig.  2.13,  and  typical  numbers  are 
inserted  for  the  dielectric  constant  of  Teflon®,  the  depth  a,  of  the  electret  line, 
the  slab  thickness  (a  +  b),  and  the  distance  (c  -  a)  to  the  FET  sensor,  we  find 
that  the  effect  of  the  dielectric  on  the  FET-measured  field  is  not  crucial.  In 
particular,  it  cannot  be  the  cause  of  the  observed  "fangs.” 


2. 2. 5. 4.  Characteristics  of  Annular  Sensing  Gates 

Calculations  above  refer  to  rectangular  gates  with  high  aspect  ratios.  Also 
constructed  (see  §2.2.4.)  were  annular  gates  which  would  have  very  different 
patterns  of  sensitivity  to  field  emitters  (see  Fig.  1.7).  To  assess  the  nature  of 
these  sensitivity  patterns,  the  annular  gate  signal  (the  induced  charge  on  the 
gate)  was  computed  in  the  case  that  the  field  emitter  is  a  long  line  electret, 
parallel  to  the  grounded  FET  plane.  The  geometry,  and  parameters,  are 
illustrated  in  Fig.  2.14. 

The  charge  density  a  at  any  point  in  the  FET  plane  is  simply  found:  The  d>  field 
is  that  of  the  line  electret  and  its  image  mirrored  in  the  FET  plane;  the  electric 
field,  and  hence  a,  follows  from  taking  the  z  =  derivative  of  d>  at  z  =  0.  The 
resulting  integration  to  get  the  gate  signal  Qind  can  be  put  in  the  form 

Qind  =~  [F(Ro)-F(Ri)] 

K 

where 

2jt 

F(R)  = 

« 

o 

The  evaluation  of  the  F  integral  is  enormously  simplified  by  the  fact  that  it,  and 
hence,  Qind.  have  a  closed-form  solution.  This  follows  from  transforming  the 
trigonometric  integral  to  an  integral  on  the  unit  circle  in  the  complex  plane. 
With  this  transformation  the  denominator  of  the  integrand  becomes  a  fourth 
order  polynomial  which  can  be  explicitly  factored.  The  value  of  the  integral  is 
given  by  the  residues  at  the  two  poles  of  the  integrand  which  lie  inside  the  unit 
circle. 

The  results  for  Qind.  as  a  function  of  geometric  parameters,  shows  sensitivity 
patterns  very  much  unlike  those  for  a  rectangular  gate.  The  surfaces  of  equal 
Qind  turn  out,  for  a  wide  range  of  parameters,  to  be  almost  spherical. 


sin2<txJ<)> 
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This  calculation  can  be  extended  to  deal  more  realistically  with  electrets.  The 
idealized  line  electret  can  be  replaced  with  a  linear  deposition  of  charge  in  a 
metal-backed  sheet  of  Teflon®.  We  have  seen  in  §2.2.5.3.  that  the  field  from 
this  emitter,  parallel  to  a  FET  plane,  can  be  adequately  represented  as  the  line 
charge  augmented  with  several  images.  The  method  of  the  present  section  can 
be  applied  separately  to  each  of  those  images,  with  the  net  result  following  by 
superposition. 

2.3  Simulation  and  Design 

The  'CED  Electrostatic  Design  System’  is  an  integrated  package  which  runs  on  the 
Digital  Equipment  Corporation  VAXstation  II/GPX  system.  It  provides  a  menu  and 
mouse  driven  system  which  can  be  used  by  the  designer  to  test  new  device 
configurations.  With  simple  commands,  the  designer  can  create  an  armature/substrate 
configuration  and  proceed  to  use  several  analysis  techniques  for  examining  issues  such 
as  stability  and  controllability. 

2.3.1  Conceptual  Framework 

Being  able  to  generate  graphs,  plots,  and  images  is  not  enough.  The  generation  of 
the  geometric  model  must  be  simple  and  straight  forward.  Once  the  model  is 
created,  the  designer  needs  to  be  able  to  pick  and  choose  from  a  menu  which  types 
of  analyses  are  desired.  The  system  should  basically  step  the  designer  through  in  a 
natural  and  logical  way.  The  system  should  cater  to  the  inherent  creativity  of  the 
designer,  yet  all  of  the  low  level  code  to  calculate  forces  and  torques  should  be  as 
mathematically  rigorous  as  possible. 

A  second  major  factor  in  developing  such  a  design  tool  is  to  make  it  flexible  and 
extensible.  The  programs  and  modules  must  be  easily  replaceable  in  the  event  the 
focus  of  research  changes.  As  an  example,  if  work  were  to  progress  away  from 
electrostatics  into  magnetostatics,  only  three  small  subroutines  would  need  to  be 
changed.  All  the  rest  of  the  display  and  dispatching  routines  would  be  unchanged. 
The  replaced  modules  comprise  only  about  1000  lines  of  code  out  of  the  over 
15,000  which  make  up  the  Design  System.  In  this  way,  major  shifts  of  direction 
(magnetics  rather  then  electrets),  geometry  (helixes  instead  of  systems  with  linear 
elements),  or  display,  can  easily  be  accommodated. 

2.3.2  The  Organization  of  the  Design  System 

Fig.  2.15  shows  the  organizational  tree  of  the  Design  System.  The  lowest  level 
routines  are  shown  at  the  bottom  of  the  page,  and  the  higher  level  user-interface 
shown  at  the  top.  Major  changes  in  the  nature  of  the  elemental  interactions,  such  as 
a  shift  to  magnetostatics,  would  involve  changing  the  lowest  level  routines,  while 
user-interface  enhancements  would  require  the  modification  of  higher  level 
routines. 

The  hierarchy  readily  separates  into  five  distinct  areas: 

1)  Low  level  routines  -  for  field  calculations  and  determination  of  forces 
and  torques  between  armature  fibers,  conductive  drive  plates,  sensors, 
and  the  substrate. 

2)  Model  input  -  generation  of  the  initial  armature  and  substrate 
configuration. 

3)  Mapping  display  manager  -  which  displays  graphs  and  plots  of  fields 
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and  sensor  sensitivities  over  space  (two  dimensional  or  six 
dimensional). 

4)  Dynamic  display  manager  -  which  displays  images  of  the  structures  is 
motion  over  time. 

5)  Static  display  manager  -  which  interacts  with  the  designer  to  generate 
several  different  'snapshots'  of  the  system  state. 

The  basic  flow  of  control  within  the  system  is  top-down.  The  designer  chooses  to 
either  input  a  new  model,  or  to  analyze  a  previous  model.  Tf  the  designer  wishes  to 
input  a  new  model  the  model  input  section  is  run,  utilizing  several  of  the  low-level 
code  shown  at  the  bottom  two  levels  of  the  chart 

Should  the  designer  wish  to  perform  an  analysis,  one  of  the  three  analysis  mode 
managers  would  execute.  These  managers  (one  each  for  snapshots  of  the  model 
system,  time  course  behavior  of  the  model  system,  or  region  analysis  of  the  model 
system)  coordinate  through  the  analysis  dispatchers  to  call  the  various  low-level 
routines.  This  hierarchical  approach  allows  each  analysis  manager  to  use  the  same 
analysis  dispatchers,  thus  reducing  the  amount  of  computer  code  and  complexity. 

2.3.3  Field  Emitter  Solutions 

The  Designer  system  is  able  to  deal  with  systems  employing  either  electret  and/or 
conductor  elements.  Some  of  the  specific  solutions  are  described  below. 

The  CED  developed  a  set  of  equations  relating  electret  fibers  of  finite  length, 
conductive  drive  plates  of  finite  size,  and  conductive  sensor  plates  embedded  in  an 
infinite  ground  plane. 

The  most  significant  assumption  required  was  that  the  substrate  drivers  and  sensor 
plates  were  co-planar  with  an  infinite,  conductive  surface.  This  allowed  the  use  of 
simple  image  charges  to  model  the  fields  generated  by  a  charged  fiber  over  a 
conductive  plane.  Further  refinement  of  the  analysis  is  necessary  to  integrate 
dielectric  effects  which  are,  for  the  moment,  ignored. 

Plates  were  defined  as  conductive  elements  co-planar  to,  yet  electrically  isolated 
from,  the  infinite  ground  plane.  A  fiber  is  defined  as  a  strip  of  fixed  charges 
implanted  in  a  dielectric  material  (usually  a  free  floating  armature).  While  there  is 
not  in  actuality  a  group  of  "fibers"  per  se,  they  are  the  mathematical  abstraction  of 
what  is  left  of  the  armature  when  the  ignored  dielectric  structures  are  omitted. 

This  fact,  along  with  the  implicit  infinite  ground  'lane,  are  the  most  serious 
assumptions  made  in  the  analysis.  The  approximation  of  an  infinite  ground  plane  is 
valid  since  most  of  the  systems  which  we  have  postulated  have  relatively  large 
ground  plane  areas  surrounding  the  substrates.  The  problem  of  dielectrics,  either 
around  the  fibers,  or  instead  Gf  the  ground  plane,  is  much  more  serious  and  is  one 
which  is  currently  under  inve  .  Jgation. 

2.3.3. 1  Fiber-Fiber  Interactions 

For  the  geometries  which  have  been  explored  so  far,  the  only  fiber-to-fiber 
interactions  were  between  a  fiber  and  its  image  generated  by  the  ground  plane. 
This  is  because  of  the  assumption  that  the  armature  would  act  as  a  rigid  body 
and  therefore  all  internal  forces  and  torques  would  cancel  out.  Using  this  fact  to 
generate  the  equations  for  a  fiber  and  its  image  charge  would  h?ve  simplified  the 
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equations.  However,  in  order  to  maintain  the  generality  of  the  program,  it  wa^ 
explicitly  ignored  and  the  equations  relating  any  two  fibers  with  six  degrees  of 
freedom  between  them  were  used.  This  allowed  the  calculation  of  image 
attractions  by  defining  the  second  fiber  as  the  image  of  the  first,  as  well  as 
systems  with  multiple  armatures. 

The  E-field  due  to  the  first  fiber  can  then  be  integrated  over  the  three  space 
definition  of  the  second  fiber  to  produce  the  forces  and  torques  on  the  second 
fiber  due  to  the  first.  These  equations  are  quite  straightforward  and  easily 
derivable. 

These  integrals,  quantifying  forces  and  torques,  were  then  solved  using 
numerical  integration.  The  integrations  were  performed  using  a  'cautious 
adaptive'  Romberg  extrapolation  technique  as  implemented  by  the  DCADRE 
subroutine  of  the  IMSL  library  .  The  major  weakness  of  the  DCADRE  routine 
is  for  periodic  functions  with  several  cycles  within  a  single  set  of  limits.  The 
form  and  solutions  to  the  integral  equations  presented  here  do  not  fall  into  this 
type  of  equation  and  the  error  checking  within  DCADRE  did  not  indicate  any 
problems. 

2. 3. 3. 2  Fiber  -  Plate  Interactions 

The  fiber  to  plate  interactions  are  significantly  more  complex.  This  is  to  be 
expected  with  the  addition  of  a  conductive  plate.  In  a  manner  completely 
analogous  to  the  fiber-fiber  interactions,  a  coordinate  system  can  be  defined 
relative  to  the  fiber  orientation.  By  integrating  over  this  coordinate  system,  the 
forces  and  torques  imparted  to  the  fiber  can  be  determined.  These  resultant 
integrals  can  then  be  solved  using  the  identical  numerical  algorithms  as  the 
fiber-fiber  case. 

2. 3. 3. 3  Sensor  Interactions 

The  sensor  interactions  were  similar  to  the  fiber-fiber  calculations.  By 
integrating  the  normal  E-field  over  a  rectangular  plate  rather  then  a  second  fiber, 
the  charge  induced  in  the  rectangular  gate  region  of  an  FET  can  be  calculated. 
Because  this  is  an  integration  over  a  surface  rather  than  a  line,  the  IMSL  routine 
DBLIN  was  used.  This  routine  in  essence  implements  a  two-variable 
integration  using  DCADRE  in  both  (orthogonal)  directions. 

Since  all  of  the  equations  assumed  the  fiber  to  be  oriented  along  the  Z-axis,  and 
centered  at  (0,0,0),  the  sensor-fiber  system  was  transformed  to  the  proper  fiber 
alignment.  This  caused  the  sensor  plate  to  be  at  an  arbitrary  orientation.  Indeed, 
the  hardest  part  of  the  sensor  interaction  code  was  performing  the  geometric 
manipulations  required  to  obtain  two  orthogonal  basis  vectors  for  the  integration 
over  the  sensor  plate.  This  orthogonal  basis  was  related  to  the  width  and  height 
of  the  sensor  plate  after  the  transformation. 

2. 3. 3. 4  Numerical  Integration  Issues 

The  IMSL  subroutines  DCADRE  and  DBLIN  both  use  'cautious  adaptive' 
Romberg  integration  techniques.  A  search  for  'faster'  or  'better'  algorithms 
showed  that  these  two  routines  offer  the  best  compromise  between 
computational  efficiency  and  resolution.  There  might  be  some  small  advantage 
in  writing  custom  numerical  integration  routines,  but  the  gains  would  be  small. 
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2.3.4  The  model  definition  programs 

There  are  two  types  of  model  generation  programs  -  one  for  the  placement  of  drive 
elements  (electrets  and  plates)  and  the  other  for  placement  of  sensors.  In  addition, 
both  radially  symmetric,  helical,  and  arbitrary  layouts  are  supported.  The  only 
restrictions  are  those  imposed  by  the  underlying  equations  as  discussed  previously. 

The  program  first  places  the  drive  plates.  The  dimensions  of  a  plate  are  requested 
and  the  program  allows  the  placement  of  an  arbitrary  number  of  such  plates.  The 
system  then  allows  the  designer  to  choose  a  different  set  of  plate  sizes  and 
subsequently  place  them.  The  system  allows  the  designer  to  specify  both  position 
and  orientation  during  the  placement  phase.  There  is  at  present  an  upper  limit  of  50 
plates  (of  various  sizes)  but  this  is  an  arbitrary  limit  and  can  easily  be  increased. 

After  all  of  the  plates  are  defined,  the  system  moves  on  to  the  electret  filaments.  In 
a  manner  identical  to  the  placement  of  the  plates,  the  designer  can  specify  the 
electret  filaments  (in  length,  position,  and  orientation)  until  the  design  is  complete. 

While  the  size  specifications  are  keyboard  driven,  the  placement  and  orientation  are 
mouse  driven  with  interactive  graphic  editing.  This  provides  a  simple  and  visually 
interactive  method  of  placement  (see  Fig.  2.16).  The  mouse  interactions  are  much 
more  natural  and  assist  the  designer  to  focus  on  the  design  in  progress  and  not  the 
keyboard.  A  future  enhancement  will  allow  dynamic  modification  of  the  sizes  via 
the  mouse.  When  the  design  is  completed,  a  name  is  specified  and  the  saved  model 
can  be  used  as  input  for  any  of  the  analysis  routines. 

2.3.5  The  mapping-display  manager 

Similar  to  the  static  display  manager,  the  mapping  display  manager  handles  the  user 
interface  to  several  map-like  display  formats.  The  mapping  display  manager 
supports  three  separate  display  formats: 

1)  two  dimensional,  single  linear  triad  field  arrow  maps 

2)  sensor  sensitivity  maps 

3)  three  dimensional  field  arrow  maps  of  full  disk  systems. 

After  selecting  the  display  format,  the  mapping  display  manager  dispatches  to  the 
appropriate  program  to  perform  the  simulation  and  display.  Several  of  these 
programs  support  saving/restoring  calculated  information  to  speed  demonstrations 
and  comparisons  to  previous  systems. 

2. 3. 5.1  Two  Dimensional  Field  Maps 

Color  was  also  an  integral  part  of  the  two  dimensional  field  maps  which  were 
generated  for  examining  fields  around  drive  strips  similar  to  those  modelled  in 
the  field  validation  and  contour  control  law  simulations.  After  noticing  'null 
control  regions'  using  the  techniques  of  the  previous  section,  maps  of  the  force 
field  in  the  null  control  regions  were  requested  to  help  explain  the  formations  of 
the  null  control  regions. 

These  simulations  used  the  general  methods  discussed  in  the  previous  chapter 
for  fiber-fiber  and  fiber-conductive  plate  interactions.  A  very  small  test  charge 
(basically  a  very  short  length  armature  fiber)  was  used  to  test  the  force  field 
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generated  by  a  pair  of  plates  while  constrained  to  be  in  the  Y-Z  plane.  This 
constraint  limited  the  simulation  to  two  dimensions  and  was  a  first  step  towards 
the  full  three  dimensional  visualizations  of  the  next  section. 

In  order  to  obtain  more  information  display  per  simulation  run,  three  windows 
were  generated: 

1 )  Control  Force  (plate  to  fiber) 

2)  Image  Force  (fiber  to  fiber  image) 

3)  Total  Force  (sum  of  the  above). 

The  control  force  was  the  force  on  the  armature  due  to  the  voltages  set  on  the 
plates  (input  to  program)  and  the  charge  density  of  the  fiber  (also  input  to  the 
program).  The  image  force  was  the  attractive  force  on  the  armature  due  to  the 
image  charge  of  the  armature  across  the  conductive  plane  surrounding  the  drive 
strips  at  z=0.  And  the  total  force  was  simply  the  sum  of  two  previous  results. 

Fig.  2.17  shows  a  sample  run  of  a  simple  triad  simulation.  The  results  were 
displayed  using  arrows.  The  centroid  of  each  arrow  is  placed  over  the  point 
were  that  arrow's  force  calculations  were  performed  and  the  direction  of  the 
arrow  indicates  the  direction  of  the  force  field.  While  this  alone  would  have 
been  useful,  what  made  these  types  of  graphs  more  useful  was,  as  usual,  the 
addition  of  color.  Using  green  to  indicate  weak  forces  and  bright  red  to  indicate 
strong  fields  simplified  interpretation  of  the  force  fields. 

2. 3. 5.2  Feature  Extraction 

Not  only  is  color  used  to  indicate  the  magnitude  of  the  field  vector,  but 
'interesting'  regions  can  be  highlighted  by  changing  the  background  color.  In  a 
process  known  as  feature  extraction,  the  two  dimension*!  region  «.a,.  be  scanned 
for  areas  with  interesting  properties.  In  this  initial  demonstration,  two  features 
were  extracted.  This  first  interesting  feature  was  defined  as  a  region  where  the 
field  magnitude  was  changing  rapidly.  These  areas  correspond  to  those  areas 
where  there  are  high  field  gradients  and  indicate  areas  of  the  field  with 
significant  texture. 

The  second  feature  extracted  was  large  changes  in  the  vector  direction  of  the 
force  field.  Areas  where  the  field  direction  was  rapidly  changing  can  indicate 
stable  points  (if  all  of  the  arrows  point  together),  unstable  points  (if  all  of  the 
arrows  are  end  to  end)  or  meta-stable  points.  Fig.  2.18  shows  the  same  force 
fields  as  before  but  now  with  highlighted  'interesting  regions'. 

By  interactively  allowing  the  designer  to  specify  the  thresholds  of  the  feature 
extraction  algorithm,  the  designer  can  attempt  to  pinpoint  regions  of  interest  for 
further  concentrated  analysis.  A  third  technique  is  being  developed  to  indicate 
capture  and  escape  regions.  By  appropriate  color  coding  of  the  background,  the 
final  position  of  a  fiber  released  at  that  point  will  be  indicated.  As  an  example, 
if  it  is  released  within  a  limit  cycle  of  the  set  point,  the  background  will  be 
green.  If  it  will  eventually  crash  down  into  the  substrate,  it  will  have  a 
background  of  red.  Various  other  colors  will  indicate  that  it  will  escape  to  the 
left,  right,  vertically,  etc. 

Initially,  this  was  conceived  as  a  'connect  the  arrows’  approach.  However,  this 
would  be  mathematically  valid  if,  and  only  if,  each  arrow  had  the  same 
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magnitude.  Since  they  do  not,  one  would  have  to  perform  a  vector  addition  of 
the  arrows,  taking  magnitude  into  account.  At  this  point,  the  algorithm  requires 
more  intelligence  than  simply  examining  the  geometry  of  the  arrows.  The 
correct  approach  actually  is  to  release  the  armature  at  the  position  to  be  sampled, 
and  then  observe  the  trajectory  over  time.  When  the  substrate  leaves  the  limits 
of  the  two  dimensional  analysis  area,  the  exit  point  is  determined  and  indicates 
the  escape  regions. 

What  makes  this  type  of  analysis  difficult  is  determining  how  long  to  run  a 
given  simulation.  If  the  armature  quickly  escapes,  it  is  easy  to  establish  the 
escape  regions.  But  how  do  you  differentiate  between  a  true  limit  cycle  and  an 
escape  trajectory  which  did  not  escape  simply  because  the  simulation  ended  too 
soon.  To  make  matters  worse,  each  start  point  requires  a  time-course  simulation 
with  full  dynamic  effects  at  a  considerable  expense  of  compute  time. 

2. 3. 5. 3  Sensor  Sensitivity  Displays 

In  many  of  the  systems  under  consideration,  the  placement  of  sensors  is  not  at 
all  obvious.  The  sensors  need  to  be  close  to  the  sensed  elements  to  allow 
acceptable  signal  to  noise  ratios,  but  the  interference  from  the  related  drive 
elements  might  indicate  a  more  distant  placement.  To  provide  some  insight  into 
this  problem,  the  designer  can  create  sensor  sensitivity  displays. 

These  displays  show  the  regions  where  sensors  can  be  used  most  effectively. 
By  displaying  a  grid  of  pseudo-sensors  with  the  color  and  intensity  conveying 
magnitude  and  sign  of  the  integrated  field,  regions  which  are  'good'  for  sensors 
can  be  quickly  identified.  An  example  of  this  type  of  display  is  in  Fig.  2.19. 

A  derivative  of  this  analysis  would  be  to  perturb  the  armature  and  obtain  maps 
which  indicated  sensitivity  to  motion  rather  then  simply  static  fields.  This  kind 
of  analysis  could  be  automated  in  order  to  optimally  place  sensors  given  a 
specified  model  and  operating  constraints. 

2. 3. 5. 4  Three  dimensional  Force  Field  Displays 

The  natural  outgrowth  of  the  previous  section  was  the  three  dimensional  field 
arrow  concept  shown  in  Fig.  2.20.  As  before  the  orientation  of  the  arrows 
indicates  the  direction  of  the  field  and  the  color/shade  is  a  function  of  the 
magnitude  of  the  field  Because  of  the  need  for  indicating  perspective  and  three- 
dimensionality,  the  length  of  the  arrow  becomes  a  function  of  distance  to/from 
the  arrow.  This  is  done  by  using  a  fixed  length  arrow  and  letting  the 
perspective  model  shorten  the  more  distant  arrows.  This  approach  causes 
regions  with  intense  fields  to  stand  out,  which,  in  turn,  creates  an  inherent 
feature  extraction  capability. 

In  this  case,  to  get  the  smooth  color  scale  needed  to  blend  in  the  various  arrows, 
a  grey  shade  color  map  was  selected  rather  than  the  green,  blue,  red  model  used 
in  the  two  dimensional  case.  An  alternative  would  be  to  use  arrows  of  various 
lengths  depending  on  the  magnitude,  and  using  the  shading  to  indicate  distance 
(bright  foreground,  dim  background).  This  approach  has  the  advantage  of 
highlighting  the  immediate  surroundings.  It  also  has  the  disadvantage  of  having 
arrows  possibly  overlapping  each  other's  'character'  cell.  This  effect  is  visually 
bothersome  and  causes  significant  viewing  problems. 
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The  best  solution,  which  has  been  partially  explored,  is  to  generate  two  stereo 
image  pairs  and  use  this  technique  to  further  the  perception  of  three 
dimensionality.  There  have  been  several  references  to  systems  that  can  display 
stereo  images  directly,  though  none  of  these  systems  were  available.  Given  the 
modularity  of  the  programs,  as  well  as  the  simplicity  of  using  the  three 
dimensional  graphics  devices,  the  adaption  of  the  Design  System  to  true  three 
dimensional  would  be  relatively  simple. 

While  no  work  with  stereo  image  pairs  was  actually  done,  the  code  for  the  three 
dimensional  perspective  transformations  is  completely  general,  and  has  the 
capability  built  in  to  do  stereo  image  pairs.  These  could  then  be  viewed  using 
standard  stereo  vision  glasses  such  as  used  currently  for  three  dimensional 
contour  maps.  In  fact,  the  major  enhancement  to  existing  simulation  and 
modelling  capabilities  offered  by  this  portion  of  the  work  was  the  development 
of  the  three  dimensional  perspective  routines. 

The  routines  allow  the  specification  of  many  viewing  parameters  (such  as  center 
of  view,  angle  of  view,  eye  location,  etc.)  and  generate  a  single  transformation 
matrix  embodying  the  entire  viewing  transform.  This  allows  the  images  to  be 
calculated  and  stored  in  full  X,Y,Z  coordinates,  and  to  reduce  the  need  to 
transform  them  until  just  before  they  are  to  be 
displayed. 

This  method  also  allows  the  viewing  transformation  to  be  modified  via  a  menu 
interaction  with  the  designer,  and  the  image  quickly  redrawn.  By  allowing  such 
immediate  redisplay,  the  designer  can  'fly'  around  within  the  system,  and 
observer  the  fields  from  various  angles.  While  quick,  the  redisplay  takes  on  the 
order  of  5-10  seconds  and  consequently  does  not  lend  itself  to  animation. 

The  limiting  factor  in  attempting  to  animate  this  type  of  display  is  the  slowness 
of  the  software  driven  redisplay.  An  intelligent  display  system  with  internal 
three  dimensional  translation/rotation  hardware,  however,  would  very 
effectively  animate  this  type  of  display. 

2.3.6  The  Dynamic  Display  Manager 

The  dynamic  display  manager  is  the  high  level  interface  to  the  dynamic  display 
routines  As  before,  it  will  prompt  for  the  model  name  and,  if  necessary,  the  initial 
conditions,  as  well  as  if  this  is  a  playback  or  a  recording  operation.  Fig.  2.21 
shows  an  example  of  the  initialization  screen  and  a  snapshot  of  the  dynamic 
simulation  image. 

2.3.6. 1  What  are  animated  displays? 

Simply  speaking,  animated  displays  are  ones  that  change  over  time.  A  critical 
feature  is  that  they  change  over  time  IN  THE  SAME  MANNER  AS  THE 
MODELLED  SYSTEM.  A  animated  display  also  requires  good  animation  so 
that  the  effect  is  of  a  smoothly  changing  image  rather  than  3  jerky  series  of  static 
'snapshots'  of  the  system. 

In  order  to  produce  good  animated  effects,  three  issues 
must  be  addressed: 


1)  solid  modelling, 
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2)  simulation  of  system  dynamics  and  the  relationship  between 
wall-clock  time  and  simulation  time 

3)  computational  capabilities. 

2. 3. 6. 2  Solid  Modelling 

The  purpose  of  an  animated  display  is  to  give  the  designer  as  realistic  a  view  of 
the  system  as  possible.  This  allows  nuances  such  as  shading,  perspective  and 
color  to  carry  information  about  the  state  of  the  system  which  would  be  hard  to 
convey  otherwise.  It  is  much  easier  to  observe  tilt  angle  in  a  solid,  shaded 
image  of  a  disk  then  in  a  simple  wireframe  drawing.  By  using  these  types  of 
display  techniques,  the  image  can  be  made  to  reproduce  inherent  visual  cues 
quite  realistically. 

To  produce  these  highlights  and  visual  cues,  an  entire  three  dimensional 
perspective  and  solid/shaded  display  package  was  developed.  This  package  can 
generate  perspective  images  which  consist  of  filled  polygons.  These  polygons 
are  scan-converted  via  an  edge  table  algorithm.  The  scan-converted  lines  are 
then  put  into  a  'pre-converted'  Z-buffer  for  hidden  line  removal. 

The  pre-conversion  is  done  to  allow  the  fixed  background  and  substrate  details 
to  be  scan-converted  once,  and  then  loaded  into  the  Z-buffer  before  scan 
converting  the  dynamic  information.  The  Z-buffer  is  then  transferred  to  the 
screen  in  a  single  bitmap  operation. 

The  package  offers  a  variety  of  shading  models  (Flat,  Gouraud,  and  Phong; 
and  lighting  parameters  (position,  color,  and  metallic  or  plastic  highlights). 
These  'extras'  consume  a  significant  part  of  the  computational  resources 
necessary  for  rendering  the  images.  However,  the  visual  cues  that  they 
generate  are  worth  the  time  and  effort. 

2.3.6. 3  'Real-time'  dynamics 

To  be  effective  as  a  tool  for  understanding  the  time  course  behavior  of  a 
proposed  system,  the  simulation  must  show  the  motion  of  the  armature  at  an 
appropriate  rate.  The  real  speeds  of  motion  for  a  micron  sized  system  would  be 
a  blur  to  a  designer.  Therefore,  the  'real-time'  display  must  in  reality  be  an 
'expanded'  time  scale. 

This  simulation  time  scale  must  be  such  that  the  motions  of  moving  objects  be 
fluid  and  easily  observable.  The  expansion  factor  must  not  be  so  great  as  to 
make  the  motions  too  slow  and  jerky,  however.  Empirically,  it  appears  that 
motions  on  the  order  of  a  single  sheet  of  paper  floating  to  the  ground  seem 
optimal. 

The  'expanded'  real-time  is  also  highly  dependent  on  the  capability  to  calculate 
the  solutions  of  the  dynamic  equations  of  motion  quickly.  In  general,  it  has 
been  observed  that  present  computational  resources  are  not  fast  enough  for 
optimal  results  when  the  additional  complexity  of  rendering  the  image  is  added. 

The  equations  of  motion  used  in  the  simulations  were  obtained  by  assuming 
constant  acceleration  over  a  given  simulation  cycle  time. 

While  these  are  approximations,  with  appropriately  small  time-steps,  the 
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accuracy  will  be  well  within  that  attainable  for  either  the  force/torque 
calculations.  A  casual  examination  of  the  "smoothness"  of  the  fields  indicates 
the  underlying  assumption  of  constant  acceleration  is  valid.  In  addition,  many 
of  the  systems  which  are  being  analyzed  are  designed  to  move  over  very  small 
ranges  of  motion  which  also  helps. 

2. 3. 6. 4  Computational  Restrictions 

Any  discussion  of  display  techniques  or  modelling  of  simulation  motion  reverts 
in  the  final  analysis  to  that  of  computational  capability.  The  more  powerful  the 
host  computer,  the  more  that  can  be  imaged,  and  the  better  the  simulation  can 
become.  While  the  facilities  used  for  this  research  were  adequate,  they  were  at 
best,  just  barely  adequate.  Because  of  this,  many  restrictions  were  imposed. 

The  first  was  that  only  flat  shading  techniques  could  be  used.  The  added 
enhancement  of  Gouraud  or  Phong  shading  simply  took  too  much  time  and 
increased  the  image  refresh  time  to  unacceptable  levels.  A  second  restriction 
was  that  only  limited  dynamic  equations  could  be  solved.  As  was  mentioned 
above,  though,  this  was  not  of  serious  concern  at  present,  but  might  be  in  the 
future. 

An  additional  major  restriction  was  that  the  dynamic  simulations  had  to  be 
performed  off-line,  and  then  redisplayed  in  simulation  time.  This  allowed  the 
visual  display  to  run  at  full  speed  without  having  to  time-share  with  the  dynamic 
simulation  itself.  This  allowed  the  generation  of  acceptable  animation  though  it 
too  was  marginal.  It  soon  became  obvious  that  high  quality  animation  would 
require  substantially  more  hardware. 

2. 3. 6. 5  Two  Degree-of-freedom  Triad-Fiber  Systems 

The  first  animated  displays  that  were  implemented  were  for  examining  the 
predicted  behavior  of  cantilevered-fiber  systems.  The  details  of  the  simulations 
and  comparison  to  reality  are  theses  unto  themselves  (see  Khanwilkar,  1987, 
and  Grace,  1987). 

While  the  system  modelling  code,  control  laws,  and  dynamics  effects  code  was 
generated  by  both  of  the  above  gentlemen,  the  basic  display  program  and 
techniques  were  generated  during  the  course  of  this  research  work.  The  display 
consisted  of  a  small  filled  circle  (representing  the  fiber  armature)  moving  in 
two-space  over  the  VS-11  display.  A  joystick  allowed  input  to  vary  several 
system  parameters  during  the  simulation  run  as  well  as  a  switch  input  used  to 
stop  and  start  the  simulation.  Fig.  2.22  shows  a  display  captured  from  the  VS- 
11  screen.  Note  the  (optional)  history  trail  of  the  trajectory  taken  by  the 
armature.  By  using  various  start  conditions  and  trajectory  colors  several 
trajectories  can  be  compared  simultaneously. 

The  VS-11  is  a  512  by  512  frame  buffer  with  associated  circuitry  to  interface  it 
directly  to  the  bus  of  our  PDP-1 1/44.  The  low  level  interface  is  very  simple  and 
allows  basically  only  the  execution  of  small  display  lists  into  bitmapped  image 
memory.  The  user  fills  a  buffer  holding  the  various  graphics  commands  and 
then  sends  it  to  the  VS- 1 1  for  execution. 

The  VS-1 1  then  executes  the  display  list  into  the  local  VS-11  image  memory. 
Since  it  discards  the  display  list  after  execution,  changing  the  displayed  image 
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via  display  list  modification  is  not  practical.  Therefore,  the  only  animation 
technique  which  works  is  to  erase  the  previous  images  and  write  new  ones. 

In  order  to  prevent  flickering,  the  time  when  there  is  no  image  displayed  must 
be  minimized.  This  required  sending  the  first  image  to  the  screen  outside  of  the 
simulation’s  main  loop.  Within  each  cycle  of  the  simulation,  a  display  list  was 
generated  to  erase  the  previous  image  using  the  system  state  information  from 
die  previous  cycle,  the  dynamic  portion  of  die  simulation  was  run  to  update  the 
system  state,  and  then  the  new  display  information  was  added  to  the  display 
list. 

When  both  the  erasure  instructions  and  the  new  drawing  instructions  were  in 
the  display  buffer,  it  was  sent  to  the  VS-11.  This  left  as  small  a  window  as 
possible  when  there  might  be  no  image  on  the  screen.  In  spite  of  this,  there  is 
still  some  observable  flicker.  But  in  general,  the  flicker  is  not  objectionable  nor 
does  it  interfere  with  the  interpretation  of  the  simulation  results. 

Since  the  simulation  was  modelled  after  a  real  running  system,  the  various 
geometrical  and  electrical  parameters  of  the  simulation  were  set  to  match  the  real 
system.  At  this  point  only  the  variables  associated  with  the  control  laws  were 
left.  By  varying  the  gains  and  frequencies  involved  in  the  control  laws  and 
observing  system  response,  the  optimal  settings  for  the  real  system  were 
determined. 

The  display  format  was  chosen  to  indicate  maximal  information.  Since  the 
armature  was  constrained  to  only  two  degrees  of  freedom,  these  were  used  as 
the  coordinates  of  the  display.  This  provided  an  'end-on'  view  of  the  armature 
which  was  similar  to  the  previously  generated  trajectory  information.  Since  the 
arrangement  of  the  modelled  systems  were  linear  triads,  the  trajectory 
information  could  also  be  stored  off-line  and  fed  into  static  display  formats  such 
as  the  'apparent'  surfaces  and  trajectory  maps. 

All  other  information,  such  as  simulation  time,  control  law  parameters,  and 
precise  position  information,  was  displayed  in  textual  form  on  the  upper  portion 
of  the  display.  This  information  was  chosen  to  be  textual  since  it  was  useful  for 
identifying  parameters  before/after  a  simulation  run,  or  special  'interesting' 
positions  after  a  simulation.  Since  it  was  not  referred  to  explicitly  during  a  run, 
the  difficulty  in  reading  some  of  the  values  without  some  effort  was  not 
considered  a  significant  drawback. 

During  the  course  of  working  with  the  actual  system,  several  coordinate 
positions  were  found  which  had  the  property  of  being  uncontrollable.  During 
several  simulation  runs,  these  points  were  precisely  identified.  Observation  of 
the  voltages  on  the  drive  strips  as  the  armature  passed  through  these  points 
showed  peculiar  patterns  of  voltage  sign  reversal.  This  observation  led  to 
several  discoveries  about  the  controllability  of  the  armature  at  certain  positions. 
This,  in  turn,  precipitated  modifications  of  the  control  systems  and  also  to 
explore  control  law  efficiency  as  was  previously  discussed. 

In  general,  it  was  found  that  the  correlation  between  simulation  and  reality  was 
quite  good.  An  additional  feature  that  assisted  in  this  comparison  was  the 
compatibility  between  the  simulation  computer  and  the  real  control  computer. 
Since  both  were  PDP-ll's,  the  same  code  which  was  analyzed  during  the 
simulations  could  be  directly  downloaded  to  the  real  system  controller  without 
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recoding  or  modifying  the  code. 


2. 3. 6. 6  Six  degree-of-freedom  Triad-Fiber  Systems 

Future  systems  would  not  necessarily  be  tethered,  cantilevered  systems  similar 
to  the  two  degree  of  freedom  system  just  discussed.  Therefore,  a  simulation  of 
a  six  degree  of  freedom  disk  system  was  performed.  The  simulation  used  the 
equations  developed  in  section  X.X  to  derive  the  forces,  torques  and  sensor 
outputs  for  a  disk  system  at  any  given  six-space  point.  These  results  were  then 
passed  through  a  controller  subroutine  which  attempted  to  stabilize  the 
armature. 

The  controller  would  modify  the  voltages  on  the  substrate’s  drive  plates.  The 
simulation  would  then  move  the  armature  as  determined  by  the  dynamics  of  the 
system  and  redisplay  the  image.  At  this  point,  the  cycle  would  repeat  with  the 
calculation  of  forces  and  torques.  Fig.  2.23  shows  an  artist's  conception  of 
such  a  system. 

The  only  variables  which  could  be  easily  set  before  a  simulation  run  were  the 
six-degree-of-freedom’s  of  the  armature  start  position  and  the  armature  set- 
point.  In  order  to  change  geometrical  configuration  information,  one  of  the 
model  definition  programs  would  need  to  be  run.  The  voltages  on  the  plates 
would  normally  be  an  internal  variable  which  varied  in  response  to  the  position 
of  the  armature  according  to  the  control  law  being  used.  In  the  present 
implementation  the  control  law  algorithm  was  a  null  subroutine  and  therefore 
the  plate  voltages  were  fixed. 

This  was  done  since  the  actual  control  law  subroutine  was  not  implemented. 
Since  the  major  focus  of  this  research  was  on  how  to  display  dynamic 
information,  not  on  what  the  dynamic  information  was,  the  controller  issue  was 
not  addressed.  In  the  two  degree  of  freedom  case,  the  controller  complexity  was 
such  that  the  solution  was  on  the  order  of  a  separate  Masters  degree  project.  In 
the  six  degree  of  freedom  case,  the  analysis  is  significantly  more  complicated 
and  was  well  beyond  the  scope  of  this  work  or  anything  less  than  a  Doctorate 
level  project. 

Therefore,  the  initial  and  desired  state  of  the  armature  were  input,  and  the 
simulation  clock  started.  Since  gravity  was  included,  the  armature  would,  in 
general,  fall  slowly  towards  the  substrate.  The  trajectory  taken  was  a  function 
of  the  voltage  on  the  fixed  voltage  plates. 

Initially,  the  only  display  format  was  an  image  of  the  system  using  shaded  solid 
models  of  the  various  structures.  By  using  three  dimensional  perspective 
modelling  techniques  along  with  the  shaded  solid  models,  the  exact  six-degree- 
of-freedom  orientation  of  the  armature  could  easily  be  observed. 

As  the  simulation  progressed,  the  screen  image  was  updated  each  cycle  time  to 
provide  animation.  Unfortunately,  the  time  required  for  computing  the  forces, 
torques,  etc.  was  too  long  and  the  images  became  jumpy  rather  than  fluid. 

At  this  point,  the  programs  were  modified  to  either  calculate  the  trajectory 
information  and  save  it  in  a  data  file,  or  to  simply  playback  a  stored  data  file. 
While  this  significantly  enhanced  the  speed  of  the  animated  display,  the  result  is 
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still  a  little  jerky  due  to  the  rendering  times  required.  The  actual  display  time  for 
a  given  Z-buffer  was  quite  fast  since  the  VAXstation  II/GPX  allows  direct 
bitmap  operations.  However,  the  scan-line  conversion  was  still  too  slow  for 
clean  animation.  Since  the  rendering  algorithms  were  already  optimized  for  the 
specific  types  of  displays  used,  any  increase  in  speed  will  have  to  involve  new 
hardware.  Either  a  faster  computer,  or  an  enhanced  graphics  display  which 
does  the  rendering  in  hardware,  will  be  necessary. 

2.3.7  The  Static-display  Manager 

The  static  display  manager  handles  the  selection  of  the  output  format  as  well  as  the 
input  of  model  name  and  initial  conditions  and  parameters. 

After  entering  the  appropriate  information,  the  static  display  manager  calls  the 
necessary  routines  to  present  an  image  of  the  system  as  well  as  a  tabulation  of  the 
forces,  torques  and  sensor  outputs  at  the  desired  (fixed)  configuration.  The  system 
then  prompts  for  a  new  position  and  continues  indefinitely  in  this  manner. 

The  second  output  display  format  is  a  static  display  of  the  trajectory  information 
using  the  same  graphics  routines  discussed  previously.  This  was  done  to  provide  a 
relatively  standard  set  of  display  formats  and  prevent  a  proliferation  of  wildly 
varying  graphs  and  plots.  Some  modifications  of  the  display  were  required  since 
there  were  now  six  degrees  of  freedom  instead  of  two. 

Another  static  display  which  was  generated  with  no  additional  work  was  a  simple 
image  of  the  structures  with  a  tabular  indication  of  torques,  forces,  and  sensor 
outputs.  This  display  used  the  exact  same  modules  as  the  dynamic  display  but  it 
cycled  through  the  display  loop  only  once  and  textually  showed  the  results  instead 
of  passing  them  to  the  controller  subroutine.  This  allowed  a  static  display  which 
could  be  easily  integrated  into  the  Design  System  with  little  overhead,  yet 
maintained  a  standard  display  format. 
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Fig.  2.1.  Four  types  of  electrets  fabricated  during  the  Microfield  project,  (a)  Sheet 
electret  patch  of  charge.  While  the  compensated  version  is  shown  here,  we  have  also 
fabricated  uncompensated  sheets,  (b)  Patterned  electret  (compensated  and 
uncompensated),  (c)  Patterned  electret  on  polymer-coated  silicon,  (d)  Uncompensated 
fiber  electret  formed  by  implanting  charge  in  polymer  coating  applied  to  140-rrucron  quartz 
fiber.  [1985] 
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Fig.  2.2.  Illustration  of  difference  between  compensated  and  uncompensated  sheet 
electrets.  (a)  Compensated  electret  with  grounded  metal  coating  on  back.  Positive  charges 
in  metal  compensate  implanted  negative  charges  in  polymer,  producing  strong  internal  E- 
field.  (b)  Uncompensated  electret,  consisting  of  nonmetallized  polymer  with  net  implanted 
charge,  producing  a  strong  external  E- field.  [1984] 
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Fig.  2.3.  Compensated  patterned  sheet  electret  fabricated  by  implanting  connected 
parallel  lines  of  charge  in  a  sheet  of  1-mil  Teflon®  aluminized  on  back.  Black  Xerox  toner, 

which  adheres  to  negatively  charged  areas  of  the  Teflon®,  show  the  pattern  of  implanted 
charge.  [1985] 


Fig.  2.4.  Profiles  of  a  patterned  electret  showing  "fang"  effect,  obtained  by  scanning 
FET  field  sensor  over  the  surface  of  electret  at  three  different  distances  (a)  100  microns; 
(b)  50  microns;  (c)  25  microns.  [1986] 
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Fig.  2.5.  Profiles  of  five  different  patterned  electrets,  showing  effect  of  dose  on  '  fang" 
effect.  Dose  is  calculated  by  multiplying  electron  beam  current  by  sweep  rate  of  beam 
across  surface  during  fabrication.  This  results  in  a  linear  dose  figure  given  in  terms  of 
incident  charge  per  unit  length,  (a)  1.8  nC/cm.  (b)  0.6  nC/cm.  (c)  0.5  nC/cm.  (d)  0.4 
nC/cm,  achieved  with  a  beam  current  of  0.15  nA  and  a  sweep  rate  of  0.33  cm/sec.  (e)  0.3 
nC/cm.,  achieved  with  a  beam  current  of  0.03  nA  and  a  sweep  rate  of  0.08  cm/sec.  Note 
that,  for  a  given  dose,  increasing  the  sweep  rate  produces  a  larger  attenuation  of  the  "fang" 
effect  than  does  lowering  the  beam  current.  [1986] 
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Fig.  2.6.  Plot  of  temperature-compensated  source-<irain  current  in  a  small  K-FET  device 
(2  microns  by  40  microns)  as  a  function  of  position,  obtained  by  passing  the  FET  over  the 
surface  of  a  patterned  compensated  sheet  electret  at  a  height  of  approximately  100  microns. 
[1986] 


Fig.  2.7.  Photomicrograph  of  FET  chip  of  "K”  design.  Four  different  sizes  o floating 
gate  FETs  are  located  dose  to  the  end  of  the  die.  The  smallest  is  approximately  2  micron 
t>v  40  microns.  Corresponding  reference  devices  (not  visible  in  this  field  of  vieu  /  an 
provided  on  the  same  chip.  1 1986) 


Fig.  2.8.  Photomicrograph  of  FET  field  sensor  of  "SG’  design.  Two  parallel  floatin 
gates  are  visible,  separated  by  about  18  microns.  The  wide  metallic  structures  are  met. 
conductors  used  to  address  the  drains  and  sources  of  the  FETs.  [  1986] 
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Fig.  2.9.  Photomicrograph  of  FET  chip  of  ”Q"  design.  Four  identical  pairs  of  annular 
floating-gate  FETs  are  visible.  Diameter  of  each  gate  is  approximately  50  microns,  f  19861 
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Fig.  2.10.  Illustration  of  six  degree-of-freedom  sensor  design  concept.  Four  pairs  of 
FETs  are  used  to  detect  small  motions  of  four  field  emitters,  each  of  which  is  located  over 
one  pair  of  gates.  [1982J 
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Fig.  2.11.  The  basic  electrostatic  configuration  for  position  sensing  of  conductive  strips 
with  an  FET.  [1982] 


Fig.  2.12.  The  electrostatic  problem  which  generates  the  Green  function  for  the 
configuration  in  Fig.  2.1 1.  [1985] 


MicroFiclds  Final  Report 


73 


Fig.  2.13.  A  line  charce  in  an  environment  with  parallel  conductors  and  dielectric 
boundaries.  [1985] 


Fig.  2.14.  Geometry  for  studying  the  sensing  of  a  line  charge  by  an  annular  FET  gate. 
|1985] 
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Fig.  2.15.  Organization  of  the  Design  System.  [1986] 


Fig.  2.16.  Placement  phase  of  model  definition.  [1987] 


Fig.  2.17.  Force  fields  due  to  a  simple  triad.  [1987] 


Fig.  2.18.  Field  arrows  with  extracted  features.  [1987] 
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Fig.  2.19.  Sensor  sensitivity  map.  [1987] 
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Fig.  2.22.  Example  of  the  VS-11  display  format.  [1986] 


Fig.  2.23.  Six-degree-of-freedom  Triad-fiber  disk  System.  [1987] 
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3.  A  Small  Cantilevered  Optical  Fiber  Servo  System: 
SCOFSS  in  Detail 

3.1.  Introduction 

Small  or  micro  electro-mechanical  systems  (MEMS)  can  often  deliver  numerous 
performance  advantages  over  their  larger- scale  counterparts.  These  typically  include  higher 
bandwidth  of  operation,  higher  energy  densities,  and  the  dominance  of  controllable 
interactive  forces  over  passive  (e.g.,  gravity,  elastic,  and  image)  forces.  In  particular,  the 
Center  for  Engineering  Design  is  interested  in  controlling  the  movement  of  micron- scale 
objects  and  in  sensing  their  positions.  As  an  intermediate-  (or  mini-)  scale  prototype  for 
testing  control  and  sensing  principles,  we  have  constructed  a  Small  Cantilevered  Optical 
Fiber  Servo  System  (SCOFSS)  as  a  precursor  to  smaller-scaled  systems  and  to  fully 
levitated  devices  (i.e.,  no  elastic  forces).  This  system  has  proven  useful  for  experimentally 
verifying  the  modelling  of  different  force  contributions,  such  as  image  forces  in 
conductors,  gravitational  forces,  controller  forces,  and  elastic  forces.  This  system  also 
involves  fabrication  issues  not  unlike  those  expected  for  other  planned  smaller  systems, 
including  implantation  of  electrons  into  polymers  (which  are  typically  not  charge- 
compensated  by  metalization),  servo  control  of  a  small  object,  micro-chip  design  and 
fabrication,  and  accurate  sensing  of  a  small  object.  Applications  of  related  devices  might  be 
in  simple  acceleration  sensors,  spatial  light  modulators,  or  actuators  [Jacobsen  and  Wood, 
1984;  Jacobsen,  1987]. 


3.2.  Experimental  Apparatus 

The  SCOFSS  system  is  illustrated  in  Fig.  3.1.  The  cantilevered  beam  is  a  quartz  Fiber¬ 
optic  of  length  L  =  5  cm  and  of  radius  a  =  70  micron  (including  cladding),  mounted  in  a 
chuck,  with  a  5  mW  laser  source  applied  at  the  proximal  (supported)  end.  The  laser,  as  it 
spreads  from  the  beam  tip  with  a  cone  angle  of  30°,  projects  a  spot  approximately  2500 
microns  in  diameter  onto  a  Lateral  Effect  Photo  Diode  (LEPD),  whose  active  sensor  plane 
is  located  about  t  =  0.45  cm  beyond  the  beam  tip.  The  LEPD  responds  to  the  centroid  of 
the  spot,  with  about  one  micron  resolution. 

The  cantilever  beam  has  an  /  =  1.5  cm  region,  which  terminates  at  b  =  4.5  cm  from  the 
supported  side,  dip-coated  with  a  charge-retentive  polymeric  material  such  as  Teflon  (FEP 
or  TFE),  polycarbonate  (PC),  polystyrene  (PS),  or  polyimide  (PI).  The  coatings, 
preferably  PC  or  TFE  depositions  (PC  is  used  in  the  experiments  herein),  are  layered  up  to 
about  3-4  microns  in  thickness.  The  fiber  and  coating  are  then  irradiated  in  an  SEM  (at 
about  30  keV),  thus  forming  an  electret  region  over  length  /.  Effective  linear  charge 

densities  on  the  order  of  X  =  -0.1  nanoCoulombs/centimeter  (nC/cm)  are  typically 
achieved. 

The  charged  region  of  the  beam  is  then  positioned  over  a  set  of  parallel  conductive  driver 
strips,  each  of  width  w  =  192  microns  (on  c  =  200  micron  centers)  and  length  s  =  2  cm. 
Thus,  if  the  shallow  angles  of  arching  are  neglected  (see  Sec.  3.3.4)  and  only  small 
deflections  are  allowed,  then  the  system  can  be  considered  as  a  rod  with  two  degrees  of 
ireedom  (y,z)  to  be  controlled. 

The  other  elements  of  the  system  include  a  MINC  (DEC)  computer  which  implements  the 
control  law,  and  specifies  the  appropriate  potentials  (up  to  ±150  volts  by  high  impedance 
amplif  crs)  to  be  applied  to  the  driver  strips. 
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Fig.  3.1.  Schematic  of  SCOFSS  system,  with  laser,  hber-optic 
beam  (with  electret  region),  LEPD,  and  driver  plates.  The  M1NC 
computer,  and  the  amplifiers  which  drive  the  plates,  are  not 
illustrated.  [1985] 


3.3.  System  Forces 

The  forces  acting  on  the  fiber-optic  beam  can  be  categorized  as  "passive"  (uncontrolled) 
and  "active”  (controlled).  In  order  to  put  all  forces  on  a  common  footing,  we  consider  all 
vector  forces  to  be  the  forces  (or  loads)  which,  if  acting  uniformly  over  the  electreted 
region  of  the  filament,  would  produce  a  deflection  equivalent  to  the  actual  force  (or  load). 
The  individual  forces  are  detailed  below. 

3.3.1.  Passive  System  Forces 

Four  forces  affect  the  beam  when  the  conductive  driver  strips  are  grounded.  These  are: 

(1)  The  elastic  force  ( Fe )•'  This  is  a  restoring  force  pointing  towards  some 
equilibrium  reference  position  (this  need  not  be  the  unstrained  position),  For 
observed  beam  deflections  away  from  the  reference  position  (y0,zo),  there  will  be 
some  force  Fe  which  can  be  considered  to  be  acting  uniformly  over  the  electreted 
region  to  produce  that  deflection.  The  force  components  are  given  by  [Timoshenko 
and  Young,  1968;  Price,  1985] 

Fe(y)  =  ~(b*)-3(E*)a*)(V-Vo) 

a 

Fe(z)  =  —  (b*)-3(E*)(I*)(C-£o) 
a 

where  c*  =  -10.9  dynes,  a  =  0.25  is  the  loading  factor  (see  discussion  below),  b* 
=  b/(l  cm),  E*  is  the  beam  modulus  normalized  by  E  =  7.2  x  10^  1  dynes/cm^,  and 
I*  is  the  beam  area  moment  of  inertia  normalized  by  I  =  na^/4  =  1.89(10'^)  cnA 

Typical  values  are  b*  =  4.5,  I*  =  1,  and  E*  =  1.  Also,  y  and  £  are  the  horizontal 
(y)  and  vertical  (z)  displacements  of  the  beam  tip,  normalized  by  100  microns. 
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Normally,  Vo  and  £o  would  be  the  normalized  reference  coordinates  of  the  beam 
attachment  point  above  the  substrate.  However,  as  discussed  in  Sec.  3.3.4,  after 

tilting  the  beam,  Vo  and  £0  become  the  initial  equilibrium  location  of  the  point 
defined  by  b*. 

A  nominal  value  of  a  =  1/4  was  used  in  the  above  elastic  deflection  equation  in 
order  to  make  a  meaningful  force  sum  with  the  image  and  control  forces  (and  the 
"corrected”  gravity  force)  which  act  on  the  electret  which  covers  the  distal  third  of 

the  beam  (a  was  calculated  exactly  for  each  coated  fiber,  based  on  the  actual  length 
and  position  of  the  electret  region). 

(2)  The  image  force  (Fi):  This  is  an  attractive  force  between  the  electreted  region  (/) 
of  the  filament  and  the  conductive  base.  It  has  components  (using  the  thin  filament 
approximation) 


Fi(y)  =  0 
k*  ( X*)2(l *) 


where  k*  =  -9.0  dynes,  and  where  X*  is  the  linear  charge  density  normalized  by 
0.1  nC/cm,  and  l *  is  the  length  of  the  electron-charged  region  along  the  filament 
normalized  by  1  cm.  Typical  values  are  X*  =  -0.5  and  l*  =  1.5,  which  at  the  usual 
reference  height  of  £  =  5.2  (see  Sec.  3.4),  yield  an  image  force  of  0.65  dynes. 

Since  the  beam  is  tilted  (see  Sec.  3.3.4),  however,  the  initial  elastic  deflection  of  the 
beam,  due  to  both  the  image  force  at  the  reference  height  and  the  gravitational  force, 
is  nulled  out.  The  effective  (eff)  image  force  is  then  given  by 

Fff(z)  =  Fi(z)-F°(zo) 

where  it  is  noted  that  the  effective  image  force  is  zero  at  the  reference  height, 
negative  bet  w  it  and  positive  above  it. 

(3)  The  gravitational  force  (Fg):  This  force,  distributed  over  beam  length  L,  is 
normal  to  the  base  in  actual  use,  with  components 

Fg(y)  =  0 


Fg(z)  =  d*  (g*)(p*)(V*) 

where  d*  =  -2.0  dynes,  g*  is  the  acceleration  due  to  gravity  normalized  by  g  =  980 
cm/sec^,  p*  is  the  mass  density  normalized  by  p  =  2.65  grams/cm^  (quartz),  and 
V*  is  the  filament  volume  normalized  by  V  =  7t  sl2  L  =  7.7  (10*4)  cm^.  Typical 

values  are  g*  =  p*  =  V*  =  i.  The  weight  of  the  polymer  (PC  electret)  deposition  is 
neglected. 
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The  effective  gravity  force  of  the  beam,  i.e.,  that  weight  distributed  over  electret 
region  l  which  gives  the  same  static  deflection  at  b  as  the  beam  weight  distributed 
over  total  length  L,  is  calc  hated  to  be,  using  standard  beam  theory, 

Fgff  =  (0-6)  Fg(z). 

However,  it  is  noted  that  part  of  the  initial  tilt  setting  used  herein  (see  Sec.  3.3.4) 
cancels  out  the  elastic  (tisplacement  of  the  beam  due  to  gravity.  The  height- 
independent  gravity  term  is  thus  not  of  concern  in  the  operation  of  the  controller.  In 
the  general  controller  discussions  (which  could  presume  an  untilted  beam),  the 

eff 

gravity  force  vector  Fg  should  be  taken  to  be  the  effective  gravity  force  Fg  .  In  the 
tilted  beam  case,  at  the  reference  position,  we  can  set  Fg  =  Fj  =  Fe  =  0. 

(4)  The  dielectrophoretic  (DEP)  force  (Fdep):  The  DEP  force  [see  Jones  and  Kallio, 
1979]  can  be  divided  into  "image-related’rand  "control-related”  components  (see 
Sec.  3.3.2).  Both  types  were  analyzed  for  this  system,  and  were  found  to  be 
relatively  insignificant  at  the  size-scale  dealt  with  herein  (at  smaller  size-scales,  the 
DEP  forces  become  more  important).  Thus  they  are  not  included  in  the  modelling. 

3.3.2.  Active  System  Forces 

Two  force  types  (DEP  and  drivers)  have  components  affected  by  the  voltages  applied  to 
the  substrate  strips.  However,  as  indicated  above,  the  control-related  components  of 
the  DEP  foice  are  ignorable,  and  are  thus  not  discussed. 

(1)  The  driver  force  (Fd):  Application  of  non-zero  potentials  to  the  driver  strips 
adds  forces  (over  region  l)  to  the  beam  in  two  controllable  directions  (y,z).  The 
dependence  of  these  driver  forces  can  be  expressed  as 


Fd  -  fiX  Dj,  Xact) 


where  Di  are  the  set  of  Five  plate  voltages  (i=l,2,...5),  Xact  is  the  actual  2D 

position  of  the  electreted  region  of  the  Filament  with  respect  to  the  plates,  and  X  is 
the  linear  charge  density  on  the  Filament.  The  plates,  modelled  as  finite  in  width 
and  infinite  in  length,  are  considered  to  be  set  in  an  infinite  (grounded)  conductive 
plane.  A  derivation  of  the  details  of  the  above  relationship  is  not  presented  [Price, 
1985],  although  it  is  noted  that  the  forces  are  a  linear  homogeneous  function  of  the 
strip  potentials,  with  terms  of  the  form 


Fp)i  -  X  D]  f  (Xact) 


where  f  is  a  nonlinear  fur.:tion  of  position. 

3.3.3.  Stability  Considerations 

In  the  absence  of  elastic  forces,  there  is  no  set  of  fixed  plate  potentials  which  can  keep 
the  Filament  in  a  stable  "levitated"  state  [Eamshaw,  1842].  However,  ii.  principle,  a 
time-varying  set  of  plate  potentials  can  be  applied  to  suspend  the  filament  stably. 
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With  the  addition  of  inherently  stabilizing  elastic  forces  (i.e.,  a  cantilevered  system), 
passive  static  stability  is  achieved  for  certain  equilibrium  heights  above  the  conductive 
base.  Different  equilibrium  heights  are  obtained  by  different  sets  of  driver  potentials. 
For  beam  equilibrium  locations  sufficiently  close  to  the  plates,  the  system  becomes 
passively  unstable  wherein,  for  small  perturbations,  the  sum  of  image,  driver  and 
gravity  forces  dominates  the  elastic  force  and  the  beam  becomes  latched  down  to  the 
base  [Petersen,  1978]. 

3.3.4  Beam  Alignment 

The  alignment  of  the  beam  over  the  driver  strips  was  important  to  proper  SCOFSS 
operation.  Lateral  alignment,  such  that  the  beam  was  parallel  to  the  strips,  was 
performed  in  two  steps.  First,  coarse  alignment  was  achieved  using  an  optical 
microscope  looking  down  on  the  beam  and  strips.  Second,  sinusoidal  inputs  were 
applied  to  the  central  three  drive  plates.  Fine  lateral  alignment  was  achieved  when  the 
beam  no  longer  showed  an  oscillatory  lateral  component  of  displacement  on  the 
oscilloscope. 

Vertical  alignment  was  more  involved  than  the  lateral  alignment.  First,  the  droop  (at  b 
=  45,000  microns)  due  to  gravity  alone  is  about  200  microns,  with  a  tip  angle  (at  L)  of 
about  0.35°.  Second,  the  droop  (at  b)  due  to  the  image  force,  for  parameters  z  =  5.2 
(for  this  choice,  see  below),  1*  =  -0.5  and  /*  =  1.5,  is  about  140  microns.  The 
combined  droop  is  thus  about  340  microns.  To  compensate  this,  the  beam  was  then 
tilted  up  (like  a  fly  rod;  see  Fig.  3.1)  such  that  the  distal  portion  of  the  beam  was 
parallel  to  the  drive  strips.  This  procedure  essentially  nulls  out  the  static  elastic  force 
due  to  gravity  and  the  image  force  at  the  reference  position,  and  removes  any  parallax 
error  at  the  LHPD  for  the  equilibrium  reference  position. 

Vertical  alignment  then  proceeds  by  adjusting  the  beam  as  it  responds  to  sinusoidal 
inputs.  Specifically,  when  the  three  central  plates  are  grounded,  and  the  two  groups  of 
three  plates  on  either  side  of  the  central  group  are  driven  at  the  same  potential,  there  is  a 

null  force  point  at  yQ  =  0  and  zo  =  (3c/2)  V3  ~  520  microns,  which  serves  as  a  useful 
calibration  point.  To  use  this,  the  outer  two  groups  of  three  plates  are  driven 
sinusoidally  while  the  height  of  the  cantilevered  end  of  the  fiber  is  adjusted.  When  the 
amplitude  of  the  beam  responses  are  minimal  or  zero,  the  beam  height  is  taken  to  be  at 
the  theoretical  value  of  520  microns.  This  point  is  also  the  equilibrium  position  of  the 
beam  for  the  normal  five- plate  operation  described  below.  Because  the  field  lines 
generated  by  the  five  plate  control  configuration  are  different  than  those  produced  for 
the  nine  plate  calibration  process,  there  is  then  no  null  condition  at  (yo,zo)- 


3.4.  Closed  Loop  Servo  Controller 

The  closed  loop  control  of  the  beam  can  be  used  to  advantage  in  both  the  passively  stable 
and  passively  unstable  regions.  The  control  loop  can  be  broken  into  four  basic 
components:  the  state  estimator,  the  control  law,  the  driver,  and  the  plant  (which  can  be 
either  a  computer  model  or  the  actual  cantilevered  filament).  The  elements  are  shown 
schematically  in  Fig.  3.2. 
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Fig.  3.2.  Schematic  of  principal  elements  of  SCOFSS 
controller.  The  input  is  Xdes-  1°  the  real  (physical) 
case,  the  actual  system  position  (Xact)  is  fed  back  to  the 
controller  via  the  state  estimator  which  produces  a 
measured  (estimated)  position  and  velocity.  For  model 
simulations,  the  output  is  instead  acceleration,  which 
must  be  integrated  twice  to  get  the  position  state.  See 
text  for  details.  [1986] 

3.4.1.  State  Estimator 

The  state  estimator  in  the  SCOFSS  system  is  accomplished  via  an  LEPD  reading  of  a 
laser  spot  emanating  from  the  beam  tip.  The  output  of  the  LEPD  is  a  (y,z)  coordinate 
reading,  approximating  the  position  of  the  beam  tip,  which  is  then  transmitted  to  the  16- 
channel  A/D  conveners  of  a  computer  controller  (a  DEC  MINC).  The  LEPD  coordinate 
reading  is  then  corrected  for  beam  tip  angle  by 

yest  =  P  Ay  +  yo 

zest  =  p  Az  +  Zq 

where  yest  and  zest  estimate  the  2D  position  of  the  beam  at  point  b,  where  p  =  0.75 

(nominal  value;  calculated  exactly  for  each  fiber),  and  where  Ay  and  Az  are  the  LEPD- 
observed  displacements  away  from  the  reference  position. 

The  corrected  position  (Xest)  signal  is  then  differentiated  (via  analog  circuitry)  to 
produce  an  estimated  velocity,  then  run  through  a  first-order  filter,  before  going  to  the 
controller  where  it  is  used  in  the  damping  term  of  the  control  law  (Vesfi  see  Sec. 
3.4.2). 

In  earlier  SCOFSS  prototypes,  some  of  the  parallel  base  strips  were  used  as  sensors, 
with  their  potentials  free-floating,  thus  reflecting  the  proximity  of  the  charged  filament. 
However,  capacitive  coupling  to  the  sapphire  base  and  adjacent  polysilicon  driver 
strips,  and  low  signal-to-noise  ratios  plagued  this  approach.  Future  devices  will 
eliminate  the  LEPD,  and  instead  use  the  output  of  sensitive  FET  gates  set  beneath  the 
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filament  (Jacobsen,  et  al.,  1987).  For  such  devices  there  would  be  equations  for 
converting  FET  currents  into  estimates  of  electret  positions  (after  the  influence  of  the 
near-by  driver  strips  had  been  subtracted  out).  With  the  LEPD,  no  such  equations  are 
needed. 


3.4.2.  Control  Law 

The  control  law  is  a  two-degree-of-freedom  (2DF)  proportional  formulation,  with 
velocity  damping  (PD  controller  with  position  reference  only).  The  controller 
•  prescribes  a  desired  total  force  (Fdes)  to  be  applied  to  the  filament,  based  on  the 

position  of  the  beam  relative  to  the  desired  position,  and  on  the  actual  velocity  of  the 
beam  (see  the  vector  representation  in  Fig.  3.3). 

The  control  law  is,  in  vector  form,  given  by 


Fdes  -  Kp(Xcjes  -  Xest)  -  Kv  Vest 


where  Xdes  is  the  desired  set  point,  Xest  is  the  beam  position  as  determined  by  the 
LEPD  (corrected),  and  Vest  is  the  estimated  2D  velocity.  Kp  and  Kv  are,  in  general, 
diagonal  matrices  [2  by  2]  of  position  and  velocity  gains,  respectively.  In  the  actual 
implementation,  they  were  simply  scalars. 


Fig.  3.3.  Schematic  of  control  law  model.  The  fiber 
diameter  is  represented  by  the  circle  at  Xact  (roughly  the 
2D  location  of  the  electreted  region  of  the  fiber).  The 
fiber  has  a  static  equilibrium  at  Xref-  Fdes  has  two 
vector  components:  Fdes(x)  which  points  towards  the 
desired  position  (Xdes)-  ard  Fdes(y)  which  is  opposite 
in  direction  to  the  actual  velocity  Vact  at  the  tip.  In 
practice,  Xact  and  Vact  are  replaced  by  Xest  and  Vest 
respectively.  [1986] 


3.4.3.  Driver 

The  driver  takes  into  account  models  of  the  image  force  (Fj),  the  gravity  force  (Fg), 
and  elastic  force  (Fe),  then  calculates  the  potentials  to  be  applied  to  the  driver  strips  in 
the  base,  so  as  to  generate  a  total  net  vector  force  on  the  filament  equal  to  the  desired 
force  (Fdes)  determined  by  the  control  law.  Since  only  two-degrees  of  freedom  are 
being  controlled,  in  principle,  only  two  conductive  strips  are  needed  to  effect  movement 
in  all  directions.  However,  we  chose  to  actuate  the  filament  with  five  driver  strips. 
This  affords  more  force  and  versatility,  but  it  also  introduces  a  third-order  redundancy. 
By  constraining  the  central  three  strips  to  be  of  equivalent  potential,  the  level  of 
redundancy  is  reduced  to  one.  This  is  resolved  (using  a  Lagrange  multiplier  approach) 
by  minimizing  the  sum  of  squares  of  the  outer  two  independent  potentials  (D]  and  D2) 
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and  the  potential  (Do)  common  to  the  middle  three  strips  (see  Fig.  3.3).  Specifically, 
we  have  the  cost  functional 


min{Do^  +  Di^  +  D2^}. 

This  was  intended  to  reduce,  with  minimal  computational  overhead,  the  occurrence  of 
saturation  of  the  amplifiers  (which  were  limited  to  ±150  volts  per  strip),  thereby 
extending  the  region  of  control. 

The  vector  force  to  be  applied  by  the  driver  strips  is  then  (see  the  force  polygon  of  Fig. 
3.4) 


Fd  =  Fdes  -  {Fg  +  F[  +  Fe}. 

This  is  a  model-based  controller  formulation,  albeit  quasi-static,  since  the  dynamics 
(inertial  terms)  of  the  beam  do  not  enter  into  the  corrective  formulation.  The  potentials 
(the  three  Dj)  generating  Fd  are  then  determined  as  linear  combinations  of  the  Fd 
components,  subject  to  the  above  cost  functional.  In  the  case  of  amplifier  saturation, 
Fd  might  not  be  achievable.  Also,  recall  that  for  the  tilted  beam  case,  Fg  =  0  and  Fj  is 
referenced  to  the  equilibrium  position.  Moreover,  Fe  and  Fi  are  basecf  on  the  actual 
(fed  back)  position  of  the  beam.  No  feedforward  control  algorithms  [An,  et  al.,  1987], 
with  inertial  terms,  were  tried,  although  they  should  be  expected  to  enhance  SCOFSS 
performance. 

The  output  of  the  driver  box  (the  three  potentials)  are  then  transmitted  through  D/A 
converters  in  the  MINC  to  the  strip  amplifiers  . 


Fig.  3.4.  Schematic  of  driver  model.  Fd  is  the  force 
to  be  produced  by  the  driver  strips  so  as  to  yield  the  net 
(desired)  force  Fdes  (see  Fig.  3.3).  The  corresponding 
force  polygon  is  also  shown.  [1986] 

3.4.4.  Mechanical  System 

The  mechanical  system  consists  of  the  cantilevered  quartz  fiber  optic  beam  which  is 
effected  by  the  model-based  controller  described  above.  The  sum  of  all  forces  acting  on 
the  filament  gives  the  dynamics 


Fact  =  Fd  +  Fg  +  Fj  +  Fe  ~  Fdes 


where  Fact  is  the  actual  force,  which  should  approximate  Fdes  (amplifier  saturation 
may  affect  the  accuracy  of  the  approximation). 
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In  the  case  of  simulations ,  Fact  (neglecting  intrinsic  and  air  damping;  Newell,  1968) 
yields  the  "actual”  acceleration  (Aact)  via 

Aact  =  Fact/meff 

which  is  then  integrated  to  get  the  new  velocity  and  position.  The  effective  dynamic 
mass  m^ff  of  the  system  is  not  the  static  mass  of  Sec.  3.3.1,  but  rather  a  mass  which  if 
uniformly  distributed  over  the  electreted  region  would  yield  the  same  natural  frequency 
as  the  first  mode  of  the  actual  beam  (see  Sec.  7.1). 

The  actual  Filament  (with  controller  off)  is  highly  underdamped  (with  a  damping  ratio  of 
£  =  0.0125,  and  thus  a  quality  factor  of  Q  =  1/(2^)  =  40),  with  a  natural  frequency  of 
roughly  fn  =  36  Hz.  With  the  controller  on,  the  effective  damping  parameters  of  the 
beam  change  significantly  (see  Sec.  3.6.3). 


3.5.  Servo  Operation 

The  cantilevered  system  can  be  run  in  trajectory  or  regulator  modes,  as  described  below. 

3.5.1.  Trajectory  Operation 

The  experiments  presented  below  were  run  in  trajectory  mode.  That  is,  the  desired  set 
point  is  changed  continuously  in  time  (the  trajectory),  while  the  beam  tries  to  follow. 
The  set  points  representing  a  desired  pattern  were  transmitted  to  the  controller  at 
differing  speeds  (trajectory  scaling)  while  the  actual  track  of  the  beam  was  monitored. 
Some  of  the  results,  with  the  set  points  kept  within  the  passively  stable  region,  are 
discussed  below. 

Set  points  were  also  input  via  a  manually  operated  2DF  joystick.  The  display  of  the 
actual  beam  trajectory  was  monitored  on  an  oscilloscope. 

3.5.2.  Regulator  Operation 

An  alternative  mode  of  operation  is  to  specify  a  fixed  set-point.  External  disturbances 
are  then  measured  by  noting  the  magnitude  of  the  plate  voltages  required  to  keep  the 
fiber  tip  at  (or  near)  the  chosen  set-point.  Such  a  mode  was  not  investigated 
extensively,  although  it  might  be  used  as  a  simple  2D  accelerometer. 

3.6.  Experiments 

The  SCOFSS  system  has  been  operated  in  several  modes,  with  a  variety  of  filament 
charges,  over  a  range  of  speeds.  These  are  detailed  below.  The  complete  control  loop  cycle 
rate  through  the  MINC  is  on  the  order  of  600  Hz. 

3.6.1.  Linear  Charge  Density 

The  effective  linear  charge  density  within  the  polycarbonate  polymer  coating  on  the 
filament  was  necessarily  calculated  for  each  experiment  involving  control  (since  Fj  and 

Fd  are  X-dependent).  This  was  done  by  first  applying  ±37  volts  to  the  central  three 
driver  plates  (in  open-loop  mode)  and  monitoring  the  up  and  down  deflections  of  the 
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beam.  Then,  by  using  the  open-loop  model  of  forces,  a  value  for  X  was  determined 
which  caused  a  theoretical  deflection  of  the  beam  which  best  matched  the  observed 

deflection(s).  This  numerical  value  of  X  was  then  assigned  to  the  controller  model. 

3.6.2  Control  Space 

The  maximum  2D  region  in  which  the  beam  tip  can  be  moved  is  limited  by  the 
amplifiers  (which  saturate)  and  the  optimization  algorithm  for  the  driver  plates.  The 

region  also  depends  on  the  filament  charge  density  (X).  A  typical  region  is  shown  in 
Fig.  3.5. 


Fig.  3.5.  Region  of  maximum  controlled  beam  tip 
excursion  (5-plate  control).  Xo  denotes  the  equilibrium 
reference  position  of  the  beam  tip  with  all  plates 
grounded.  The  linear  charge  density  of  the  filament  was 

X  =  -0.069  nC/cm.  The  maximum  plate  voltages  were 
±150  volts,  with  the  optimization  algorithm  in  effect. 
The  lower  bound  represents  the  transition  to  the  unstable 
region  (feedback  on).  At  a  height  (measured  from  the 
plates  to  the  center  of  the  beam)  of  about  440  microns  the 
beam  transitions  to  latchdown.  The  corresponding  plate 
voltage  is  about  +70  volts  applied  to  the  three  central 
strips.  The  controller  could  have  applied  increasingly 
positive  potentials  to  increase  the  attraction  to  the  strips, 
but  this  just  puts  the  filament  deeper  into  the  unstable 
zone.  [1987] 


The  tapered  shape  of  the  region  is  probably  attributable  to  the  drop  in  producible  lateral 
field  strength  at  increasing  heights.  The  effects  of  the  controller  optimization  algorithm 
on  region  shape  were  not  analyzed. 

Increasing  X  increases  the  lateral  extent  of  deflection  and  the  maximum  height  which 
can  be  achieved  (it  produces  more  repulsive  force  when  the  plates  are  maximally 
negative)  while  at  the  same  time  increasing  the  minimum  height  which  can  be  stably 

achieved  (since  the  destablizing  image  force  increases  by  X^,  while  the  control  force  is 
proportional  to  only  X).  The  opposite  applies  for  a  decrease  in  X.  If  X  is  sufficiently 
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low  (of  order  0.05  nC/cm  or  less)  the  beam  cannot  be  forced  by  the  controller  into  a 
passively  unstable  region  (because  of  amplifier  limitations).  In  fact,  there  may  be 
passively  stable  regions  which  are  inaccessible. 

3.6.3  Step-input 

A  typical  step-input  response  of  the  system,  with  the  controller  on,  is  shown  in  Fig. 
3.6.  The  overshoot  of  the  step  is  about  30%,  which,  if  the  system  is  considered  to  be 

second-order-linear,  yields  an  effective  damping  ratio  of  £  =  0.36. 


Fig.  3.6.  Step-input  response.  The  controller  position 
gain  was  Kp  =  3.8  x  10*5  dynes/micron,  while  the 
velocity  gain  was  Kv  =  4.0  x  10'^  dynes-sec/micron. 
The  initial  height  of  the  filament  was  520  microns  with  a 
20  micron  vertical  step  input  applied.  Data  taken  from  an 
oscilloscope  photograph.  [1987] 


It  is  noted  that  there  is  no  noticeable  steady-state  error  for  step  inputs.  This  is  because 
the  reasonably  accurate  models  of  gravity,  image,  elastic  and  voltage-dependent  forces 
have  obviated  any  need  for  residual  position  error  to >.npensate  load  variations. 

The  smallest  steps  which  produce  a  noticeable  displacemt.it  on  an  oscilloscope  were  on 
the  order  of  one  micron  (this  is  also  the  resolution  limit  of  the  LEPD). 

3.6.4.  Sinusoidal  input 

In  this  test,  sine-waves,  discretized  into  twelve  steps  per  period,  were  transmitted  at  a 
range  of  frequencies,  from  0  to  25  Hz.  (this  is  well  below  the  observed  natural 
frequency  of  the  beam;  fn  ~  36  Hz.).  This  data  basically  demonstrates  the  tracking 
ability  of  the  model-based  controller  stored  in  the  MINC.  Typical  error  traces  are 
shown  in  Fig.  3.7. 
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Fig.  3.7.  Sinusoidal  errors  versus  frequency.  The 
input  amplitude  is  ±10  microns.  The  %  error  is  the 
difference  in  actual  and  desired  amplitudes,  divided  by 

the  desired  amplitude.  The  linear  charge  density  was  X  = 
-0.0465  nC/cm.  Data  taken  from  an  oscilloscope 
photograph.  [1987] 


3.7.  Higher  Modes 

In  order  to  assess  the  importance  of  higher  flexural  modes  on  the  dynamics  and 
controllability  of  the  beam,  computer  simulations  of  a  multi-mode  beam  were  investigated. 

3.7.1.  Mode  amplitudes 

Beam  vibration  is  a  problem  handled  by  classical  texts  [Timoshenko,  et  al.,  1974],  In 
short,  a  mode  of  free  vibration  is  presented  for  every  root  (e)  of  the  equation 


cos(e)  cosh(E)  +  1=0. 


These  roots  in  turn  yield  the  natural  modal  frequencies 


The  following  investigation  addresses  only  the  first  two  modes,  while  ignoring  higher 
modes.  Using  the  "assumed  modes  method"  presented  by  Craig  [1981],  the  mass, 
stiffness  and  forcing  function  matrices  can  be  calculated,  and  then  inserted  into  the 
mass-spring-driver  model  of  the  beam  dynamics,  which  yields: 


[mi  0  1  [d2vi/dt2l 

El 

'(El)4  0 

L  0  m2j  Ld2v2/dt2J 

+  U 

.  o  (E2)4. 

where  vl  =  vl(t)  and  v2  =  v2(t)  are  the  time  varying  mode  displacement  functions  to  be 
determined.  Note  that  no  damping  terms  are  included.  P  =  P(t)  is  the  magnitude  of  the 


MicroFields  Final  Report 


93 


distributed  load  (of  about  1  dyne/cm,  over  /)  as  determined  by  the  control  law  for  a 
sinusoidal  input,  El  =  1.358  (10^)  dynes-cm2,  and  m]  =  m2  =  m  =  pAL  where  A  = 
1.539  (10‘4)  cm2. 

The  effective  dynamic  mass  is  me^  which  is  that  mass,  uniformly  distributed  over 
length  /,  which  yields  the  same  natural  frequency  in  the  first  mode  as  that  of  the  actual 
beam  (with  the  total  mass  distributed  over  L).  Numerically,  using  a  computed  fn  ~  40 
Hz,  have 


meff  _  o .44  m 


where  m  =  pAL  is  the  static  mass. 

3.7.2.  Stability  considerations 

The  above  2-mode  model  was  investigated  using  a  computer  simulation.  The 
magnitude  of  the  tip  deflections,  for  modes  1  and  2,  over  a  range  of  frequencies,  are 
shown  in  Fig.  3.8.  It  is  noted  that  the  amplitude  of  mode  2  is  insignificant  compared  to 
mode  1,  and  thus  mode  2  can  be  discounted  as  a  significant  source  of  beam  instability. 
Thus,  a  simple  first-mode  model  of  beam  dynamics  suffices. 


Fig.  3.8.  Simulation  of  tip  deflections,  for  modes  1 
and  2,  versus  the  sinusoidal  input  frequency  applied  to 
the  driving  plates.  [1987] 

3.7.3.  Controller  delays 

The  role  of  delays,  such  as  exist  between  the  time  the  MINC  computer  receives  the  state 
estimate  (from  the  LEPD)  and  the  time  the  actual  potential  appears  at  the  plates  (about  a 
1.6  ms  delay),  was  investigated  using  the  above  2-mode  model.  Only  for  significant 
delays  (order  of  one  loop  period),  did  the  system  go  unstable.  And,  when  the  system 
did  go  unstable,  because  of  gains  and/or  delays,  the  second-mode  vibrations  remained 
small  (less  than  4%)  compared  to  the  first-mode  displacements. 
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3.8.  Discussion 

The  SCOFSS  system  has  served  as  a  useful  test  bed  for  addressing  issues  of  micro  electro¬ 
mechanical  system  fabrication,  analysis  and  control.  It  is  shown  that  sufficient  forces  can 
be  generated  between  electreted  regions  and  conductive  elements  driven  at  reasonable 
voltages  to  control  at  least  a  cantilevered  system.  Clearly,  further  analysis  can  be  given  to 
this  system.  However,  it  is  the  intention  of  the  CED  to  start  designing  and  fabricating  new 
MEMS  devices. 
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4.  OSCILLATORY  STABILIZATION  OF  A  SMALL  CANTILEVERED 
OPTICAL  FIBER  (OSSCOF) 

4.1.  Introduction 

Although  stable  electrostatic  configurations  are  possible  in  principle,  as  a  practical  rule 
electrostatic  equilibrium  will  be  unstable  [Wood,  et  al.,  1987;  Price,  1988;  Eamshaw, 
1842;  Denner  and  Pohl,  1982;  Jones  and  Kallio,  1979].  If  applications  of  electrostatic 
devices  (sensors,  actuators,  modulators,  etc.)  are  to  be  fully  explored,  some  means  of 
imposing  stabilization  on  a  passively  inherently  unstable  system  must  be  employed.  The 
use  of  active  feedback  is  widely  applicable,  versatile,  and  powerful.  But  it  is  complex;  it 
requires  sensing,  control  logic,  etc.  An  alternative  approach  to  stabilization,  "oscillatory 
stabilization",  does  not  share  these  weaknesses  of  feedback  stabilization.  We  discuss  here 
what  other  weaknesses,  and  what  strengths  are  inherent  in  oscillatory  stabilization. 

This  section  is  organized  as  follows:  Sec. 4. 2  describes  the  principles  of  oscillatory 
stabilization,  including  the  analytical  approach  (the  Mathieu  equation[McLachlan,  1947; 
Meirovitch,  1970])  and  a  useful  approximate  analysis.  In  Sec.  4.3  a  brief  description  is 
given  of  a  laboratory'  configuration  [Wood,  et  al.,  1987],  "SCOFSS",  used  to  test  many 
aspects  of  micro  electro-mechanical  systems.  The  mathematical  modelling  of  oscillatory 
stabilization  with  SCOF  (that  is,  SCOFSS  without  the  Servo  System  turned  on,  or 
OSSCOF  when  oscillatory  stabilization  is  applied)  is  described  in  Sec.  4.4,  along  with  the 
results  of  experimental  attempts  to  demonstrate  oscillatory  stiffening  and  stabilization.  One 
of  the  insights  gained  from  the  work  on  OSSCOF  was  an  improved  understanding  of  the 
stabilization  of  a  multimode  system.  This  is  described  in  Sec. 4. 5. 

4.2.  Oscillatory  Stabilization 

The  idea  of  oscillatory  stabilization  can  be  seen  in  the  following  simple  model.  Suppose  a 
particle  of  mass  m  moves  in  one  dimension  under  the  influence  of  two  springs,  one  of 

spring  constant  koc,  the  other  a  spring  which  has  a  spring  constant  kAccoscot  that  oscillates 

in  time  at  frequency  co.  The  position  of  the  particle  is  then  governed  by  the  equation 

md2x/dt2+(kDc+kAccoscot  )x  =  0.  (1) 

If  to  is  large  compared  to  the  natural  "DC  frequency"  of  the  system,  V kpc/m,  and  if 
kAc«  kDC,  then  intuition  would  suggest  that  the  only  effect  of  the  AC  spring  would  be  to 
add  high  frequency  jitter.  The  significant  motion  would  be  the  motion  at  the  DC  frequency, 
and  that  motion  would  be  stable  (unstable)  for  koc>0  (<0). 

In  fact  it  turns  out  that  intuition  fails.  The  oscillatory  spring  does  affect  the  stability.  This 
can  be  seen  with  a  simple  heuristic  analysis  [Wuerker,  et  al.,  1959].  We  view  the  motion 

as  the  sum  of  a  "slow"  part  xo(t),  with  characteristic  timescale  »l/co,  and  a  small  fast  part 

at  frequency  to,  so  that  we  can  write 


x(t)  =  xo(t)  +  Dcostot, 


(2) 
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where  D  is  the  (small)  amplitude  of  the  fast  oscillations.  We  then  substitute  this  into  (1) 

and  separately  set  to  zero  the  "fast"  (i.e.,  frequency  to)  part,  and  the  slow  [i.c.,  frequency 
characteristic  of  xo(t)]  part  of  the  resulting  equation.  We  find  that  the  slow  motions  are 
described  by  the  approximate  equation  for  xo(t) : 

md2xo/dt2  +  k^ff  xq=  0  ,  (3) 


where 


^  =  kDC+SS(^)  '  (4) 

The  oscillatory  spring  increases  the  effective  spring  constant  for  slow  motions  of  the 
system  and  hence,  contrary  to  intuition,  stiffens  the  system.  Furthermore,  the  condition  for 
the  motions  to  be  stable,  k*ff>0,  can  be  met  even  if  koc  is  negative. 

The  "fast  jitter"  approximation  above  can  make  predictions  about  stability  only  for  small 

^ac  (more  precisely,  kAc/m«to2).  For  more  general  results,  we  can  use  the  fact  that  (1) 
has  the  form  of  the  Mathieu  equation,  which  has  been  extensively  studied[McLachlin, 
1947;  Meirovitch,  1970],  Fig.  4.1  shows  regions  in  parameter  space  for  which  the 
general  solution  of  the  Mathieu  equation  is  stable  and  unstable.  For  small  kAC,  and  small 
negative  koc,  the  boundary  between  the  stable  and  unstable  regions  takes  the  form  of  a 
parabola  corresponding  to  k«ff  =  0  in  (4). 


Fig.  4.1.  Regions  of  stability  (shaded  and  instability 
for  solutions  of  the  Mathieu  equation. 

4.3.  SCOFSS:  In  Oscillatory  Mode  (OSSCOF) 

To  gain  familiarity  with  small  scale  electromechanical  devices,  the  Center  for  Engineering 
Design  constructed  a  miniature  (100  micron  scale)  device  in  which  an  optical  fiber  mounted 
as  a  cantilever  was  actuated  by  electrostatic  forces.  This  Small  Cantilevered  Optical  Fiber 
Servo  System,  or  "SCOFSS  "  is  shown  in  Fig.  3.1.  Details  are  presented  elsewhere 
[Wood,  et  al.,  1987].  The  driver  strips  shown  are  narrow  (200  micron  x  2  cm)  conducting 
plates,  insulated  from  each  other  by  thin  Si02  strips.  Electrical  fields  produced  by  the 
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driver  strips  interact  with  electrical  charge  embedded  in  the  fiber  (typically  1010  Coul 
spread  over  1.5  cm).  The  geometry  of  the  device  is  such  that  the  fiber  always  remains 
nearly  parallel  to  the  base,  and  remains  aligned  with  the  long  dimension  of  the  strips.  As  an 
adequate  approximation  we  treat  the  system  as  having  two  degrees  of  freedom,  shown  as  y 
(lateral  position)  and  z  (height)  in  Fig.  4.1.  In  oscillatory  operation  (OSSCOF),  no  LEPD 
feedback  is  employed. 

The  forces  acting  on  the  fiber  [1,2]  are  (i)  the  weight  of  the  fiber,  (ii)  the  elastic  force 
caused  by  the  bending  of  the  fiber,  (iii)  the  force  of  attraction  between  the  electrically 
charged  portion  of  the  fiber  and  its  image  in  the  conducting  base,  and  (iv)  control  forces, 
the  forces  between  the  charge  in  the  filament  and  the  (controllable)  potentials  applied  to  the 
driver  strips. 

As  a  test  bed  for  oscillatory'  stabilization  on  a  miniature  electrostatic  system,  SCOFSS  has 
two  very  desirable  features.  The  first  has  to  do  with  the  way  in  which  the  system  mixes 
electrostatic  and  (stabilizing)  elastic  forces.  When  the  fiber  is  positioned  high  (several 
hundred  microns)  over  the  base,  elastic  forces  tend  to  dominate  and  the  system  is 
’  passively"  stable  (i.e.,  stable  without  feedback  or  oscillations).  For  fiber  equilibrium 
positions  close  to  the  base,  electrostatic  forces,  especially  the  image  forces,  tend  to 
dominate  and  to  make  the  system  unstable.  By  choosing  the  location  of  the  equilibrium 
point  we  can,  therefore,  control  the  passive  stability  of  the  system. 

A  second  very  important  feature  of  SCOFSS  is  that  each  of  the  ten  driver  plates  can  be 
supplied  with  a  different  voltage  (in  the  range  ±150  volts).  As  we  shall  show  presently, 
this  flexibility  is  essential  for  creating  electric  fields  that  produce  oscillator)'  stabilization. 

To  use  SCOFSS  for  oscillator)’  stabilization,  a  set  of  DC  voltages  is  applied  to  the  strips  to 
establish  an  equilibrium  point,  at  which  elastic,  electric,  and  weight  forces  balance,  though 
not  necessarily  stably.  The  simplest  choice  is  to  center  the  negatively  charged  filament  over 
a  plate,  to  apply  to  that  plate  some  positive  voltage  Vo,  and  to  leave  all  other  plates 
grounded.  The  position  of  the  equilibrium  point,  of  course,  depends  on  Vo.  As  the  voltage 
is  increased  the  equilibrium  point  moves  downward.  Oscillations  can  now  be  added  by 
applying  AC  voltages  to  some  of  the  plates.  For  the  frequency  range  in  which  we  are 
interested,  the  resulting  fields  are  essentially  electrostatic,  though  time  varying,  i.e., 
induction  and  radiation  effects  are  negligible. 

Unless  the  voltage  configuration  is  chosen  rather  carefully  the  oscillations  cannot  be 
stabilizing  in  the  sense  of  Sec. II.  Corsider,  for  example,  what  happens  if  we  simply  make 
Vo  a  superposition  of  a  DC  and  an  AC  voltage.  As  Vo  oscillates  the  equilibrium  point  will 
oscillate,  and  the  fiber  will  necessarily  oscillate  at  the  imposed  AC  frequency  (possibly  in  a 
stable  limit  cycle).  The  motion  of  the  fiber  will  not  be  described  by  (1)  which  requires  that 
the  equilibrium  point  (xo  =  0)  be  fixed.  This  means  that  at  the  equilibrium  point  the  AC 
forces,  as  well  as  the  static  forces,  must  vanish. 

In  principle  we  can  arrange  this  by  making  Vo  static  and  adjusting  it  so  that  the  equilibrium 
point  is  at  a  height  V3/4  times  the  width  of  a  driver  plate.  At  this  special  position  the  fiber 
turns  out  to  be  unaffected  by  voltages  applied  symmetrically  to  the  two  driver  strips 
adjacent  to  the  central  strip.  If  the  voltage  Vi,  shown  in  Fig.  4.2a,  is  oscillatory,  we  will 
then  have  satisfied  the  requirement  that  the  equilibrium  point  be  fixed,  i.e.,  that  there  be  no 
oscillatory  force  at  the  equilibrium  position.  Small  deviations  from  equilibrium,  both  in  the 
y  and  the  z  directions,  will  be  described  by  (1).  The  constant  kAC  is  proportional  to  Vj  and 
depends  on  the  spatial  variation,  near  equilibrium,  of  the  control  force  due  to  the  symmetric 
pair  of  plates.  The  constant  k^c  will  be  determined  by  the  details  of  the  elastic  force,  the 


MicroFields  Final  Report 


98 


image  force,  and  the  control  force  due  to  Vo,  near  equilibrium.  (A  static  component  of  Vj 
would  also  influence  koc  ) 

In  practice  this  three-plate  configuration  is  too  inflexible.  It  is  preferable  to  have  an 
adjustable  equilibrium  point.  This  can  be  accomplished  with  the  five-plate  configuration 
shown  in  Fig.  4.2b.  In  this  arrangement  the  equilibrium  height  is  set  through  the  choice  of 
the  voltage  Vo-  The  ratio  of  to  V2  is  then  determined  analytically,  and  fine  tuned 
experimentally,  so  that  there  is  no  control  force  at  the  (adjustable)equilibrium  point  caused 
by  these  two  pairs  of  plates.  If  Vj  and  V2  are  made  oscillatory,  the  dynamics  of  small 
deviations  in  both  the  y  and  the  z  directions  are  described  by  (1).  Nonoscillatory 
components  could  also  be  added  to  V[  and  V2  to  alter  the  lateral  position  of  the  equilibrium 
point.  This  five-plate  configuration  was  used  for  experimental  testing  with  SCOFSS. 


Fig.  4.2.  Voltage  configurations  for  oscillatory 
stability,  a)  A  three-plate  configuration,  b)  A  five-plate 
configuration.  [1986] 


4.4.  Results 

4.4.1  Computer  Simulations 

The  first  question  investigated  for  SCOFSS  was  how  effective  the  oscillator)' 
stabilization  could  be  in  principle.  To  answer  this  question  the  forces  acting  on  the 
fiber  were  modelled  mathematically,  with  the  fiber  idealized  as  perfectly  rigid  and  free- 
floating,  i.e.,  not  cantilevered  but  constrained  mathematically  to  remain  orthogonal  to 
the  y-z  plane.  Elastic  forces  were  ignored  so  that  stability,  if  achievable,  could  be 
credited  to  oscillatory  stabilization. 

The  flavor  of  the  mathematical  problem  can  be  seen  in  the  form  of  the  equation  for 
lateral  motion: 


d^/dt2  =  V0fo(y,z)  +  Vjfi(y.z)  +  V2f2(y,z). 


(5, 
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Here  y  and  z  are  the  lateral  and  vertical  coordinates  of  the  fiber,  the  V,  are  the  voltages 
shown  in  Fig.  4.2,  and  'he  f;  are  algebraic  functions.  The  expression  for  the  vertical 
motion  is  similar  but  terms  also  appear  representing  weight  and  image  forces.  Let  us 

choose  V]  and  V2  to  have  time  dependence  cos  cut  and  to  have  the  correct  ratio  so  that 
there  is  no  oscillatory  force  at  the  equilibrium  point.  For  convenience  let  us  also  choose 
y  and  z  to  be  zero  at  the  equilibrium  point.  Equation  (5)  does  not  have  the  form  of  ( 1 ); 
it  takes  that  form  only  if  it  is  linearized  in  the  deviations,  i.e.,  if  it  is  expanded  for  small 
y  and  z,  and  if  only  terms  linear  in  these  variables  are  kept.  Equation!  1),  the  Mathieu 
equation,  describes  the  dynamics  only  very  close  u»  equilibrium.  To  get  more  general 
results  about  stability,  and  to  see  how  valid  the  predictions  for  stability  are  from  the 
(linearized  ap;  roximation)  Mathieu  equation  theory,  we  need  (5)  and  the  corresponding 
equation  for  vertical  motion. 

These  equations  w'ere  studied  on  a  computer  with  a  standard  differential  equation  solver 
(DVERK  from  the  scientific  subroutine  package  IMSL).  The  voltage  Vo  was  varied  to 
select  equilibrium  points  at  various  heights  in  the  median  plane  of  the  five-plate 
configuration.  The  point  representing  the  fiber  was  then  started  from  rest  several 
microns  from  the  equilibrium  point.  The  trajectory  of  the  fiber  point  was  monitored  as 
its  position  was  evolved  numerically ,  and  it  was  determined  whether  the  motion  was 
stable.  At  each  equilibrium  point,  several  combinations  were  tried  for  the  magniude 
and  frequency  of  the  oscillatory'  voltages.  Since  the  initial  deviations  from  equilibrium 
were  small  (several  microns)  the  predictions  of  linearized  theory'  were  expected  to  apply 
and  were  found  to  apply.  The  motion  w;as  stable  or  unstable  in  accordance  with 
predictions  based  on  the  Mathieu  equation  (Fig.  4.1). 

To  investigate  the  effectiveness  of  oscillatory  stabilization  another  set  of  simulations 
was  made.  The  oscillating  voltage  available  was  taken  to  be  1500  volts  and  the  charge 
on  the  fiber  to  be  10' 10  Coul  spread  over  1  cm.  The  maximum  initial  deviation  was 
determined  for  which  stability  could  be  achieved  for  these  parameters.  The  results 
indicated  that  the  effectiveness  of  oscillatory  stabilization  depends  on  height.  At  an 
equilibrium  height  of  200-300  micron  (1 -11-plate  widths)  the  motion  can  be  stabilized 
for  initial  deviations  of  ~30  microns;  the  range  of  stabilizable  initial  deviations  decreases 
with  both  greater  and  smaller  height. 

4.4.2  Experimental  Results 

To  stabilize  the  motion  of  a  charged  SCOF  the  following  experimental  protocol  was 
tried:  (i)  With  all  driver  strips  grounded,  the  fiber  was  set  up  to  be  parallel  to  the 
substrate,  under  the  combined  action  of  elastic,  weight,  and  image  forces,  at  an  initial 
calibration  height  of  520  microns,  (ii)  The  control  voltage  Vo  was  increased,  thereby 
lowering  the  equilibrium  point,  (iii)  A  transition  height  zo  (typically  around  150 
microns)  was  determined  beiow  which  the  equilibrium  is  unstable,  (iv;  Oscillatory' 
voltages  Vi  and  V2  are  applied  with  the  correct  V1/V2  ratio  so  that  there  is  no 
oscillatory  force  at  the  equilibrium  point;  the  stabilizing  oscillations  are  first  applied  at 
an  equilibrium  point  above  zo;  as  the  equilibrium  point  is  lowered  (by  increasing  Vo) 
the  V 1/V2  ratio  is  adjusted  to  insure  that  there  is  no  oscillatory  force  at  equilibrium. 

In  principle  this  approach  should  have  allowed  us  to  keep  the  fiber  stable  as  we  lowered 
it  below  zo-  In  practice,  this  could  not  be  done  in  a  predictable  and  repeatable  manner. 
The  experimental  obstacle  was  noise  due  to  mechanical  vibration  of  the  SCOFSS 
apparatus,  even  after  the  apparatus  had  been  mounted  on  vibration  absorbing  pads.  At 
an  equilibrium  \  oint  jusi  above  the  transition  height  zo,  the  fiber  was  observed  to 
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undergo  vibration-induced  oscillations  of  around  30  microns.  Mathieu  equation 
theory  (Fig.  4.1)  was  used  to  determine  how  far  below  zo  the  motion  could  be 
stabilized.  The  answer,  20-30  microns  shows  that  30  micron  vibrational  noise  makes 
consistent  stabilization  impossible. 


Fig.  4.3.  The  stiffening  of  a  fiber  by  an  AC  driving 
force.  Curve  a)  shows  the  natural  oscillation  of  the  fiber 
with  no  AC  force  applied.  Curve  b)  represents  the  same 
experimental  conditions  except  that  an  AC  force  has  been 
added;  the  AC  force  slightly  increases  the  natural 
oscillation  frequency.  [1987] 


A  demonstration  of  the  stiffening  effect  of  oscillations,  however,  is  possible.  Fig.  4.3 
shows  the  vertical  free  vibrations  of  the  fiber  (tip  height  as  a  function  of  time).  The 
equilibrium  position  for  both  curves,  roughly  175  micron  above  the  middle  of  the 
center  strip,  is  slightly  above  the  stable/unstable  transition  height  zo-  Curve  a  shows 
the  motion  in  the  absence  of  oscillatory  voltages.  The  frequency  (-35  Hz)  is 
significantly  lower  than  that  (~77  Hz)  of  the  Fiber's  fundamental  elastic  mode,  due  to 
the  destabilizing  (i.e.,  frequency  lowering)  effects  of  the  electrostatic  interactions. 
Curve  b  represents  the  same  conditions  but  with  oscillatory  voltages  V!  and  V2,  applied 
in  the  correct  ratio  and  with  the  maximum  strength  allowed  by  the  150  volt  limitation. 
The  decrease  (-10%)  in  frequency  of  curve  b  compared  to  curve  a  shows  that  the 
oscillations  stiffen  the  system,  but  not  very  dramatically. 

4.5.  Multimode  Systems 

The  SCOF  has  an  infinite  number  of  degrees  of  freedom  associated  with  the  infinite 
number  of  its  dynamical  modes.  The  treatment  of  oscillatory  stabilization  given  so  far  treats 
the  fiber  as  a  system  with  only  two  (decoupled)  degrees  of  freedom  corresponding  to 
vertical  and  lateral  motions.  This  treatment  is  justified  only  if  all  forces  are  distributed 
along  the  fiber  in  such  a  way  that  only  a  single  mode  is  ever  excited,  or  if  the  oscillation 
frequencies  are  very  small  compared  to  the  second  elastic  mode  of  the  fiber.  Neither  of 
these  conditions  was  satisfied  in  the  SCOFSS  or  OSSCOF  experiments.  To  clarify  the 
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meaning  of  those  experiments  we  present  here  a  treatment  of  the  effect  of  oscillations  on  a 
linear  system  with  many  coupled  degrees  of  freedom,  a  subject  which  seems  not  to  have 
previously  appeared  in  the  literature  on  driven  oscillations  of  multimode  systems 
[Takahashi,  1979;  Takahashi,  1981;  Sinha  and  Chou,  1979;  Bolotin,  1964]. 

We  assume  that  an  N  mode  system  (where  N  may  be  infinite)  can  be  described  by  the  basis 
vectors  Uj  (i=  1  to  N).  For  an  elastic  beam  the  vectors  Uj  =U;(x)  would  be  a  set  of  mode 
shapes  for  the  beam.  The  state  of  the  system,  as  a  function  of  time,  is  given  by 

N 

Iu,(t)Ui .  (6) 

i=0 

Let  us  suppose  that  we  have  chosen  the  Uj  to  be  modes  of  the  system  which  would  be 
normal  in  the  absence  of  oscillatory  forces,  then  only  the  oscillatory  force  can  couple  the 
modes.  To  write  the  multimode  analog  of  (1)  we  define  the  mass  associated  with  every 

mode  to  be  unity  and  we  denote  by  kf  the  (time  independent)  spring  constant  associated 

with  mode  i.  We  denote  by  k *  the  coupling  constant  that  describes  the  oscillatory 

influence  on  mode  i  of  an  excitation  in  mode  j.  The  general  equation  of  motion  then  has  the 
form 


N 

d2Ui/dt2+k|Xu1  +cosG)t^Tk^CUj=0.  (7) 

j=l 

If  the  Uj  were  normal  modes  for  the  oscillatory,  as  well  as  for  the  time  independent 

Af'' 

interactions,  then  k*  would  be  a  diagonal  matrix,  but  this  is  not  generally  the  case.  When 
a  "fast  jitter"  approximation 


Uj  =  u^ +Djcosox 


(8) 


like  that  of  (2)  is  used,  we  find  that  the  "slow"  motions  are  described  by 

N 


AC,  AC  0 


d2u?/dt2  +  kfcu?  +  -i-VkXu;  =0-  (9) 

’  20)2^  ,J  J 

j,k=l 


The  analog  here  of  kefy  in  (3)  and  (4)  is  the  set  of  normal  mode  frequencies  for  (9).  These 

are  the  values  of  Q  for  which  the  determinant  of  the  matrix 

N 

Aik=(kfc-D2)6ik+J-^kACkAC  (10) 

j=l 


AC 

vanishes.  If  we  are  to  be  consistent  with  the  approximations  already  used  we  must  take  k;j 
to  be  small  and,  in  computing  the  eigenvalues,  we  must  keep  terms  only  up  to  second  order 
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in  this  small  quantity.  With  this  simplification  the  i^  eigenvalue  has  the  relatively  simple 
form 


N 


.AC.  AC 

Kij  Kji  • 


>1 


(ID 


In  the  absence  of  driven  oscillations  (i.e.,  with  k^  =0)  the  frequencies  are  the  obvious 

normal  mode  frequencies  for  the  system  described  by  (7).  The  summation  term  in  (11) 
gives  the  changes  to  those  frequencies  induced  by  oscillations.  The  i111  mode  will  be  stable 

2 

or  unstable  according  to  whether  Q*  is  positive  or  negative.  As  a  check  on  the  validity  of 

the  approximation  scheme,  the  behavior  of  a  two  mode  linear  system  [N=2  in  (7)]  was 
numerically  simulated  to  determine  under  what  conditions  oscillations  would  add  stability. 

AC 

The  results  were  in  excellent  agreement  with  (11)  except  when  k;j  was  significantly  larger 

than  k  ,  in  which  case  the  approximation  should  be  expected  to  fail.  We  now  turn  from 

the  general  case  of  oscillatorily  coupled  normal  modes,  to  the  specific  case  of  an  elastic 
beam.  The  typical  set  of  equations  for  an  N  mode  linear  elastic  system  is 

£  N 

XMijd2Uj/dt2  +  £KyUj  =  0,  (12) 

j=l  j=I 

where  My  is  the  mass  matrix  and  Ky  is  the  stiffness  matrix  (Craig,  1981].  Equation  (7) 
differs  from  this  general  form  only  in  that  My  =  8y  (by  the  choice  of  basis  modes)  and  in 
that  there  are  two  stiffness  matrices.  One  matrix,  Ky  =  ki  8y,  represents  the  time 

AC 

independent  stiffnesses.  The  other  stiffness  matrix  Ky=coscotkjj  has  an  unfamiliar 

oscillatory  multiplier,  but  aside  from  this  it  is  computed  in  the  same  way  as  any  other 
stiffness  matrix.  Like  any  other  stiffness  matrix  it  must  therefore  be  symmetric  and  hence 

kjj  ~kj|  .  With  this  symmetry  condition  the  summation  in  (1 1)  can  be  seen  to  be  a  sum  of 
2  DC 

squares  so  that  Qj  is  always  greater  than  kj  .  The  effect  of  oscillations  is  always 
stabilizing. 


DC 

In  the  OSSCOF  experiments,  the  lowest  kj  of  the  fiber,  that  corresponding  to  i=l  (the 

fundamental  mode),  became  negative  at  heights  below  zo-  The  other  kj  ,  corresponding  to 

stiffer  modes,  remained  positive.  If  the  i=l  summation  in  (1 1)  were  large  enough  to  make 
2 

£>i  positive,  then  the  fiber  would  be  stable.  The  oscillations  cannot  destabilize  the  other 
modes.  This,  of  course  is  true  of  any  elastic  system  (more  generally,  any  system  with  a 

A /"* 

symmetric  kjj  ).  That  is,  if  only  one  mode  of  such  a  system  is  unstable  we  need  not 
concern  ourselves  with  the  other  modes. 
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This  result  means  that  we  can  view  the  problem  of  stabilizing  the  SCOF  as  qualitatively  the 
same  as  stabilizing  a  single  mode  system.  There  are,  however,  quantitative  differences.  In 
particular,  the  summation  in  (1 1)  contains  stabilizing  contributions  from  all  modes,  not  just 
the  fundamental  mode.  In  practice  this  is  of  minor  importance  in  the  OSSCOF  experiments 
for  two  reasons.  First,  simulations  showed  that  even  at  fairly  high  driven  frequencies, 
fiber  motion  in  the  higher  modes  was  much  smaller  (except  close  to  a  higher  mode 
resonance)  than  motion  in  the  fundamental  mode.  Second,  the  prevalence  of  vibrational 
noise  in  the  experimental  apparatus  made  it  difficult  in  any  case  to  test  quantitative 
agreement  of  theory  and  experiment. 

4.6.  Conclusions 

The  application  of  oscillatory  stabilization  to  a  small  electro-mechanical  system  has  been 
studied  analytically  and,  with  the  SCOFSS/OSSCOF  apparatus,  experimentally. 
Oscillatory  stabilization  requires  that  an  AC  signal  can  be  applied  to  a  system  without 
shifting  the  system's  equilibrium  point  and  this,  in  tum,  requires  that  there  be  considerable 
flexibility  in  the  way  in  which  electrical  fields  can  be  shaped.  In  the  OSSCOF  apparatus 
this  flexibility  resided  in  the  large  number  of  separately  adjustable  conductor  plates. 

In  the  SCOFSS/OSSCOF  apparatus  the  length  scale  for  variation  of  the  electrical  fields  was 
on  the  order  of  100  microns.  This  length  scale  sets  an  approximate  upper  limit  on  the 
deviations  from  equilibrium  which  can  be  stabilized  with  oscillations.  Both  experimental 
and  laboratory  investigations  showed  that  the  maximum  deviations  (20-30  microns)  that 
could  actually  be  stabilized  were  rather  smaller  than  the  upper  limit.  The  situation  for 
oscillatory  stabilization  is  to  be  contrasted  with  stabilization  through  model-based  feedback 
control,  which  has  no  a  priori  length  scale  dictated  limitations  on  the  maximum  deviations 
that  can  be  stabilized.  In  applications,  the  relative  merits  of  oscillatory  stabilization  and 
feedback  control  may  depend  on  how  homogeneous  the  fields  are  and  on  how  noisy  the 
system  is. 

One  advantage  of  oscillatory  stabilization  in  dealing  with  miniature  mechanical  systems 
involves  the  multiple  modes  of  an  elastic  element.  An  analytic  investigation  shows  rather 
generally  that  the  motion  of  all  modes  of  the  element  are  oscillatorily  stabilized  if  the  least 
stable  mode  of  the  element  is  stabilized. 
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5 .  ACTUATORS 
5.1  Background 

The  initial  motivation  for  the  Microfield  Project  was  the  development  of  actuators. 
Over  the  past  few  years,  we  have  generated  concepts  for  a  number  of  alternatives  and 
studied  their  potential  performance.  Although  the  force  calculations  are  difficult  (as 
shown  in  the  following  sections)  we  are  nevertheless,  very  optimistic  about  the 
possibility  of  developing  actuators  with  very  interesting  properties.  Scaling  arguments 
imply  that  high  performance  systems  can  be  developed. 

The  actuator  project  is  very  interesting,  but  it  must  be  emphasized  here  that  efforts  must 
proceed  carefully  and  at  a  realistic  pace.  Actuators  are  the  most  complicated  systems 
proposed  and  must  integrate  the  function  of  all  elements  considered  including  sensors, 
controllers,  force-generating  elements,  etc.  Therefore  our  plan  is  to  continue 
developing  our  analytical  capabilities  so  that  systems  can  be  designed  which  will  have  a 
good  probability  of  successful  function.  We  also  plan  to  continue  developing  our 
capability  to  fabricate  micro  systems,  and  as  confidence  is  gained,  we  will  generate 
physical  components  for  evaluations  purposes. 

We  have  considered  several  approaches  for  the  investigation  of  force  production.  In 
the  cases  below,  force  is  produced  by  the  interaction  of  electrically  active  structures 
composed  of  conductors,  dielectric  materials,  electrets,  ferroelectrics,  and  other 
materials. 


5.2  Forces  on  Conductors,  Dielectrics  and  Ferroelectric  Materials 

In  many  of  the  devices  previously  analyzed,  or  proposed,  it  is  necessary  to  calculate  a 
force,  such  as  for  the  deflection  of  a  cantilevered  beam,  or  for  an  actuator.  However, 
even  though  the  potential  functions  and  fields  within  a  device  may  be  known,  the 
calculation  of  forces  is  yet  a  higher  level  of  integration,  and  hence  analytical  difficulty. 
Moreover,  the  forces  are  often  attenuated  by  dielectric  or  charge-compensating  effects. 

Below  are  force  calculations  for  a  variety  of  geometries  and  materials  pertinent  to  micro 
systems  (actuators  in  particular).  In  some  cases,  such  as  for  ferroelectrics,  the  forces 
are  theoretically  quite  high.  However,  experimentation  has  shown  these  forces  to  be 
substantially  less,  although  not  useless.  Clearly,  the  analysis  and  experimentation 
should  converge  in  order  to  make  more  rational  predictions  about  the  proposed  devices. 


5.2.1  Force  Between  Parallel  Conductive  Plates 

In  the  most  elementary  form,  the  longitudinal  attractive  force  between  two  off-set 
conductors  (of  width  w  and  spacing  s;  see  Fig.  5.1)  is  given  by  [Zahn,  1979,  p. 
253] 


Fx  =  Eo  (w/s)  V02/2 

which  is  seen  to  be  independent  of  the  amount  of  overlap  x  (provided  there  is  some) 
and  length  of  the  conductive  plates. 
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Fig.  5.1.  Overlapping  conductor  plates.  [1986] 


If  potential  difference  V0  is  determined  by  the  breakdown  field  of  air,  Eb  =  Vo/s, 
then 


Fx  =  £o  (ws)Eb2/2  . 

But  Ac  =  (ws),  which  is  the  "cross-sectional"  area  of  the  plates.  Thus,  the  stress 
produced  by  this  configuration  is 

=  Fx/Ac  =  £oEb2/2 

Numerically,  Eb  =  3  x  106  V/m,  which  yields 
cx  =  6  x  10’ 3  psi  (40  Pa) 

which  is  quite  weak  compared  to  human  muscle  which  can  generate  about  102  psi. 
Note  that  the  breakdown  field  can  be  increased  to  about  50X106  V/m  (or  even 
l.OOOxK^V/m)  by  putting  either  compressed  gas  (such  as  hydrogen)  or  a  vacuum 
between  the  plates. 

However,  the  force  of  the  conductor  system,  since  being  an  "end-effect,"  can  be 
increased  as  shown  in  the  configuration  below  (Fig.  5.2). 


Fig.  5.2.  Repeated  overlapping  metal  plates  added 
mechanically  in  parallel.  [1986] 


A  "grain-size"  A  can  be  defined,  which  then  gives  a  force  per  volume  of 


Fx/Vol  =  EoEb2/(4A) 
which  is  plotted  in  Fig.  5.14. 
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Obviously  making  the  elements  smaller  increases  the  actuator  strength.  This 
geometry  can  be  extended  to  an  array  of  planes  (MIPA)  as  shown  in  Fig.  5.3,  and 
discussed  in  section  6. 


Fig.  5.3.  Multiple-interleaved-plane  actuator  (MIPA).  [1982] 


The  force  per  volume  for  this  structure  should  be  essentially  the  same  as  for  the 
two-plane  array  above  since  adding  more  plates  in  parallel  does  not  increase  the 
effective  cross-sectional  area.  However,  putting  the  planes  closer  together  does 
have  the  benefit  of  decreasing  the  required  voltages. 

It  should  be  noted  that  in  addition  to  the  longitudinal  force,  there  are  voltage- 
dependent  transverse  attractive  forces  between  elements,  in  addition  to  image  forces 
between  active  elements  and  the  supporting  substrates.  The  image  forces  are 
always  attractive  and  will  thus  tend  to  compress  the  actuator  planes.  This  effect  can 
be  minimized  by  making  the  substrates  of  a  low-dielectric  material,  such  as 

Teflon®.  In  fact,  the  transverse  image  forces  are  roughly  the  same  as  the 
longitudinal  force.  Thus,  with  a  coefficient  of  friction  less  than  unity,  the 
transverse  forces  should  dominate,  and  thus  provide  contractile  actuation. 

5.2.2  Forces  on  Dielectric  Materials 

One  can  postulate  an  actuation  system  consisting  of  parallel  conductive  planes,  at  a 
potential  difference  Vo,  which  have  a  dielectric  material  (with  dielectric  constant  e) 
sandwiched  between  them.  Since  there  are  two  types  of  dielectric  materials,  there 
are  two  types  of  force  analyses.  The  first  analysis  is  for  a  "linear"  dielectric 
material  where  in  the  dipoles  of  the  material  are  ordered  depending  on  the  field 
strength  and  direction.  The  second  analysis  is  for  a  permanently  polarized  material 
which  has  a  polarization  which  changes  only  when  the  material  is  "coerced"  by  a 
sufficiently  strong  field. 


For  the  linearly  polarized  material,  the  geometry  would  look  like  Fig.  5.4. 
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Fig.  5.4.  Linear  dielectric  material  sandwiched  between 
conductor  plates.  [1986] 


The  force  in  the  longitudinal  direction  is  [Zahn,  p.219] 

Fx  =  (£  -  £o)  V02  (w/2s)  =  (e  -  Eo)  Eb2  (ws)/2 

where  Vb  =  EbS  is  the  breakdown  potential  in  air.  Dividing  by  the  cross-sectional 
area,  Ac  =  (ws),  yields  the  "stress" 

Ox  =  (e  -  Eo)  Eb2/2 

For  a  material  with  a  relative  dielectric  constant  of  say  9,  such  as  for  sapphire,  we 
have  ax  =  4.6  x  10'2  psi.  This  is  stronger  than  for  plates  without  a  dielectric,  but 
still  weak. 

If  the  above  "single-cell"  geometry  is  repeated  as  per  Fig.  5.2,  the  force  per  volume 
scales  with  grain  size  as 

Fx/Vol  =  (E  -  Eo)  Eb2/(2A) 

This  is  plotted  in  Fig.  5.14,  and  is  seen  to  be  stronger  than  for  the  pure  conductor 
plates.  This  is  essentially  due  to  the  increased  permeability  of  the  dielectric  over  air. 
Similar  analyses  hold  for  cylindrical  linear  dielectrics  between  two  cylindrical 
conductors,  such  as  shown  in  Fig.  5.5,  which  might  approximate  fiber-in-hole  type 
actuators. 
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Fig.  5.5.  Linear  Dielectric  between  Conductive  Tubes.  [1986] 


5.2.3.  Forces  on  Ferroelectric  Materials 

For  a  permanently  polarized  material  placed  between  parallel  conductor  plates,  as  in 
Fig.  5.6,  the  analysis  is  different,  yielding  two  terms  [Zahn,  p.  218], 

Fx  =  P0V0w  +  P20sw/£o  . 

Note  that  the  First  term  is  "active”  (i.e.  a  function  of  the  applied  voltage)  while  the 
second  term  is  "passive”  (independent  of  voltage). 


Fig.  5.6.  Permanently  polarized  material  sandwiched 
between  conductor  plates.  [1986] 


The  corresponding  stresses  are,  with  V0  =  Ebs, 

cx  =  PoEb  +  P02/£o  =  Ox(act)  +  Cx(pass) . 

A  strongly  polarized  material  such  as  lithium  niobate,  with  P0  =  3  C/m2  =  3  x  105 
nC/cm2,  yields 


Ox(act)  =  1.3  x  103  psi 
cx( pass)  =  1.5  x  108  psi . 
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Obviously,  the  "passive"  force  swamps  the  "active"  force.  The  same  holds  for  a 
weaker  ferroelectric  material,  such  as  PZLT,  with  P0  =  0.1  C/m2,  which  yields 

ox(act)  =  4.3  x  101  psi 

Cx(pass)  =  1.7  x  105  psi . 

For  still  weaker  ferroelectrics,  the  active  force  could  overcome  the  passive  force, 
thus  allowing  "pushing"  as  well  as  "pulling"  forces. 

The  force  per  volume  is, 

Fx/Vol  =  (PoEb  +  P02/£o)/2A 
with  the  results  for  PZLT  plotted  in  Fig.  5.14. 

Similar  analyses  hold  for  cylindrical  ferroelectrics  between  two  cylindrical 
conductors,  such  as  shown  in  Fig.  5.7. 


Fig.  5.7.  Permanently  polarized  material  between 
conducting  plates.  [1986] 


Obviously,  by  inspection  of  Fig.  5.14,  the  force  densities  for  ferroelectrics  could 
rival  those  of  muscle  for  sufficiently  small  elements.  However,  due  to  the 
limitations  of  breakdown  and  charge  compensation,  it  is  not  likely  that  such  high 
forces  would  be  obtained  without  going  to  still  smaller  structures. 

5.3.  Forces  Generated  by  MultLelement  Planar  Actuators 


5.3.1  Introduction 

A  number  of  actuator  systems  we  have  considered  utilize  different  geometries  to 
produce  the  same  result.  That  is  the  deployment  of  multiple  and  controllable 
charged  areas  in  armature  and  stator  structures  such  that  a  sustained  force  can  be 
generated.  Fig.  5.8  schematically  shows  that  situation.  Unfortunately,  this 
geometry  is  exceedingly  difficult  to  analyze,  yet  it  is  important  to  generate  at  least 
some  analytical  technique  which  can  predict  the  force  levels  which  can  be  expected 
from  these  machines. 
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In  order  to  make  force  estimates,  an  approximate  method  which  transforms  the 
system's  structure  into  a  form  more  suitable  for  analysis  will  be  described  in  the 
following  figures.  Fig.  5.9  shows  a  system  with  fewer  active  elements  which  are 
connected  by  symbolic,  non-interacting,  links.  Fig.  5.10  shows  those  conceptual 
links  switched  such  that  the  force  generated  should  be  identical  to  the  previous  case 
except  that  the  system  is  incapable  of  moving  more  than  one  element  width.  Figure 
5.1 1  shows  the  charged  regions  now  blended  into  planar  sections.  An  analysis  has 
shown  that  in  this  circumstance  such  an  assumption  is  reasonable.  That  is,  the 
force  predictions  made  will  be  accurate  to  better  than  one  order  of  magnitude.  It  is 
now  possible  to  use  this  geometry  for  two  designs.  The  first  design  utilizes 
interdigitated  conductor  plates  whose  potential  is  varied  such  that  the  system 
produces  force.  The  second  design  uses  parallel  conductor  plates  interspersed 
within  charged  polarized  planes. 

Fig.  5.12  illustrates  a  series  of  conductive  plates  each  charged  to  specific  potentials. 
This  system  can  be  analyzed,  including  the  effects  of  dielectrics  and  image  forces, 
to  gain  an  estimate  of  force  which  might  be  produced  by  a  system  constructed  of 
conductors  only.  Fig.  5.13  illustrates  a  similar  system  except  that  the  stator 
elements  are  now  replaced  by  radically  charged  materials  with  the  armatures 
composed  of  conductors  which  can  be  adjusted  to  appropriate  potentials.  These 
cases  are  discussed  below. 


r-i-»i -t -^1 -Ni-»  i  ■*  i  — 1 
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Fig.  5.8.  MIPA  consisting  of  charged  segments  of 
alternating  sign.  [1982] 


Fig.  5.9.  Idealized  MIPA.  [1986] 
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Fig.  5.10.  Transformed  MIPA  such  that  each  plane  now 
has  the  same  sign  throughout.  [1986] 


Fig.  5.11.  Transformed  MIPA  with  the  discretized  planes 
now  approximated  by  continuous  conductor  sheets.  [1986] 
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Fig.  5.12.  Transformed  MIPA  representation.  [1986] 
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Fig.  5.13.  Transformed  MIPA  with  the  stator  plates  taken 
to  be  fixed-charge  electret  elements  (of  alternating  sign),  and 
the  armature  plates  are  conductors.  [  1986] 


5.3.2  Force  Analysis 

In  consideration  of  the  previous  transformation,  it  is  noted  that  the  force  per  volume 
for  the  configuration  of  Fig.  5.12  ;s  approximates  twice  that  for  the  oiiginal  M1PA- 
Fx  (Fig.  5.8).  Discretizing  the  planes  (Fig.  5.10)  lowers  the  force  by  at  most  a 
factor  of  two,  thus  making  for  a  good  approximation  of  the  MIPA 

The  force  per  volume  for  the  conductor-conductor  (c-c)  configuration  of  Fig.  5.12 
is,  assuming  air  breakdown, 

■yVol  =  eoEb2/2A 

which  is  seen  to  be  twice  that  of  the  overlapping  plate  configuration. 

The  force  per  volume  for  the  electret-conductor  (e-c)  configuration  of  Fig.  5.13  is, 
assuming  air  breakdown, 

Fx/Vol  =.  s  Eb/A. 

Each  of  the  different  types  of  interactions  presented  above  are  plotted  in  Fig.  5.14. 
There  is  obviously  a  progression  in  increasing  force/densities  as  one  moves  from 
conductor-conductor  interactions,  to  conductor-dielectric  interactions,  to  conductor- 
electret  interactions,  to  electret-electret  interactions,  to  conductor-ferroelectric 
interactions. 
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Fig.  5.14.  Comparison  of  actuation  forces  versus  grain  size  for  a  variety 
of  configurations  (note  that  these  represent  only  rough  approximations  to 
what  the  actual  forces  may  be  in  a  real  system)  [1ST86]: 

(a)  "Electret"  type  DOC  actuator  of  Fig.  5.17.  The  surface  charge 

density  is  a  =  100  nC/cm2. 

(b)  Repeated  overlapping  conductors  of  Fig.  5.2. 

(b)  Sheets  of  conductors,  as  per  Fig.  5. 1 1. 

(c)  Repeated  dielectric  "sandwiches,"  as  per  Fig.  5.4.  Assumes  a 
dielectric  material  of  £r  =  9. 

(d)  Repeated  permanently  poled  dielectric  (PZLT)  "sandwiches,”  as 
per  Fig.  5.6.  Note  one  line  is  for  the  "passive"  force  component, 
while  the  other  line  is  "active"  component.  The  PZLT  is  assumed 
to  have  a  polarization  density  of  P0  =  30,000  esu. 

(e)  Sheets  of  conductors/electrets,  as  per  Fig.  5. 1 3. 
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Note  that  data  points  for  muscle  force  densities  are  plotted,  for  muscle  lengths  of  1 
micron,  1  cm  and  1  m. 

Note  that  all  actuator  models  are  assumed  to  be  homogeneous,  in  contrast  to 
muscle,  which  consists  of  sarcomeres  in  series. 

Note  that  a  breakdown  of  air,  with  Eb  =  3X106  V/m,  is  assumed  to  be  limiting.  If  a 
vacuum  is  used,  then  all  curves  [except  (a)  and  (d:pass)]  will  shift  up  three  orders 
of  magnitude. 

Numerical  examples  are  given  below. 

Case  1  -  Actuation  system  with  conductive  plates  only. 

The  force  per  volume,  for  a  system  of  alternately  charged  conductor  sheets  (Fig. 
5.14  or  Fig.  5.12)  is  Fx/Vol  =  EoEb2/(2A). 


C 


Fig.  5.15.  Conductor-conductor  actuator  [1986]. 


To  achieve  a  force/density  of  100  psi,  or  equivalently,  an  actuator  which  is  a  one 
inch  cube  which  generates  100  lbs  of  actuation  force,  a  grain  size  of  0.1  microns 
would  be  necessary,  assuming  air  breakdown.  These  sizes  may  be  untenable. 
However,  assuming  vacuum  breakdown,  a  grain  size  of  30  microns  would  be 
necessary,  which  is  possible  to  build  with  conventional  IC  techniques.  Of  course 
other  factors,  such  as  the  smallest  radius  of  conductive  structures,  affect  the  actual 
result. 

Case  2  -  Actuator  using  a  combination  of  electrets  and  conductor  planes. 

The  force  per  volume  for  a  system  of  conductors  sheets  interspersed  with  electret 
sheets  (Fig.  5.16)  is  Fx/Vol  =  s  Eb/A. 
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Fig.  5.16.  Conductor-electret  sheet  actuation.  [1986] 

To  achieve  a  force  density  of  100  lbs  out  of  a  one  cubic  inch  actuator  would  require 
a  grain  size  of  10  microns,  assuming  a  useable  surface  charge  density  of  100 
nC/cm2  on  the  electret  and  the  breakdown  of  air.  Assuming  a  breakdown  limitation 
in  a  vacuum  (of  50x1  CfiV/m)  requires  a  grain  size  of  130  microns. 

Case  3  -  MIPA  of  electret  elements. 

For  purposes  of  comparison,  we  repeat  the  results  presented  in  the  DARPA  report 
of  15  December  1984,  for  an  actuator  system  consisting  of  discretized  planes  of 
electrets  interacting  as  shown  in  Fig.  5.17. 
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Fig.  5.17.  MIPA  consisting  of  alternating-sign  electret 
elements.  [1982] 

The  grain  size  necessary  to  achieve  100  lbs  per  cubic  inch  was  on  the  order  of  30 
microns,  which  are  sizes  within  the  realm  of  charge-trapping  materials,  such  as 

Teflon®. 

In  summary,  the  limit  of  force  per  volume  or  stress  per  length  for  actuators  will  be 
determined  by  the  breakdown  in  air  (lower  limit,  3xl06V/m)  or  the  surface 
emission  voltage  (higher  limit,  109V/m).  As  a  function  of  grain  size  (A),  these  are 
plotted  in  Figure  5.14  as  heavy  lines.  The  vertical  dotted  lines  illustrate  the  size- 
scale  of  interest  for  the  microfield  devices. 
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5.4  Fiber-based  Actuators 

5.4.1  Introduction 

Fiber-based  actuators  represent  an  alternative  which  possesses  a  number  of 
advantages  over  planar  systems.  For  example: 

1 .  A  large  interactive  area  is  achievable  between  stator  and  armature 
elements  fra-  the  attachment  of  force  producing  elements. 

2.  Lower  precision  is  required  since  the  alignment  of  fibers  is  less  critical 
than  that  of  large  planes.  Fiber  flexibility  produces  self  aligning 
properties. 

3.  The  possibility  exists  for  self  assembly.  (A  number  of  schemes  can  be 
generated  which  seem  plausible  for  the  self-induced  orientation  of  fibers 
in  stator  and  armature  structures.  These  approaches  use  the  field 
structures  on  the  fibers  themselves  to  align  bundles.  A  discussion  of  the 
approaches  is  relatively  complicated  and  will  not  be  presented  here.) 

5.4.2  Fiber-to-Fiber  Systems  (FTF) 

One  possible  representation  for  a  fiber-based  actuator  system  is  shown  very 
schematically  as  the  fiber-to-fiber  geometry  depicted  in  Fig.  5.18.  It  uses  two 
hexagonally-arrayed  sets  of  fiber  engaged  as  shown  in  Fig.  5.19.  Along  the 
fibers,  helically  arrayed  static  and  variable  fields  are  produced,  which,  if 
appropriately  controlled,  cause  a  contraction  force  to  be  generated.  Conductors  are 
used  to  produce  controlled  fields  which  are  spatially  complex  and  whose  placement 
can  be  altered  by  commutation.  The  fields  produced  by  other  structures  are 
passive.  For  example,  Fig.  5.20  shows  an  insulated  fiber  on  which  a  double  helix 
of  conductive  material  has  been  deposited.  The  potential  on  the  helices  can  be 
controlled  from  one  end  of  the  fiber  such  that  a  desired  longitudinally-rotating 
dipole  is  produced  along  the  fiber.  Note  that  multiple  helices  can  be  wrapped  on  the 
same  fiber  as  shown  in  Fig.  5.23. 


F 


(a) 


0.0.0 


(b) 


Fig.  5.18.  Multiple  filament  actuator:  Hexagonal  array;  (a)  cross-sectional  view,  (b) 
end  view  of  Filaments,  (c)  isometric  view. 
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In  all  cases,  FTF  systems  will  include  armature  and  stator  fibers  in  hexagonal 
arrays  interdigitated  to  act  as  sliding  filaments  in  a  geometry  roughly  similar  to  that 
of  natural  muscle.  Fig.  5.19a  shows  an  end-on  view  of  armature  and  stator  Fibers 
of  equal  size.  Fig.  5.19b  shows  the  case  where  the  armature  and  stator  fibers  of 
different  size. 

Helical  FTF  actuator  systems  include  stator  fibers  with  classes  (from  3  up  to  6 
pairs)  of  helically  wound  conductors  interdigitating  with  armature  fibers  which  also 
have  helically  wound  regions  (1  or  2  helices)  of  implanted  charge. 

Geometric  considerations  which  define  these  systems  are  complex,  very  interesting, 
and  in  more  comprehensive  cases,  difficult  to  visualize.  Therefore  a  number  of 
figures  are  presented  in  order  to  "walk-through"  a  number  of  simple  cases. 

Fig.  5.21  shows  a  single  helix  crossing  seven  planes  (six  segments)  perpendicular 
to  the  axis  subdividing  one  cycle  by  angles  of  1/3  n.  Fig.  5.22  shows  a  double 
helix  together  with  an  important  straight  line  marked  a,  b,  c.  Line  a,  b,  c  lies  in  the 
cylinders'  surface,  parallel  to  the  fiber's  axis  and  periodically  intersects  the  helix  at 
points  a,  b,  and  c.  Figure  5.23  shows  a  fiber  with  six  helices  which  advance  with 
the  slope  mark  "S". 

If  a  fiber  is  surrounded  by  six  other  fibers,  as  shown  in  Fig.  5.24,  then  regions  of 
close  proximity  between  the  cylinder's  surfaces  occur  as  shown  in  Fig.  5.25 
(marked  as  points  a,  b,  c,  d,  e  and  f).  Looking  down  and  inside  of  a  hexagonal 
array,  as  shown  in  Fig.  5.26,  it  can  be  seen  that  along  lines  similar  to  those 
described  in  Fig.  5.22  the  helices  of  adjacent  fibers  come  into  close  proximity.  In 
regions  of  these  "spatial  nodes",  force  can  be  produced  if  helices  are  appropriately 
charged.  In  Fig.  5.26  the  nodes  are  shown  as  points  a  and  b. 


# 
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A  more  comprehensive  case  is  shown  in  Fig.  5.27  where  adjacent  nodes  a,  b,  and  c 
(between  fibers  0  and  1)  are  staggered  with  respect  to  nodes  a^,  b^,  and  c^ 
(between  fibers  0  and  4). 

A  number  of  very  interesting  cases  can  be  generated  where  complex  nodal 
interactions  (spatial  beats)  result  between  populations  of  fibers.  These  relationships 
are  governed  by  a  selection  of  simple  parameters  which  define  each  fiber.  The 
parameters  are: 

1 .  The  initial  angle  of  the  helix,  0,  on  a  fiber  (relative  initial  angles  can  be 
observed  by  passing  a  plane  perpendicularly  through  the  population  of 
fibers  and  observing  the  intersection  of  fiber  helices  with  that  plane). 
Examples  are  shown  in  Figs.  5.28,  5.30,  and  5.31. 

2.  The  lead  rate  of  the  helices,  defined  as  "S”,  as  in  Fig.  5.23. 

3 .  The  number  of  helices  on  each  fiber. 

Fig.  5.28  shows  a  very  simple  case  where: 

1 .  All  lead  rates  "S"  are  of  the  same  magnitude  and  sign  for  Fibers  0  up  to 

6. 

2.  The  initial  angle  of  a  zero  fiber  is  0  and  the  initial  angle  for  fibers  1  to  6 
is  7t. 

In  this  case  the  1  to  6  fibers  interact  with  the  central  fiber  at  bulkheads  a,  b,  c,  d,  e, 
f,  and  g  as  shown  in  Fig.  5.29. 

Figure  5.30  shows  another  case  where: 

1 .  The  zero  fiber  and  1  through  6  fibers  have  equal  magnitude  but  opposite 
sign  lead  rates. 

2.  The  initial  angles  of  the  outer  fibers  1, 2,  3,  4,  5,  and  6  are  as  shown  in 
Fig.  5.30  (7t,  57t  /6,  rt/3,  n,  5n/6,  n/3). 

This  reversal  in  slope  causes  individual  regions  a  through  f  of  Fig.  5.25  to  pass 
close  to  each  other  in  a  roughly  parallel  fashion  rather  than  at  significant  angles. 
This  case  is  interesting,  however,  the  choice  of  parameters  does  not  permit  its 
extension  to  large  arrays  such  as  shown  in  Fig.  5.19. 

An  even  more  interesting  case  is  shown  in  Fig.  5.31  where: 

1 .  The  lead  rates  "S"  for  fibers  1  through  6  are  twice  as  large  and  opposite 
in  sign  to  the  lead  rate  of  the  zero  fiber. 

2 .  Initial  angles  are  shown  in  Fig.  5.3 1  as  either  0  or  7t. 

In  this  case,  the  six  stator  fibers  lock  together  (their  helices  cross)  at  points  marked 
as  a,  b,  c,  d,  e,  f.  Furthermore,  this  interaction  occurs  in  a  repeatable  way  such  that 
the  property  extends  across  large  arrays.  This  fact  provides  a  first  but  incomplete 
view  of  self-assembly  possibilities  for  these  systems.  That  is,  the  fiber  contains 
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intrinsic  properties  which  are  self-organizing  by  structural  design.  Also,  in  this 
case,  the  helices  of  the  zero  fiber  interact  with  surrounding  helices  along  points 
marked  g,  h,  i,  j,  k,  and  1  in  Fig.  5.31  at  various  points  along  the  fiber  as  shown  by 
Fig.  5.29. 

Discussions  of  these  systems  are  most  compelling  and  can  go  on  and  on  with 
respect  to  geometry  alone.  In  addition  to  geometrical  considerations,  numerous 
other  issues  exist.  For  example:  (1)  force  generating  capabilities,  (2)  fiber 
stability,  (3)  commutation,  (4)  wiring,  (5)  sensing,  (6)  assembly,  and  (7)  use  of 
dielectric  fluids  to  surround  fibers  and  promote  stability. 

We  intend  to  continue  studying  these  systems;  however,  significant  fabrication 
efforts  will  be  deferred  until  more  simple  actuation  systems  are  studied  and 
understood.  As  suggested  by  sponsors,  we  intend  that  a  very  small  actuator  will 
produce  forces  equivalent  to  the  weight  of  a  quarter. 
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Fig.  5.21.  Helically  patterned  field-emitting  filament. 


Fig.  5.22.  Double  helix  field-emitting  pattern  on  a  filament  showing  reference  points  a, 
b,  and  c. 
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Fig.  5.25.  A  fiber  surrounded  by  six  others  showing  regions  of  close  proximity. 


Fig.  5.26.  View  down  a  hexagonal  array  showing  regions  of  spatial  nodes. 
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Fig.  5.29.  Central  fiber  showing  bulkheads  a  through  g. 
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